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FOREWORD 


This  document  describes  work  performed  at  Philco-Ford  under  Task  2,  "Wake 
Structure  Measurements"  of  Advanced  Penetration  Problems,  ARPA  Order  888, 
monitored  by  SAMSO  under  contract  F04  701-68  -  0  0032,  in  the  period 
15  October  1967  to  15  October  1968.  Certain  subtasks  of  this  work  have, 
upon  completion,  been  issued  under  separate  cover  as  SAMSO  technical 
reports  for  the  purpose  of  documenting  information  in  detail  beyond  that 
required  contractually  in  monthly,  semiannual,  or  iinal  reports.  These 
technical  reports,  already  supplied  to  the  appropriate  SAMSO  and  ARPA 
personnel  according  to  the  pertinent  distribution  channels,  are  for 
brevity  referred  to  frequently  in  the  present  document.  They  are: 

(1)  "Turbulent  Front  Structure  ot  an  Ax isymme tr  ic 
Compressible  Wake",  Philco-Ford  Publication 
Number  UG-4259,  SAMSO  TR  68-44,  15  November  1967. 

(2)  "Mean  Flow  Measurements  in  a  Self-Preserving 
Turbulent  Plasma  Jet"  SAMSO  TR  68-166, 

February  1968. 

(3)  "Turbulent  Mean-Flow  Measurements  in  a  Two- 
Dimensional  Compressible  Wake",  Philco-Ford 
Publication  Number  UG-4431,  SAMSO  TR  68-369, 

Oct  iber  1968, 

inis  technical  report  ha.--  been  reviewed  and  is  approved. 

R.  W.  PadtieJd,  Lt.,  L’SAF 
Advanced  Penetration  Problems 
Project  Officer 


ABSTRACT 


This  report  summarizes  experimental  work  done  in  the  period  15  October  1*167 
to  15  October  1968  on  the  structure  of  compressible  turbulent  wakes  and 
turbulent  plasmas.  Detailed  point  measurements  in  the  two  dimensional 
wake,  made  very  tar  from  the  body,  have  fully  confirmed  the  predictions 
of  the  Dynamic  Equilibrium  Hypothesis  regarding  the  asymptotic  values  of 
the  velocity  and  temperature  fluctuation  and  the  statistics  of  the  inter¬ 
face.  For  flight  at  angle  of  incidence  radical  changes  in  the  wake  struc¬ 
ture  have  been  observed,  while  heat  transferred  from  the  model  to  the  flow 
was  found  to  delay  greatly  the  onset  oi  '  ans  i t ion  to  turbulence.  In  the 
plasma  jet  transition  zone  thr  remnants  oi  the  laminar  1 1 u  id  cieale  a 
turbulent  fluid  oi  highly  heterogeneous  electron  density  with  hitherto 
unexpected  statistics.  The  longitudinal  scales  oi  the  electrons  are 
larger  than  those  of  the  temperature  but  numerically  not  much  d  A  teren‘ 
from  low-speed  turbulence  scales.  The  electron  density  t  Luc  tuat  ions  alv.a 
predominate  the  temperature  fluctuations  especially  at  the  higher  tic- 
quencies.  Spectral  decay  obtains  a  (-5/3)  range  for  the  temperature  but 
this  behavior  is  not  observed  for  the  electrons. 
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SECTION  I 


SUMMARY  OF  MAJOR  ACCOMPLISHMENTS  AND  CONCLUSIONS 


1.1  OBJECTIVES 

The  objectives  of  this  work  can  be  summarized  as  follows: 

(1)  The  general  objective  is  to  predict  the  turbulence 
structure  of  the  wake  behind  a  vehicle  moving  in 
the  atmosphere  at  hypersonic  speeds. 

(2)  The  specific  objectives  are  (a)  to  measure 
experimentally  the  neutral  gas  dynamical  features 
of  turbulent  compressible  wakes  such  as  the 
mean,  intermittent,  turbulent,  and  eddy  magni¬ 
tude  distributions  of  velocity,  density,  and 
temperature  a>d  to  map  the  statistical  features 
of  these  variables,  such  as  correlation  and 
spectral  functions;  (b)  to  assess  the  importance, 
on  these  turbulence  features,  of  reallistic  re¬ 
entry  phenomena  such  as  angle  of  attack  and  heat 
and  mass  transfer;  (c)  to  measure  the  corresponding 
features  of  turbulent  plasmas  with  the  aim  of 
establishing  functional  relationships  between 
gasdynam  ica 1  and  electronic  phenomena  such  as 
electron  density,  and  gas  velocity  and  temperature; 
and  (d)  to  unity  these  results  into  laws  n't 
contradicting  other  available  experimental  results 
and  to  extend  these  scaling  laws  for  hypersonic 
wake  predictions  thr  nigh  the  basic  flight  parameters 
such  as  drag  and  altitude. 

1.2  GUIDE  TO  THE  PRESENT  REPORT 

This  report  covers  work  toward  the  c  ited  objectives  performed  from 
15  October  1967  to  15  October  1968,  and  discusses  live  major  experiments 
at  various  stages  of  completion.  In  Section  l,  the  status  of  these  exper¬ 
iments  is  given  in  tabular  torm  and  major  conclusions  drawn  from  them  are 
discussed  as  affecting  hypersonic  wake  calculations. 

Section  II  deals  with  turbulence  measurements  in  a  two-dimensional  com¬ 
pressible  wake  without  heat  transfer  or  incidence  effects  (designated 
KED) ,  and  Section  III  with  the  intermit tency  properties  of  the  same  wake. 
The  WED  mean  flow  was  measured  in  the  1966-1967  period,  analyzed  in  the 
current  period  and  reported  in  a  separate  SAMSO  technical  report 
(Reference  1).  Since  the  latter  report  has  been  printed  and  distributed 
to  SAMSO,  Aerospace,  and  ARPA,  it  is  referred  to  herein  only  in  abstract 
form  (Paragraph  1.3)  for  brevity. 

Section  IV  discusses  the  turbulence  and  intermittency  measurements  in  the 
plasma  jet  experiment  (designated  JEA) .  The  JEA  mean  flow,  measured  in 
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the  1966-1967  period,  was  also  analyzed  in  the  current  period  and  a  SAMSO 
technical  report  (Reference  2)  to  this  effect  has  been  distributed  as  above. 
It  is  abstracted  in  Paragraph  1.3. 

Section  V  discusses  the  wake  at  angle  of  attack  (designated  WEF)  and 
Section  VI  the  wake  with  heat  transfer  (WEG)  ,  and  finally,  Section  VII 
presents  recent  developments  in  techniques  and  facilities  devoted  to  this 
work  and  accounts  for  other  such  "logistics"  of  the  project. 

1.3  MAJOR  ACCOMPLISHMENTS  AND  CONCLUSIONS 

In  the  current  period,  the  WED  mean  flow  has  well  substantiated  the  earlier 
conclusions  (Reference  3)  that  the  mean  axial  velocity  and  the  ambient  wake 
temperatures  are  a  Gaussian  function  of  the  wake  radius.  The  abstract  of 
Reference  1  states:  "Excellent  correlation  of  the  lateral  (radial)  dis¬ 
tribution  of  axial  velocity  is  obtained  when  the  lateral  distance  is 
divided  by  the  transverse  wake  scale  which  is  formed  from  the  measured 
velocity  defect.  The  latter  is  found  to  decay  as  the  inverse  square  root 
of  distance  for  the  entire  range  of  about  1850  virtual  model  thicknesses 
mapped;  the  transverse  scale  thus  increases  as  the  square  root  of  distance. 
In  the  latter  half  of  this  range,  the  lateral  d  istribution  of  static  tem¬ 
perature  also  appears  to  correlate  in  the  same  lateral  coordinate,  and  the 
temperature  defect  also  decays  as  the  inverse  square  root  of  distance.  It 
is  demonstrated  that  these  results  are  accurately  predictable  from  a  basic 
similarity  analysis  beginning  with  Townsend's  measured  velocity  decay  on 
the  axis;  the  turbulent  Reynolds  number  of  13.0  agrees  closely  with 
Townsend's  12.5.  The  corresponding  Prandtl  number  found  lies  in  the  range 
from  0.65  to  0.70."  Implications  of  these  statements  are  that  the  com¬ 
pressible  wake  is  predictable,  in  the  mean,  on  the  basis  of  the  local 
(inviscid)  wake  conditions  and  independently  of  the  body  shape  (onl\  the 
drag  matters).  The  turbulent  Reynolds  number  is  now  well  settled,  and  for 
the  axisymmetric  wake  the  Prandtl  number  is,  of  course,  higher  than  0.70 
(it  is  0.83,  according  to  Reference  3). 

The  expected  pronerties  of  the  turbulent  front  (boundary)  have  been 
abstracted  in  Reference  4:  "Intermittency  measurements  have  been  performed 
within  the  first  one  hundred  virtual  diameters  of  an  axisymmei r ic  com¬ 
pressible  wake.  The  primary  measurements,  performed  with  the  hot-wire 
anemometer,  concerned  the  intermittency  factor  and  the  frequency  of  zero 
occurrences.  The  intermittent  flow  covered  most  of  the  wake  profile  but 
for  a  narrow  region  about  the  axis,  and  was  fc und  to  be  normally  distributed 
around  the  average  front  position.  The  latter,  along  with  the  extent  of 
the  front  standard  deviation,  was  found  to  agree  numerically  with  expecta¬ 
tions  based  on  low  speed  wakes  and  to  grow  as  the  1/3  power  of  axial 
distance.  A  week  periodicity  of  the  front  was  detected  at  a  wavelength 
about  nine  times  greater  than  the  longitudinal  scale  of  turbulent  velocity 
fluctuations.  This  periodicity  affects  the  turbulent  spectra."  In  the 
two-dimensional  wake  experiments  done  in  the  current  period,  these  con¬ 
clusions  were  soundly  reinforced  (Section  III)  regardless  of  the  character¬ 
istic  growth  (as  the  1/2  power  of  distance)  of  the  front  scales  for 
two-dimensional  geometries. 
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Extremely  important  conclusions  affecting  the  hypersonic  wake  turbulence 
were  drawn  in  the  current  period  from  the  WED  turbulence  studies  detailed  in 
Section  II.  The  Dynamic  Equilibrium  Hypothesis  had  already  been  demonstrated 
to  apply  during  the  axisymmetric  wake  experiments  (designated  WEB)  done  in 
196b- 196 7  (see  References  5  and  6),  and  the  same  conclusion  is  now  drawn 
from  Section  III.  Specifically: 

(1)  The  Strong  Reynolds  Analogy  holds  (i.e.,  given  the 
velocity  fluctuations  one  can  compute  the  density 
or  temperature  fluctuations). 

(2)  The  mass-flux  and  total  temperature  fluctuations 
are  uncorrelated;  the  velocity  and  temperature  are 
almost  everywhere  perfectly  anti-correlated. 

(3)  On  the  axis,  the  velocity  fluctuations  become 
asymptotically  38  percent  of  the  axis  mean  velocity; 
the  temperature  fluctuations  become  45  percent  of 
the  temperature  "deficit"  across  the  wake. 

There  have  also  been  phenomena  observed  differing  substantially  from 
previous  expectations: 

(1)  The  radial  position  of  maximum  velocity  fluctuations 
(at  a  given  distance  behind  the  body)  does  not 
correspond  to  the  maximum  shear  region  (in  contrast 
to  WEB;  see  Reference  5). 

(2)  The  radial  distribution  of  the  velocity  fluctuation 
extends  much  farther  out  from  the  axis  than  the 
comparable  low-speed  fluctuations. 

(3)  It  is  possible  to  obtain  fluctuation  intensities  which 
are  larger  in  magnitude  than  the  local  velocity  or 
density  "deficits".  The  latter,  therefore,  do  not 
indicate  the  absolute  maximum  attainable  by  the 
fluctuations  as  previously  thought. 

Although  the  angle-o f-attack  and  heat-transfer  data  are  not  fully  reduced, 
three  important  conclusions  are  drawn  from  Sections  V  and  VI: 

(1)  Even  small  angles  of  attack  at  moderate  Mach  numbers 
make  the  wake  resemble  a  shear  layer  and  create 
severe  departures  from  the  above  conclusions. 

(2)  The  boundary  (front)  structure,  measured  by  the 
local  wake  thickness,  is  unaffected  by  heat  transfer. 

(3)  Heat  ing  the  body  stabilizes  the  wake  and  greatly 
displaces  (downstream)  the  transition  region.  This 
fully  verifies  stability  predictions  (Reference  7) 
which  also  predict  destabilization  by  cooling. 
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The  conclusions  drawn  from  work  done  in  the  current  period  with  the  turbulent 
plasma  are  now  entering  the  phase  of  practical  application.  Reference  2 
summarizes  best  the  conclusions  regarding,  first,  the  mean  flow  of  the 
reacting,  highly  heated  flow:  "Measurements  of  the  temperature,  density, 
axial  velocity  and  electron  density  show  that  the  overall  gasdynamical  jet 
structure  attains  the  so-called  self-preserving  state  within  the  region 
investigated.  In  the  radial  direction  a  similarity  coordinate  based  on  the 
Howarth-Dorodnitzyn  radius  correlates  tiie  distributions  of  temperature  and 
axial  velocity,  and  it  is  shown  that  the  same  coordinate  should  also  corre¬ 
late  the  relative  concentration  of  electrons  though  not  the  absolute  electron 
density.  The  measurements  show  that  the  electron  relative  concentration 
does  not  in  fact  become  similar,  and  it  is  conjectured  that  the  axial  decay 
of  electron  density  must  be  controlled  by  recombination  reactions.  There- 
tore.  it  is  demonstrated  that  dynamical  se 1 f-equil ibrat  ion  can  occur 
independently  of  chemical  relaxation  if  the  ionization  level  is  not  too 
high.  In  Paragraph  4.2  it  is  further  shown  that  the  front  structure  is 
independent  of  the  chemistry  and,  very  significantly,  the  "electron  inter- 
mittencv"  is  indistinguishable  from  the  regular  intermittency.  This  applies 
even  to  the  "weakly  periodic"  turbulent  trout  whose  wavelength  scales  are 
the  same  as  for  the  compressible  wake. 

From  the  host  of  turbulent  plasma  phenomena  observed  and  measured  (Section  IV), 
the  following  conclusions  are  drawn: 

(1)  Two  distinct  phases  of  turbulence  exist  in  the  jet: 

A  heterogeneous  portion  apparently  associated  with 
transition  and  a  homogeneous  portion  resembling 
the  familiar  turbulence  phenomenon. 

(2)  Nearer  the  heterogeneous  region  the  gas  temperature 
spectra  have  a  very  clear  -5/3  slope  which  change 
into  a  much  faster  spectral  decay  (with  frequency) 
in  the  homogeneous  region.  The  electrons  never 
attain  the  -5/3  slope,  and  in  the  homogeneous  region 
they  decay  about  as  the  -4.b  power  of  frequency. 

(3)  In  the  homogeneous  region  the  off-axis  maxima  in 
the  e  lec  t  ron  dens  it  y  fluctuations  occur  where  the 
flow  shear  is  a  maximum. 

(4)  In  the  heterogeneous  portion,  abnormally  high  "spike" 
electron  fluctuations  make  the  rms  electron  density 
much  higher  (by  a  factor  ot  5)  than  the  mean.  These 
spikes  do  not  exist  in  the  temperature  and  are 
unrelated  to  intermit tencv ,  which  also  makes  An/n 
very  high  near  the  jet  edge. 

(5)  Agreement  of  the  gas  velocity  and  temperature  fluctu¬ 
ations,  scale  magnitudes,  etc.,  with  values  obtained 
with  non-plasma  jets  in  other  laboratories  is  generally 
very  good.  Specifically,  the  longitudinal  integral 
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scale  A-]  of  the  temperature  seems  to  bear  the  same 
relat  ion  to  the  so-callec  "transverse  scale" 

(A^  =  0.75L)  that  it  does  in  other  jets  and  wakes. 

(6)  The  longitudinal  integral  scale  of  electrons  is 
related  to  Ax  as  follows:  Ae  =  0.6  A-p. 

(7)  The  frequency- integrated ,  normalized  electron 
density  fluctuation  An/n  is  greater  everywhere 
than  the  corresponding  temperature  fluctuation 
AT/T.  At  ea- h  po^nt,  it  is  also  greater  than 
AT/T  across  the  spectrum. 

1.4  PRESENT  STATUS 

The  present  status  of  the  work  is  '’hown  on  Table  I. 
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SECTION  II 


TURBULENCE  MEASUREMENTS  IN  A  TWO-DIMENSIONAL 
COMPRESSIBLE  WAKE 


2.1  INTRODUCTION 

The  simplicity  of  the  two-dhmensional  model  ^tifa'few" 

turbulence  measurements  at  large  d^tance  ^  wakfis  carne  £rom  experiments 

years  aRo  the  bulk  o£  all  knowledge  on  A  number  ol 

with  two-dimensional  models  immersed  ™  “  “‘o[  cyUnders  placed 

investigator  reported  detaile  wor  in  rk  is  tyiat  0f  A.  A.  Townsend 

mT'ence's)  “ITttie 'present' series  of'exper iments  detailed  compressible 
turbulence  information  has  been  obtained  with  an  a*Ja^*)tr“a^‘t'rs  Irom 

(Reference  5)  where  data  were  taken ^o^about^lOOJvirtual^  ^  ^  ^ 

the  model.  All  evide.  ce  c  ion  within  that  distance,  but  the 

l^SroTt^.-n.  T  MS- 

men^of^compress  ib^le  Tur^Ience  “b^ior  with  a  ^-dimensional  model,  where 
such  distances  could  be  easily  attained. 

As  with  the  mean-flow  results,  some  important  ^differences^xtst  between^ 

turbulence  behavior  tn  a  two-dimens  tona  h  £ornier  case  as 

example,  the  rms  axis  iluctuatLon  U.el  decreas^^in 

X-l/2,  whereas  in  the  latter  it  decreases  as  X  •  J*  ^  ^  constant 

normal izi  IS  with  the  mean  Pa°per  a®^  anticipated  by  Townsend  and  already 

(i.e.,  independent  of  X)  and,  in  fac  ,  Two-dimensional  studies 

x&sx&r.  r  ;£  i.  .  *rsru«. 
set;  •<  — »« 

appl  icat  ions  . 

The  turbulence  measurements  described  ““ported  install  in 

the  two-dimensional  model  whose  wa  e  me  information  on  the  tunnel 

Kssri  rsiSE "  -rrr. 

their  method  of  acquisition. 

2.2  FLOW  FACILITY  AND  MODEL 

stagnation  pressure  of  730  mm  g  h  id  stainless  steel 

38°C.  The  model  was  a  0.004  inch  thick,  O.^inch  wid^  ^ 

ribbon  stretched  at  zero  angle  of  decrees  and  the  trailing 

leading  edge  was  beveled  to  a  half-angle  of  about  _7  degree. 
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edge  was  square  (i  e. ,  the  base  was  flat).  No  heat  was  exchanged  between 
the  model  and  flow,  i.e.,  the  experiment  was  adiabatic.  The  relevant 
Kevnolds  numbers  were  as  follows: 


(1) 


(2) 


(1)  The  Reynolds  number 

U^Cph 

Re  *  =  530 

ao  V 

00 

is  based  on  conditions  external  to  the  wake  (see 
Reference  1)  and  the  equivalent  ribbon  "base  height" 

Cph  was  -n  reality  the  wake  momentum  thickness  and 
was  measured  by  integrating  the  measured  momentum 
flux  in  the  wake. 

(2)  The  Reynolds  number 

R  .  »  .  158 

W  fr'  00 

IS  based  on  the  velocity  deficit  U*  -U(0)  between 
the  axis  and  external  velocities  U*,  and  U(0), 

respectively,  and  on  the  So-called  transverse 
scale  L. 

These  and  further  findings  of  the  mean  Mow  field  will  be  discussed  later. 

2 . 3  INSTRUMENTATION 

«obeU>rMihnCe  Tr™"CS  "'ere  PO'Co™*  «Uh  a  single  hot-w.r.  anemometer 
probe  vhrch  could  be  traversed  continuously  in  the  wake  and  could  thus 

thl  ,lkhaT?Vn  the  turbulcnce  over  distances  mm  h  smaller  than 

the  wake  thickness.  An  elaborate  description  of  the  technique,  which  .s  in 

many  respects  classical,  is  provided  in  Reference  5;  only  a  few  points  of 
general  importance  will  theretore  be  mentioned  here.  The  hot-wire  is  oper- 
ated  at  a  number  of  predetermined  electric  heating  currents  (constant 

ot  nuctnat  “"I  an<i  tl,US  providop-  at  each  Point,  a  corresponding  number 

of  fluctuating  voltage  outputs.  The  latter  are  then  inserted  as  forcing 

whic^rhe  lnt°  an  ovar-decermined  system  of  linear  algebraic  equations  in 

ft  lclent®  ot  the  unknowns  are  numerirallv  fixrj  bv  the  local 

Tnvokiirfh‘°PCrl;leS  (measurod  oarlier  and  reported  alrcadv  In' Reference  1). 
Invoking  the  constant-pressure  assumption,  the  system  is  solved  bv  the 

least-squares  method  for  the  local  fluctuations  in  the  flow  velocity, 

cient  y  -Jhis' temperatura  •  a"d  the  temperature-velocity  correlation  c^effi- 
h„^  -.1  procedure  ts  normally  performed  lirst  for  a  very  wide  frequency 

band  within  which  the  total  (integrated)  rms  is  obtained.  Afterwards  a 

each  auctua°ttral'raH°dal  a"glySlS  U  P"'°™d  »>•  the  contribution  of 

each  fluctuation  mode  is  obtained  at  eacli  point  in  the  spectrum. 

detaileV!CCeriS5iCf  °f  the  hot‘“iro-  t'"1  calibrations  performed,  and  a 
e tailed  journal  of  ns  usage  are  discussed  in  pages  bl  and  b?  oi  Reference  6. 
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2.4  GENERAL  FEATURES  OF  THE  WAKE  FLOW  FIELD 


The  two-dimensional  wake  of  the  ribbon  grew  in  a  regular  manner  behind  the 
model  for  a  distance  of  about  1850  virtual  body  diameters.  It  is  worthwhile 
to  summarize  here  some  general  features  of  this  wake  before  the  turbulence 
measurements  are  detailed: 

(1)  The  uniformity  of  the  wake  along  its  span  is 
quite  satisfactory, 

(2)  No  shock  waves  or  other  disturbances  emanating 
from  the  model  intersect  or  distort  the  wake. 

The  bow  shock  system  affects  the  wake  only  beyond 
(downstream  of)  the  region  studied. 

(3)  It  was  demonstrated  conclusively  that  the  angles 

of  attack  and  twist  were  sufficiently  close  to  zero. 

(4)  Measurements  of  transition,  fluctuation  intensity, 
spectral  density,  and  intermittency  demonstrated 
that  the  wake  was  indeed  turbulent. 

(5)  Transition  to  turbulence  occurred  within  the  first 
few  X-Stations  along  the  wake  (details  will  be 
furnished  when  the  fluctuation  measurements  are 
descr ibed)  . 

(6)  Day-to-day  reproducibility  of  the  wake  flow  field 
was  excellent. 

A  summary  is  given  below  of  the  mean-flow  results  for  the  turbulent  wake 
at  hand.  These  experimental  results  are  shown  correlated  into  simple 
algebraic  formulas  which,  at  the  same  time,  are  the  general  scaling  laws 
for  two-dimensional  compressible  wakes.  The  nomenclature  has  the  usual 
aerodynamic  connotation. 

(1)  Proper  axial  coordinate: 

X-X 

X  =  — (X  is  virtual  origin)  (3) 

VcDh 

(2)  Vel  oc ity  de feet : 

W  =  (Rt  =  13.0)  (4) 

Uoo  y  10  x  1 

(3)  Transverse  scale: 

L  =^rp-  =  0.22  VT  (5) 

4W 


9- 


(4)  Late  ral  coordinate: 


fy  ~ 

Jo  P  00 


dY  (The  liowarth-Donodn itzyn  radius) 


(5)  Proper  non-dimensional  lateral  variable: 
V  =  Y/L 

(6)  Lateral  variation  of  axial  velocity: 


u  = 

u„  -U(0) 


=  e 


-0.697]  2 


(7)  Lateral  variation  of  static  temperature: 


T-Too 

T(0)-Too 


-0.457)  1 


(8)  Axial  variation  of  temperature  defect: 

e  g  IffitS.  ,^pr(y  -  ,)  ,  0.336  H?, 


(10) 


These  laws,  which  will  be  used  frequently  to  understand  the  turbulence 
behavior  at  each  point  in  the  wake,  are  applicable  primarily  in  its  self¬ 
preserving  portion  so  that  they  are  only  approximate  in  the  lirst  20  percent 
of  the  wake  length  studied,  and  exact  thereafter. 

The  pertinent  Reynolds  numbers  which  are  relevant  to  this  wake  are  shown  in 
Paragraph  2.2.  The  quality  and  behavior  of  turbulence  can  be  prejudged  on 
the  basis  of  two  more  Reynolds  numbers,  added  to  those  indicated  above, 
lirst,  consider  the  Reynolds  number  of  turbulence: 

Au*  A 


where  Au  is  the  rms  velocity  fluctuation,  A  the  integral  scale  of  turbulence, 
and  v  the  local  kinematic  viscosity  of  the  gas.  Now,  asymptotically  Au  is 
about  0 .4(  j-U(O))  (on  the  axis)  whereas  A  is  of  the  order  of  Lhe  transverse 
scale  u  anu  v  is,  tar  in  the  wake,  of  the  order  of  the  free-stream  kinematic 
viscosity  Therefore,  Rerp  is  about  0.4  Rew  where  the  latter  is  the  wake 

Reynolds  number,  already  noted  in  Paragraph  2.2,  to  be  constant  at  about  160; 
tints,  Ko'j  -  64.  Secondly,  consider  the  dissipation  time  for  the  ver\  small 
eddies  which  Kovasznay  (Reference  9)  has  computed  in  typical  wind  tunnel 
situations.  On  the’  basis  ot  itis  computation,  a  corresponding  limit  tor  the 
Reynolds  number  ot  t lie  dissipating  eddies  has  been  derived  in  Reference  6; 
assuming  that  the  smallest  eddy  is  that  which  dissipates  over  a  converted 
distance  equal  to  its  size,  then  this  Reynolds  number  (called  Rc^)  is  40. 
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These  computations  put  bounds  on  the  expected  structure  of  turbulence 
because  they  imply  that  the  macroscale  and  the  microscale  are  about  equal. 
This  wake  is  therefore  coarse-grained  in  the  sense  that  no  eddies  smaller 
than  the  integral  scale  are  expected.  Also,  Reference  1  has  already  shown 
that  the  mean  flow  in  this  wake  is  in  all  respects  almost  exactly  analogous 
to  those  obtained  at  a  much  higher  Reynolds  number.  The  important  additional 
conclusion  is  drawn  that  the  arrangement  of  the  mean  (average)  properties, 
and  therefore  the  important  mixing  processes,  result  from  the  action  of 
the  larger  eddies  alone.  It  remains  to  see  what  the  power  spectrum  of  this 
coarse-grained  wake  is. 

2.5  EXPERIMENTAL  TECHNIQUE  AND  PROCEDURE 

2.5.1  PURPOSE  OF  THE  INSTRUMENTATION 

The  hot-wire  anemometer  was  used  to  measure  the  following  statistical 
properties  of  the  turbulent  wake:  (1)  the  root-mean-square  fluctuation  of 
the  axial  component  of  fluid  velocity  (where  "axial"  will  henceforth  refer 
to  a  direction  parallel  to  that  of  the  tunnel  stream)  resulting  from  all 
active  Fourier  components,  i.e.,  integrated  across  a  very  broad  frequency 
band;  (2)  the  root-mean-square  fluctuations  in  the  static  temperature 
similarly  integrated  across  the  spectrum;  (3)  the  spectral  "density"  of  the 
above-mentioned  velocity  fluctuation;  (A)  the  spectral  density  of  the 
temperature  fluctuation;  (5)  the  integral  scale  (correlation  length)  of 
the  velocity  fluctuations  in  the  axial  direction;  (b)  the  integral  scale 
of  temperature  fluctuations;  (7)  the  "Taylor  microscale"  of  velocity 
fluctuations;  and  (8)  the  same  microscale  of  the  temperature  fluctuations. 

It  was  furthermore  desired  to  map  these  quantities  at  each  point  in  the 
turbulent  flow. 

Some  points  should  be  stressed  regarding  these  objectives.  First,  certain 
combinations  of  the  sought-after  parameters  (such  as  the  mass-flux  and  total 
temperature  fluctuation)  are  also  directly  obtainable  by  the  hot-wire. 
Secondly,  certain  assumptions  justified  by  earlier  experiments  are  needed 
in  reducing  the  data.  For  example,  it  is  necessary  to  have  no  pressure 
fluctuations  (sound  waves)  in  order  to  extract  the  velocity  fluctuations 
from  the  data.  Third,  there  is  no  need  to  resort  to  questionable  assump¬ 
tions  (e.g.,  isotropy  or  homogene  ity)  in  order  t)  obtain  these  results. 

2.5.2  CRITERIA  FOR  RESOLUTION  AND  STATISTICAL  ACCURACY 

A  number  of  crrteria  must  be  satisfied  so  that  the  measurement  will  be  valid. 
Some  of  these  involve  a  precise  knowledge  of  the  measured  quantities  and 
thus,  can  be  checked  only  a  posteriori.  However,  with  a  precise  knowledge 
of  the  electronic  circuits  and  a  fair  estimate  of  the  turbulence  properties 
drawn  from  the  mean  measurements  (Reference  1)  and  exploratory  work,  the 
criteria  were  also  checked  a  priori. 
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a.  Space  Resolution.  In  the  direction  of  flow  (x  direction),  the 
0.00005-inch  diameter  of  the  wire  assures  excellent  resolution.*  Taking  a 
typical  wake  width  of  0.100  inch  and  estimating  a  typical  small  eddy  size 
of  0.010  inch,  one  obtains  a  resolving  power  ot  1  to  200,  applying  equally 
to  the  x-  and  the  y-direction  (i.e.,  t he  lateral  direction).  However,  in 
the  z-direction  the  hot-wire  length  is  0.010  inch,  that  is  equal  to  the 
estimated  eddy  size.  Resolution  in  this  direction  is  therefore  not  as 
good . 


b.  Time  Resolution.  Even  if  adequate  space  resolution  in  the  x-direc 
tion  is  assured,  the  question  arises  whether  the  hot-wire  and  its  circuits 
are  capable  oL  exploiting  this  resolution.  The  frequency,  f,  corresponding 
to  an  eddy  of  size  d  =  0.01  inch  =  0.025  cm  and  its  convection  velocity 
u  =  b0,000  cm/sec  past  the  probe  is  f  =  u/d  =  2,400,000  Hz.  Formally,  this 
is  the  estimated  maximum  frequency  to  which  the  turbulence  spectrum  extends 
The  hot-wire  system  must  there  Lore  operate  properly  at  these  frequencies. 
Now  the  hot-wire  employed  in  this  work  lias  a  time-constant  of  about 
0.15  msec,  that  is,  an  intrinsic  "cutoff"  at  7000  Hz;  this  is  much  lower 
than  needed.  To  bridge  the  gap  between  the  required  and  the  available 
frequency  limits,  two  supplementary  techniques  were  used,  one  based  on 
real-time  electronic  compensation  of  the  high  frequency  response,  and  the 
other  on  computer  corrections  during  the  reduction  of  data.  These  two 
schemes  have  been  described  in  Paragraph  2.b  of  Reference  5,  so  that  only 
a  brief  review  will  be  given  here. 


The  constant-current  amplifier  used  includes  an  adjustable  RC  network 
which  has  a  gain  of  unity  at  low  frequencies  and  amplifies  pre Lerent  ia  1  ly 
the  higher  frequencies.  The  transfer  function  of  this  network  is  of  the 
type 


—  ^  /l  +  03  2) 

in  '  ' 
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where  0)  is  the  circular  frequency  and  eout  and  ej_n  are  the  rms  voltage  out 
of  and  into  this  network,  respectively,  in  a  narrow  passband  about  a).  The 
circuit  time  constant  tj  is  fully  adjustable  at  the  front  panel  to  an 
accuracy  of  1  microsecond.  The  advantage  of  this  transfer  function  is 
that,  of  course,  it  is  algebraically  inverse  to  that  oL  the  hot  wire  itseli 
where 


out 


in 
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with  t  being  the  time  constant  of  the  wire  arising  from  thermal  lag.  Just 
prior  to  each  measurement,  t  is  measured  and  tj  is  set  equal  to  it.  Then 


The  wire  is  actually  not  stretched  taut  between  the  needle  supports  but 
with  a  curvature  which  is  oi  order  20  percent  of  its  length,  i.e.,  about 
0.002  inch.  Thus,  for  eddy  sizes  of  order  0.010  inch,  one  still  obtains 
good  resolution. 
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the  compensator  output  eout  is  equal  to  the  "flow"  input  a  ' in  into  the  hot 
wire  (by  multiplying  Equation  (12)  with  (  .')  as  would  be  expected  of  an 
ideal  system. 


The  upper  frequency  at  which  Equatioi  (12)  applies  ai  o  marks  the  upper 
frequency  limit  where  the  compen^a  just  indicated  u, ,  ,ps.  This  limit 

depends  on  ti.  To  find  th  iwl  lorm  of  Equation  (12),  the  a  1  fier  is 
calibrated  before  e  h  experiment.  In  the  WED  experiment  the  ca  ibiatLon 
was  least-square.c  i  "ed  by  a  trivariate  function  and  used  in  the  'EB-11 

program: 


(14) 


The  departure  of  Equation  (12)  from  the  desired  behavior  (Equation  (14))  is 
such  that  at  tr  =  0.1  msec  the  upper  compensating  frequency  is  800  kHz  and 
at  t1  =0.2  msec,  it  is  400  kHz,  which  appears  to  great iv  improve  the 
7  kHz  intrinsic  limit  of  the  wire  itself. 


Unfortunately,  there  are  other  transfer  functions  in  the  circuit  to  con¬ 
sider,  and  another  limitation  arises  from  the  "straight"  amplifier  itself. 
For  the  WED  experiment,  the  transfer  function  of  the  latter  was  calibrated 
and  curve- fitted  into: 


1  +  6.144  x  10"^F  -  7.718  x  10"bf'2  +  1.085  x  10 

(15) 
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-4.610  x  10  F 

(with  F  in  kHz)  which  indicates  a  roll-off  frequency  beginning  at  about 
80  kHz  and  an  attenuation  of  eout/ein  =  0.43  at  500  kHz.  Thus,  the 
straight  amplifier  response  is  generally  inferior  to  that  of  a  perfectly- 
matched  (t  =  tj)  wire-compensator  system.  It  is  estimated  that  the  com¬ 
plete  system  begins  degenerating  severely  around  300  kHz. 

To  extend  this  range  further  to  the  2  MHz  estimated  range  of  turbulence, 
the  "response  restoratiot"  technique  described  in  pages  13  to  29  of 
Reference  5  was  used.  This  is  based  on  the  spectrum  of  a  signal  being 
obtained  with  a  non-ideal  system  (i.e.,  a  system  with  a  non-flat  response), 
then  the  correct  spectrum  can  be  reconstructed  provided  that  the  frequency 
response  of  the  entire  system  (i.e.,  its  transfer  1  unction)  is  known 
accurately.  To  do  this,  the  frequency  response  of  the  chain  of  experimental 
equipment  is  obtained  component-by-component  so  that  the  so-called  iverall 
transfer  function"  of  the  system  is  known  at  each  point  in  the  wake.  Also, 
at  each  point,  the  spectrum  of  each  admitted  signal  is  measured,  and  a 
computer  program  is  then  used  to  restore  the  latter  spectrum  to  the  form 
it  would  have  if  detected  by  an  ideal  electromechanical  system. 

The  response-restoration  technique  is  obviously  limited  by  the  electronic 
noise  present.  In  this  experiment,  spectra  were  taken  to  1.2  MHz  so  that 
slightly  more  than  half  the  presumed  range  of  turbulence  activity  is 
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covered.  However,  it  should  be  considered  that  the  majority  nf  the  turbu¬ 
lence  signal  at  the  hot-wire  arises  at  the  lower  frequencies.  The  a  priori 
conclusion  was  that  with  the  present  system,  it  should  be  possible  to 
measure  the  rms  velocity  and  temperature  fluctuations  integrated  over  the 
entire  band  of  turbulence  activity,  but  that  perhaps  the  details  of  the 
spectrum  above  1  MHz  are  elusive.  Further  discussion  of  these  points  will 
be  made  during  the  presentation  of  experimental  results. 

c*  Statistical  Accuracy.  A  statistical  measurement  can  become  very 
inaccurate  if  the  domain  of  measurement  is  not  much  larger  than  the  size 
of  eddies  lying  within  it.  Otherwise  stated,  the  hot-wire  must  remain  at 
the  point  of  measurement  for  a  time  which  is  so  long  that  a  large  number 
of  eddies  of  all  sizes  go  by  it.  In  the  course  of  earlier  experiments 
(the  WEB),  it  was  seen  that  the  integral  scale,  for  example,  was  only 
4  times  smaller  than  the  statistical  wake  diameter.  Using  the  integral 
scale  as  a  useiul  measure  of  the  larger  eddies,  and  a  local  wake  velocity 
of  order  6  x  10  cm/sec,  the  time  of  passage  of  such  an  eddy  past  a  point 
is  estimated  typically  as  2  microseconds .  For  proper  statistics,  the 
residence  time  of  the  probe  at  a  point  should  be  long  enough  to  record 
the  passage  ot,  perhaps,  1000  such  eddies;  this  time  should  then  be 
2  milliseconds.  For  the  cont inuous- f low  tunnel  used  here,  this  time  is 
safely  exceeded  by  orders  ot  magnitude.  Even  with  the  method  of  tracersing 
the  probe  continuously  through  the  jet  the  res  idence- 1  ime  requirement  is 
fulfilled  easily,  s’nce  it  takes  several  seconds  for  the  probe  to  traverse 
a  lateral  distance  equal  to  one  correlation  length. 

2.5.3  FULFILLMENT  OF  THE  "FROZEN  TURBULENCE"  CRITERION 

Another  criterion  to  be  Lul filled  for  these  measurements  to  be  valid  is 
the  so-called  frozen  turbulence"  criterion,  i.e.,  the  Taylor  hvpothesis. 
This  states  that  in  a  locally  uniform  convection,  the  time  and  space 
coordinates  can  be  interchanged  provided  that  the  time  needed  to  change 
the  turbulence  structure  is  not  shorter  than  the1  averagmg  time  of  the 
signal.  Otherwise  stated,  the  eddies  must  remain  unchanging  as  they  are 
being  counted  by  the  probe.  Estimates  usually  made  of  the  validitv’of 
Taylor's  hypothesis  are  limited  to  changes  made  in  the  tine  (small -eddy) 
structure  of  turbulence  by  viscosity.  Kovasznay  (Reference  9)  estimates 
the  viscous  decay  time,  tv,  for  small  eddies  to  be  of  order  (vk2)-1,  where 
V  is  the  molecular  kinematic  viscosity  and  k  is  the  edd\  wake  number, 
k  =  2tt/X  (where  X  is  the  eddy  diameter).  Obviously,  it  the  "lile  distance" 
d  =  u t v  of  such  an  eddy  is  much  larger  than  its  own  size,  the  Ta>lor 
hypothesis  is  valid;  this  criterion  (see  also  Reference  6,  page  123)  can 
be  expressed  by  the  following  alternate  forms: 

utv»X,  Rey-  f»40,  X»^>;  (1S) 

Since  the  wake  conditions  will  be  close  to  ambient  over  most  o 1  the  wake, 
the  freestream  unit  Reynolds  number  Maa/v^  =  80,000  per  cm  (at  the  conditions 
used)  can  be  substituted  in  the  latter  inequality.  The  result  is  that  the 
minimum  eudy  size  which  can  be  measured  in  the  framework  of  Taylor's 
hypothesis  is  40/80,000  =  0.0005  cm  in  size.  This  rough  estimate  thus 
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shows  that  the  space-time  transformation  is  ver y  ^  ^  0’°25  C" 

eddies  estimated  to  make  up  the  fine  structure  of  the  wake  at 

2.5.4  EXPERIMENTAL  PROCEDURE 

The  hoL-wire  wa.  positioned  parallel  to  the  z  axis  (Figure  1)  at  the 

(y  -  0)  olane  *  Twenty-four  positions  along  the  x-axis 

center-span  (z  -  0)  plane.  iwe  >  r  ,  m1tion  o  l,  2  . 23. 

✓  j  " v-St it  ions")  were  chosen  and  designated  X-Statio  * 

^%aMru0ndit8edachtposeitionithetwtire0wI;  jy  -toss 

the  wake  “'the  Y  direction  and  its  output  versus  position  was  plotted  on 

function  plotter. 

c  •  9  The  hot-wire  is  connected 

^al^et^s  SSM  cons  t  ant -cur  rent 2  anemometer  set  wh  .ch  includes 

rhS  circuit  ac  amplifier  and  f^  ^Hn 

the  wire  time-constant,  ciiaracu  nsn  '  ,  a  Hewlctt- 

tne  wire  f  5  For  measurements  of  spectral  density,  a  Hewlett 

Packard  wave  analyzer  Model  310A  -s  -dj^th  a ^band^ 
tor  the  integrated  rms  measurements  a  B‘“an“n*  displa>  analyzers, 

true-rms  voltmeter  was  emp  oyed.  Asso  d  -  Uoacopu..^  |(a,>pli,ttc(1  by 
etc.,  complemented  this  system. 
a  Moselcv  Mood  7001A  function  plotter. 


2.5.5  MEAN  VOLTAGE  MEASUREMENTS 


Since  all  data  is  handled  and  utilized  in  dimensionless  ^Qvk 

(average)  wire  resistance  variation  acro^s  current  through  the  wire, 

^^that^ar^ac  curate  leas  t- square  “de  tormina  t"^^ o  t^hc^rcsulls^could  be 

made.  These  currents  are  shown  in  Table  ^  C^' 7/ irst  measured 

S  ^"hStaJrr  Sn’-.isx*:  - - 

shown  in  Figure  3. 

2.5.b  FLUCTUATING  VOLTAGE  MEASUREMENTS 

Another  set  of  fifteen  traverses  at  each  X-Sta.ion  was  made  hy^again  ^ 
setting  the  dc  current  at  the  lcvc  so  hiuh-nassed  above  1000  Hz  to 

sas  szsx  wraMgSMBw  •££& 


*No  detailed  measurements  were 
of  the  flow  at  z  =  +0.5  and  z 


taken  at  z  ft  0;  however,  the  good  quality 
=  -0.5  was  confirmed  by  preliminary  tests 
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FIGURE  1.  SCHEMATIC  OF  WED  EXPERIMENT  GEOMETRY  SHOWING  COORDINATE 
NOMENCLATURE.  MEASUREMENTS  WERE  MADE  WITHIN  THE  RECTAN¬ 
GULAR  PARALLELEPIPED  SHOWN. 


F09813  U 

FIGURE  2.  ELECTRONIC  CIRCUIT  USED  K>R  THE  WED  TURBULENCE  MEASUREMENTS 
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TABLE 

II.  AXIS 

VALUES  OF 

MEASURED 

TURBULENCE 

PROPERTIES 

X-S  ta  tion 

X 

Vu 

Pu 

At 

0 

To 

Vi 

u 

Q 

r 

mt 

r 

a  r 

0 

0.0477 

0.0624 

0.0117 

0.0357 

-1.23 

-1.09 

1 

0 

0.519 

0.0038 

0.0469 

-3.94 

2 

84 

0.131 

0.0525 

0.0962 

-0.95 

3 

168 

0.200 

0.0758 

0.0668 

0.121 

0.65 

-0.88 

4 

251 

0.152 

0.0459 

0.0436 

0. 1035 

0.38 

-0.83 

5 

335 

0.122 

0.0338 

0.031 

0.0926 

-0.05 

-0.76 

6 

419 

0.121 

0.0307 

0.0268 

0.0957 

-0.25 

-0.80 

7 

503 

0.119 

0.0367 

0.0214 

0.0933 

-0.08 

-0.84 

8 

587 

0.113 

0.0293 

0.0057 

0.0860 

-0.09 

9 

670 

0.103 

0.035 

0.0155 

0.077b 

0.14 

-0.81 

10 

754 

0.091 

0.0281 

0.0149 

0.0702 

0.03 

-0.82 

11 

838 

0.081 

0.0241 

0.0145 

0.0626 

0.01 

-0.83 

12 

922 

0.077 

0.0219 

0.017 

0.0581 

0.08 

-0.77 

13 

1006 

0.068 

0.024 

0.012 

0.0522 

0.13 

-0.74 

14 

1089 

0.068 

0.020 

0.014 

0.052 

0.09 

-0.76 

15 

1173 

0.064 

0.019 

0.015 

0.048 

0.18 

-0.76 

16 

1257 

0.059 

0.014 

0.013 

0.046 

-0.08 

-0.79 

17 

1341 

0.056 

0.016 

0.013 

0.044 

-0.08 

-0.71 

18 

1425 

0.051 

0.013 

0.0066* 

0.043 

-0.41 

-0.92 

19 

1508 

0.051 

0.013 

0.0091* 

0.041 

-0.26 

-0.80 

20 

1592 

0.045 

0.012 

0.0056* 

0.038 

-0.46 

-0.98 

21 

1676 

0.047 

0.014 

0.0049* 

0.038 

-0.20 

-1.09 

22 

1760 

0.045 

0.013 

0.0056* 

0.037 

-0.32 

-0.86 

*Uncorrec  ted 
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TABLE  III. 

HOT-WIRE  HEATING 

CURRENTS  FOR 

TURBl LLNCL  MEASUREMENTS 

Hi  pH  Ranee 

Medium  Range  (ma) 

Low 

Range  (ma) 

n  i  ju,  » •  i'a  “n'' 

7.404 

6.437 

5.320 

7.192 

6.273 

4.892 

b.988 

6.115 

4.292 

6.795 

5.825 

3.650 

6.614 

5.562 

3.176 

TABLE  IV. 

CHARACTERISTICS  OF 

HOT-WIRES  USED  IN  THE  WEI) 

EXPERIMENT 

Tre  No. 

Measuring 
Pronert  ies 

R(°C) 

a (per  "Cl 

£ 

( inch) 

i/d 

6-7/1 

Turbulence 
(all  X-Stations) 

39.10 

0.0162 

0.011 

220 

6-7/1 

Spec  tra 

X-Stat ion 

39.10 

0.00162 

0.011 

220 

7-12/1 

Spec  tra 

X-Station  3-9 

37.43 

0.00158 

0.0105 

210 

7-19/1 

Spec  tra 

X-Station  11-21 

36.37 

0.00160 

0.102 

204 
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UGl'RK  3. 


lYJ’ICAL  PROFILES  OF 


..  ...  •  HOT  VIRK  VOLT  U  K 

U  aKK  AT  CONSTANT  DISTANCE  BUI  INI)  THK  BODY 
WIRK  CURRENT  IS  CONSTANT 


U  ROSS  TVO-DIN!  NS  ION  VI 
TOR  KACII  CURVE  THK 


first  positioned  just  outside  the  wake;  its  time  constant  was  read  and 
the  compensator  dial  was  set  at  that  level.  Typical  time  constants  at 
each  point  varied  from  150  to  110  microseconds,  depending  upon  dc  current. 

The  reference  (ambient)  turbulence  level  and  electronic  noise  lev- 1 
(independent  of  Y)  were  also  recorded  at  this  point.  Then,  for  eat 
current,  the  wire  was  traversed  through  the  wake  at  a  compensation  setting 
equal  to  the  free-stream  value  corresponding  to  that  current.  During  this 
traverse  the  ac  signal  is  continuously  converted  to  dc  by  the  CL^uit  oi 
the  rms  meter  proportional  to  the  mean-square  of  the  ac  level.  Finally 
this  dc  voltage  is  plotted  versus  position.  Five  of  a  typical  set  of 
fifteen  such  "turbulence"  traces  are  shown  in  Figure  4. 

2.5.7  MEASUREMENTS  OF  THE  SPECTRA 

In  the  previous  paragraphs  it  was  shown  that  a  total  of  15  +  15  -  30  traces 
were  taken  at  each  X-Station;  for  the  total  of  24  stations  involved  this 
adds  up  to  720  traces.  These,  in  principle,  contain  all  the  information 
necessary  to  deduce  the  turbulent  fluctuations.  As  already  noted,  however, 
the  insufficient  frequency  response  of  the  system  made  it  necessary  to 
also  measure  the  spectra  of  turbulence  signals  (at  each  wire  current  and 
point  in  the  flow).  Now,  it  Well  be  recalled  (Reference  5  page  2-  • '  ^at 
the  objective  of  the  spectra,  besides  yielding  the  spectral  densitj 
tion,  is  to  provide  an  "error  ratio"  J  which  is  used  to  multiply  the 
fluctuation  intensity  discussed  in  Paragraph  2.5.6;  this  brings  “p ^ 
latter  to  the  correct  value  if  recorded  by  an  ideal  system.  This  would 
in  turn,  require  that  spectra  are  taken  at  closely  spaced  points  for  ea. 
wire  current .  Since  this  is  clearly  impractical*  the  following  was  done: 
spectra  were  taken  for  each  of  6  currents  at  each  of  10  radial  points  at 
alternate  X-Stations,  beginning  with  X-Station  1,  lor  a  total  of  7-0  spect  . 
Extensions  of  these  J's  so  obtained  to  other  currents  and  points  tould  be 
later  made  by  interpolation.  The  electric  currents  used  spanned  the  range 
of  those  used  for  the  "turbulence"  measurements  (sec  Table  III),  and  were 
3  174  4.891,  5.820,  6.432,  6.981,  and  7.401  mil  1  iampens .  hire  No.  6-7/1 
had  by  ^h is  time,  failed  and  other  wires  closely  resembling  it  were  used 
for  the  spectral  measurements.  Their  characteristics  are  shown  in 

Table  IV. 

Each  wire  was  first  positioned  outside  the  wake  where  its  time-constant 
was  measured  and  compensated  for.  The  wire  was  then  ^elT  ^ 

point  in  the  wake,  and  with  the  compensation  unchanged  from  the 
value  (for  that  particular  current),  the  spectrum  was  taken.  Amplifier  gain 
was  650,  the  passband  was  from  1  to  320  kHz,  and  the  wave  analyser  was  set 
aVl  volt  range,  normal  mode  at  -10  dB  input.  A  3-kHz  analyzer  bandwrd  h 
was  used  and  was  swept  across  the  spectrum  at  a  speed  of  about  7  to  8  ^llz/sec 
by  means  of  the  mechanical  sweep  drive.  The  signal  in  the  3  kHz  bandwidth 


For  example,  if  20  radial  points  were  chosen  for  the  24  x-stations  and 
15  currents,  one  would  require  7200  spectra  -  the  equivalent  of  about 
1000  hours  of  testing. 
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was  converted  into  dc  across  a  1000  ohm  resistor  and  Led  into  the  Y  axis 
o£  the  Moseley  plotter.  The  sweep  drive  also  produced  a  dc  voltage  pro¬ 
portional  to  frequency,  which  drove  the  other  (x)  axis  of  the  plotter.  A 
typical  spectrum  is  shown  in  Figure  5.  Since  it  is  necessary  to  know  the 
spectral  distribution  of  electronic  noise,  that  spectra  oL  the  noise  were 

also  taken. 


2.6  RESULTS 

2.6.1  RESPONSE  RESTORATION  AND  SYSTEM  PERFORMANCE 

Data  reduction  began  with  the  restoration  of  the  frequency  response  of  the 
system,  accomplished  by  feeding  the  spectra  information  (Paragraph  -.5.7) 
into  the  WEB- II  program  (SRS  No.  0073).  The  main  objective  here  was  to 
extract  the  error  ratio,  J,  at  each  point  in  the  wake. 

The  frequence  response  of  the  main  system  components  is  shown,  L  a 
typical  situation,  in  Figure  6.  The  ordinate  of  this  plot  denotes  the 
logarithm  of  "gain  squared",  that  is,  of  the  mean- square  voltage  cm  g 
of  each  compment  divided  by  the  mean-square  voltage  mpuU  The 
shows  results  at  a  single  point  in  the  wake  (X-Station  I,  *  =  *  ,/> 
id  for  one  particular  current  (7.401  ma)  .  Note,  first,  the  rapi 
deterioration  of  the  hot-wire  itself,  the  attenuation  (in  the  square  being 
about  a  million  at  600  kHz.  The  compensate  has  a  L  ime -c  ons  Lant  sett  g 
not  quite  adequate  because  of  mismatched  t  ime-c  onstants ,  and  its  gam 
never  exceeds  100,000.  The  ampl  if  ier*  a  1  so  deteriorates,  and  its  Lrequenc> 
response  together  with  that  of  the  hot-wire  and  compensator,  produce  the 
turve  labeled  "overall  transfer  faction"  (OTF)  which  is  seen  to  produce 
an  attenuation  u  about  100,000  at  1  MHz.  The  computer  program  is  thus 
set  up  to  multiply  the  wire  output  by  1/OTF  and  thus  bring  it  to  th  level 

it  would  have  ideally. 

The  hot-wire  response  is,  of  course,  controlled  by  its  time  constant. 

Figure  7  shows  tl'ie  axial  aud  lateral  variation  ot  the  hot-wire  time  con- 
stant  at  the  two  extremes  of  current  empl  >  ed:  3.17  and  7.41  ma.  c 
orientation,  it  should  be  recalled  that  the  zero-current  time  constant 
limit  for  the  type  of  wire  employed  is  about  0.09  nulliseconc  w  i in -i  ^ 
consistent  with  the  data  of  Figure  7.  More  importantly,  note  that  data  ar 
displayed  from  three  different  wires.  Except  at  the  lowest  overhea  s 
(currents)  it  is  evident  that  the  time  constant  variation  is  small 
one  wire  to  another.  it  is  also  evident  that  the  variation  o  time  con¬ 
stant  with  x  and  Y  becomes  quLte  small,  as  it  should,  at  tai  cistunns 
although  it  depends  strongly,  of  course,  on  the  heating  current  va  ut  . 

A  crucial  test  of  the  effect  of  using  different  wires  (unavoidable  because 
of  wire  breakage)  insofar  as  the  error  ratio  J  is  concerned  is  shown  in 
Figure  8.  Here  can  be  seen  the  very  large  "corrections  necessary  to  bung 


,?The  zero- frequency  gain  of  the  amplifier  ls  not  shown 
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FIGURE  5.  TYPICAL  WIRE  OUTPUT  SPECTRUM.  VERTICAL  LINES  MARK  SCALE  CHANGES 


FIGURE  b. 


FREQUENCY  RESPONSE  OF  EACH  COMPONENT,  OF  THE  COMPLETE  ^™°NIC 
EQUIPMENT  ("INSTRUMENT  TRANSFER  FUNCTION")  AND  01 
SYSTEM  ("OVERALL  TRANSEFR  FUNCTION  ) 


CURRENT: 


FIGURE  7.  VARIATION  OF  MEASURED  HOT-WIRE  TIME  CONSTANT  WITH  AX  I  *L  BU11U 

( ABO\M  )  \ND  LATERAL  (BELOW)  DISTANCE 


FIGURE  8.  VARIATION  OF  ERROR  RATIO  .1  ALONG  WAKE  AXIS.  DASHED  LINE 
ALTERNATE  SMOOTHING  (FOR  I  7.401)  WHICH  WOULD  ELIM IN ATI 
OBSERVED  ANOMALIES  BEYOND  X- STATION  17. 


the  data  to  the  Level  produced  by  a  hypothetical  ideal  system.  The  J  is 
as  high  as  15  on  the  axis  very  near  the  body  and  at  the  higher  currents. 

By  X-Station  6,  J  has  attained  asymptotically  a  constant  value  bracketed 
between  4.7  (at  high  currents)  and  2.8  at  the  lowest  current.  The  curves 
marked  "smoothed"  are  an  empirical  attempt  to  reduce  the  data  scatter  (in 
the  WEB- IV  and  WEB-VII  outputs)  by  reducing  the  scatter  in  the  J.  As  will 
be  seen  later,  the  break  in  the  high-current  curve  occuring  around 
X-Station  16  is  significant.  Equally  significant  is  the  rather  gross  depar¬ 
ture,  at  the  lowest  current,  between  the  actual  and  the  smoothed  J's  past 
the  same  point.  In  the  meantime,  it  is  clear  that  there  is  no  obvious 
break  in  the  curves  around  X-Station  10  where  the  change  from  wire  7-12/1  to 
7-19/1  occured.  As  is  seen  from  Table  IV,  these  two  wires  are  structurally 
similar;  the  conclusion  is  that,  with  this  similarity,  uniform  J's  are 
obtainable . 

The  lateral  variation  of  J  for  I  =  7.401  ma  is  shown  in  Figure  9  for 
various  X-Stations.  It  is  noteworthy  that  off-axis  J  reaches  values  even 
larger  than  the  axis  values.  Finally,  Figure  10  shows  the  variation  of  J 
with  current  for  two  lateral  positions  in  the  wake.  In  all  these  figures, 
the  disconcerting  fact  is  noted  that  J  does  not  attain  the  low  values 
(about  1.5)  in  the  free-stream  which  it  attained  in  the  axisymmetric  wake 
(WEB)  experiment  (see  Reference  5,  page  19).  At  this  juncture,  the  most 
likely  explanation  is  that  there  is  much  more  high-frequency  content  in  the 
freestream  turbulence  of  this  experiment  than  in  the  WEB  experiment.* 
However,  no  turbulence  spectra  were  taken  in  the  free-stream. 

With  the  large  J's  shown  on  Figures  8  and  10  and  the  deteriorated  frequence 
response  of  the  system  shown  in  Figure  6,  it  is  clear  that  the  output 
spectrum  of  the  measuring  system  will  be  severely  distorted  at  each  current 
and  flow  point.  The  printed  output  of  WEB-II  lists  both  these  actual 
(distorted)  spectra  and  the  "ideal"  spectra;  these  are  of  course  voltage 
output  spectra,  at  this  point  still  unresolved  moda1 iy  (i.e.,  not  indicating 
the  velocity  and  density  spectral  densities  separately).  A  comparison 
of  actual  and  ideal  spectra  are  typically  shown  in  Figure  11.  They  differ 
greatly,  as  they  should.  The  differences  lie  not  only  at  the  high  fre¬ 
quencies  where  the  actual  spectra  are  greatly  attenuated,  but  also  in  shape. 
The  latter  very  important  point  will  be  discussed  further  in  Paragraph  2.6,4, 

A  final  important  point  regarding  the  output  of  WEB-II  concerns  the  two 
alternate  methods  of  measuring  the  (integrated)  flunctuations  intensity. 

This  is  formally  measured  by  the  wide-band  Ballantine  voltmeter  (see 
Paragraph  2.5.4),  but  is  also  independently  available  from  the  spectral 
measurements  (Paragraph  2.5.7)  since  the-  total  rms  intensiu  is,  alter 
all,  the  integral  of  the  spectral  density.  The  latter  integral  is  listed 
in  WEB-II  as  VTVM2  (in  units  ul  volts  squared)  while  a  total  rms  measurement 


* 

The  WEB  experiment  was  done  at  P0  --  508  mm  Hg  where  the  tunnel  side-wall 
boundary- layer  capable  of  radiating  onto  the  wake  was  laminar.  In  WED 
(PQ  =  730  mm  Hg)  the  layer  was  turbulent,  with  possible  enhancement  ot  the 
stream  turbulence  by  radiation. 
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FIGURE  10.  VARIATION  OF  J  WITH  IIKATINt  CURRENT  FOR  TWO  LATERAL  POSITIONS 


-30- 


(1)  (ABB.  UNITS ) 


RESTORED 


FIGURE  11. 


EFFECT  OF  COMPUTER-AIDED  RESPONSE  RESTORATION  ON  THE  MODALLY 
UNRESOLVED  TURBULENCE  SPECTRA  (WED  EXPERIMENT) 


made  in  conjunction  with  the  spectra  is  labeled  E02.  Representative 
values  for  the  ratio  E02/VTVM2  are  shown  in  Figures  12  and  13.  This  ratio 
should  be  unity  for  perfect  correspondence,  but  it  is  in  tact  almost 
everywhere  constant  near  1.5  (i.e.,  a  difference  of  about  20  percent  in  the 
root-mean-square)  meaning  that  the  directly  measured  output  is  alwa>s  higher 
than  that  reconstructed  by  integrating  the  spectrum.  Because  ot  the 
multitude  of  operations  necessary  to  the  latter,  the  20  percent  difference 
is  actually  much  better  than  expected  and  is  in  tact  an  endorsement  of 
the  accuracy  of  the  measurement. 

2.6.2  THE  WEB- III  PROGRAM 


The  next  step  in  the  data  reduction  process  combines  the  J  values  trom 
the  WEB- I I  output  and  the  turbulence  data  (Paragraph  2.5.6)  into  the 
uEB-III  program  (SRS  Program  0083)  from  w’hich  the  integrated  rms  wire 
voltage  is  obtained  at  each  point  oi  the  wake,  and  for  each  heat.ng  current, 
tree  from  frequency-response  limitations  of  the  system. 

Turbulence  information  was  obtained  from  22  equidistant  points  spaced 
0.004  inch  apart,  beginning  at  the  axis  and  extending  laterally  to 
'  =  0.084  inch  The  total  points  involved,  therefore  were  (22  Y's)  x 
(15  currents)  x  (23  X-Stations)  =  7590,  since  the  24th  X-Station,  that 
numbered  23,  was  not  used.  At  the  WEB- II  output  the  number  of  J's  available 
were  (10  Y's)  x  (6  currents)  x  (12  X-Stations)  =  720,  since  the  spectra  were 
taken  sparsely  to  save  time.  The  required  J's  were  then  obtained  trom 
those  available  by  interpolation.  In  the  latter  process  a  certain  amount 
of  smoothing  of  the  data  was  performed  in  order  to  decrease  experimental 
scatter  in  the  final  outcome.  An  indication  ot  this  type  of  smoothing 
is  already  shown  in  Figure  8.  Certain  parts  of  this  process  are  thought 
to  have  caused  artiticial  perturbations  ot  the  final  data,  a  point  which 
will  be  discussed  further. 
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FIGURE  12.  CONSISTENCY  CHECK  WITHIN  THE  WEB- II  PROGRAM:  AXIAL  VARIATION 


XSTATION 
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FIGURE  13.  CONSISTENCY  CHECK  WITHIN  THE  WEB- II  PROGRAM:  LATERAL  VARIATION 


2.6.3  EXPERIMENTAL  RESULTS  IN  THE  FLUCTUATION  INTENSITY 

The  WEB-IV  program  gives  the  intensity  of  turbulent  fluctuations  at  each  of 
the  22  lateral  positions  and  the  23  X-Stotions  (X-Station  0  through  22). 

This  paragraph  discusses  the  nomenclature,  method  of  analyzing  the  data, 
and  the  results  themselves.  The  latter  deal  only  with  the  so-called  modal 
analysis,  i.e.,  the  fluctuation  intensity;  the  spectra,  as  well  as  the  scaling 
rules  by  which  the  results  are  applicable  to  high  speeds,  are  discussed  in 
subsequent  paragraphs. 

In  the  text  and  figures  found  below,  the  following  nomenclature  is  used: 

(1)  A  u/u  is  the  local  root-mean-square  fluctuation  in 
the  axial  velocity  u,  normalized  by  the  local  mean 
(average)  axial  velocity,  and  integrated  over  all 
frequencies .* 

(2)  AP/p  is  analogous  to  (1),  above,  for  the  local  gas 
density. 

(3)  AT/T  is  analogous  to  (1)  for  the  local  static  temper¬ 
ature.  Note  AT/T  -  Ap/p  for  constant  pressure. 

(4)  APu/pu  and  ATq/T0  are  analogous  to  (1)  for  the  local 
mass-flux  pu  and  total  temperature  T0. 


(5)  v 


mt 


(Apu/pu)  (AT  /T  ) 

_  1  ‘ _ Q  Q 

(Apu/pu)  (AT  /T  ) 
o  o 


:i7) 


is  the  cross-correlation  coefficient,  at  the  same 

point,  between  pu  and  TQ  fluctuation  (the  barred 

numerator  obviously  is  the  averaged  product  of  the 

instantaneous  counterparts  of  Apu/pu  and  Al  /T  )  . 

^  r  o  o 

(6)  r  -£Wu)  (AT/T) 

K  '  OT  (Au/u)  (AT/T) 


(18) 


is  the  cross-correlation  coefficient,  at  the  same 
point,  between  u  and  T  fluctuations.  Note  that  although 


AP  _  A  T  (rms  sense) 
P  T 


(19) 


ins  tantaneously 

=  .  (to) 

\  P  /INST.  \T  / 


INST. 


(20) 


*  In  this  work,  and  in  Reference  5,  A  is  used  as  a  root-mean-square  operator 
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so  that  the  following  relation  exists  between  the 
velocity-density  and  the  velocity- temperature 
correlations : 


^&F/P)(^u/u)  =  (AT/T)(Au/u)  _  (21) 

(&P/P)(£ u/u)  "  £t/T)(Au/u)  "  "r0T 

(7)  In  the  WEB-IV  output  the  following  alternate  desig- 


nations  are  used: 

M  *  Apu/pu 

(22) 

T  ^  AT0/T0 

(23) 

RMT  =  r 

mt 

(24) 

0  =  Sigma  = 

Ap/p 

=  AT/T 

(25) 

T  =  Tau  = 

Au/u 

(26) 

RSTAU  =  r 

6  r 

(27) 

U  PRIME  h  Au/j 

S- 

u(o)) 

(28) 

RHO  PRIME  ^  Ap/  | 

(p„- 

P(o)) 

(29) 

where  in  the  latter  two  cases,  00  and  (0)  refer  to 
conditions  outside  of  and  on  the  axis  of  the  wake, 
respectively,  and  where  Au  and  A  P retain  their  rms 
meanings . 


(8)  Eta,  or  T)}  is  the  normalized  radial  distance  (see 
Equation  (7)). 


The  equation  relating  the  fluctuation  intensities  to  the  local  rms  voltage 
output  e  of  the  hot-wire  can  be  written  in  the  following  alternate  ways: 

2  2  2  2  2 

e  -e  M  +e_x  ■  2  MT  e  e  r  (30) 

m  T  m  T  mT  v  ' 


2  2  2  2  2 

e  -  e  T  +  e  no  +  2  t a  ea  c  r  (31) 

T  °  t  v  ' 

with  the  symbols  as  given  above  and  with  e  ,  e j  and  e a  being  the  appropriate 
"sensitivity  coefficient"  to  the  various  fluctuation  modes.  These  coefficients 
are  in  principle  functions  only  of  the  local  mean  flow  properties,  which  were 
measured  earlier  (see  Paragraph  2.4  and  Reference  1)  and  of  the  wire  character¬ 
istics,  provided  by  the  wire  calibration  and  inputed  into  WEB-IV.  Equation  (30) 
is  exact,  but  Equation  (31)  is  approximate  in  that  contributions  to  e^  due  to 


-36- 


sound  waves  are  assumed  ni.*  Either  equation  may  be  used,  depending  on  which 
of  the  unknowns  are  to  be  obtained.  Separate  computations  are  performed  in 
WEB-IV,  one  for  each  of  Equations  (30)  and  (31). 


In  the  present  program,  Kovasrzay's  original  suggestion  (Reference  9)  is 
followed  reducing  these  two  formulas  into  convenient  second-degree  equations 
with  unknown  coefficients.  For  example,  Equation  (31)  is  written: 


which  is  of  the  second  degree  in  the  known  and  controllable  parameter  (^t • 
Since  at  each  heating  current  setting  a  pair  of  number  (e-^/e,,)  and  e  /a ~o 
are  measured,  three  such  settings  result  in  a  system  of  three  F.quations  (31) 
in  the  three  unknown  t,  n  and  r^.  The  accuracy  is  much  improved  in  the 
present  system:  fifteen  settings  are  used  and  the  computer  least-squares 
fits  a  second-degree  curve  to  the  fifteen  points. 


a.  Sensitivity  Coefficients.  Figure  14  shows  typical  variation  of  the 
hot-wire  sensitivity  coefficients  at  a  point  in  the  wake.  The  true  indication 
of  wire  sensitivity  is  its  so-called  "overheat"  A'  ^vhich  is  roughly  pro¬ 
portional  to  the  percentage  increase  of  wire  resistance  over  its  equilib¬ 
rium  unheated  value).  A  range  of  A'  from  perhaps  0.05  to  0.5  is  desired 
usually  in  order  to  maintain  measurement  accuracy,  the  wire  melts  around 

A^  =  0.6.  In  Figure  14,  A^  varies  from  0.02  to  0.40.  In  the  example 
shown,  note  that  Ay  =  0.335  (corresponding  to  I  =  6.7  ma)  the  hot-wire 
sensitivity  to  veloci^"  fluctuations  is  nil,  and  the  wire  is  therefore 
sensitive  to  density  fluctuations  onlv. 

b.  Mass -Flux-Total -Temperature  Cor re  1 at  ion.  Figure  15  shows  the  axial 
(along  Y  =  0)  and  Figure  16  a  typical  lateral  variation  of  rmt,  the  cross¬ 
correlation  coefficient  between  pu  and  TQ.  No  correlation  was  found  in  either 
direction,  as  also  was  noticed  in  the  axisymmetr  ic  wake  (Reference  5,  page  5^.). 


c.  Mass-Flux  and  Total  Temperature  Fluxtuation.  Typical  lateral  vari¬ 
ations  of  \pu/p u  and  \To/T0  appear  in  Figures  17  and  18,  and  their  axial 
variation  (along  Y  =  0)  in  Figure  19.  The  \To/T0  fluctuations  are  always 
much  smaller  than  ,\pu/pu  by  a  factor  varying  from  3  to  6. 

The  small  size  of  ,\T0/T0  is  caused  by  the  fact  that  the  flow  is  adiabatic, 
and  it  should  be  zero  if  the  effective  Prandtl  number  of  the  wake  was  unity. 
As  will  be  seen  later,  with  very  small  AT0  the  so-called  "Strong  Reynolds 
Analogy"  holds  and  enables  one  to  predict  the  temperature  (or  density) 
fluctuations  from  known  velocity  fluctuations. 


*Formally,  inclusion  of  sound  waves  introduces  three  more  terms  in  Equation  (31) 
and  three  more  unknowns:  The  pressure  fluctuation  and  the  pressure-velocity 
and  pressure-density  correlations.  It  should  be  also  noted  that  both 
Equations  (30)  and  (31)  are  linearized,  that  is,  applicable  for  small  magni¬ 
tudes  of  the  unknowns,  M,  T,  a  and  T. 
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X  STATION  15 


FIGURE  14.  EFFECT  OF  HEATING  CURRENT  (AT  A  FIXED  I’OINT  IN  THE  WAKE)  ON 
VARIOUS  SENSITIVITY  COEFFICIENTS  OF  THE  HOT  WIRE 
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FIGURE  15.  VARIATION  OF  THE  MASS -FLUX  AND  TOTAL- TEMPERATURE  CROSS - 
CORRELATION  COEFFICIENT  ALONG  THE  WAKE  AXIS 
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FIGURE  16.  VARIATION  OF  THE  MASS -FLUX  AND  TOTAL- TEMPERATURE  CROSS - 

CORRELATION  COEFFICIENT  IN  THE  LATERAL  DIRECTION  (FIXED  X) 
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Laterally  both  types  01  tluctuations  show  the  off-axis  na  •.  inun  characteristics 
ol  symmetric  free  turbulent  flows;  this  will  be  discussed  later. 

\’e  1  jc  it;>  F  1  -ic  tujti  uvs  .  The  as  Lai  \ariat  Lon  of  velocity  tluctuations 
\u/u  is  shown  in  Figure  20.  As  with  the  axial  behavior  ot  \T0/T0  and  \pu/pu, 
the  tluctuations  in  velocity  initially  increase  until  they  reach  a  value  of 
about  0.0b  at  a  distance  o'  pprox imately  200  diameters  trom  the  body  and 
thereatter  decrease.  In  the  end  of  the  range  probed  (beyond  X-Station  17) 
the  tluctuations  become  quite  anomalous.  It  is  conjectured  that  this  rather 
abrupt  change  is  due  to  incorrect  smoothing  of  the  J  lactors,  that  is,  to 
the  artitically  decreasing  J  shown  beyond  that  distance  on  Figure  8.  To  test 
this  conjecture  the  following  test  was  performed:  the  J  values  (lor  each 
heating  current)  at  X-Station  21,  Y  =  0  were  taken  as  inputted  into  WEB- IV 
program  but  prior  to  the  change  noted  around  X-Station  17.  Actually,  the 
J  values  at  X-Station  15,  Y  =  0  were  used.  The  results  were  as  follows: 

Using  J  from  X-Station  21: 

T  =  \u/u  =  0.0049 

a  =  \p/p  =  0.038 

Using  J  from  X-Station  15: 

T  =  0.0113 

o  =  0.039 

Two  very  important  conclusions  are  drawn  from  this  computition:  (1)  if  the 
J’s  used  in  the  "anomalous"  region  (X-Station  17  to  22)  are  made  identical 
to  those  prior  to  that  region,  then  the  \u  fluctuation  levels  come  back  on 
what  seems  to  be  the  level  properly  expected;  (2)  the  density  fluctuation  \/C 
are  almost  completely  unaffected  by  changing  J.  In  changing  J,  \u  changes 
by  about  250  percent,  whereas  \p  by  about  3  percent  only.  The  explanation 
of  this  behavior  is  a^  follows:  In  the  computational  method  used  (which) 
finds  the  fluctuations  from  the  shape  parameters  of  a  curve  fitting  experi¬ 
mental  points),  the  velocity  fluctuations  result  from  the  slope  of  a  curve 
(specifically,  the  e^/e£  versus  (eT/ea)  curve  of  Equation  (32)),  whereas 
the  temperature  fluctuations  result  from  the  intercept  of  this  curve  with 
the  e-j-/ea  axis.  It  turns  out  that  the  slope  of  the  curve  is  very  sensitive 
to  J,  whereas  the  intercept  is  not.  At  this  juncture,  the  anomaly  in  J  is 
studied  so  that  its  cause  can  be  located. 

Lateral  distributions  of  the  velocity  fluctuations  are  typically  shown  on 
Figure  21.  This  lateral  profile  follows  the  general  trend  of  the  fluctua¬ 
tions  in  \ou/pu  and  \T0/T0:  a  region  of  relatively  low  turbulence  on  the 
axis,  with  a  considerable  increase  off  axis  and  the  natural  decline  toward 
zero  at  the  wake  edge.  Note  that  the  maximum  in  the  fluctuation  intensity 
occurs  at  7)  definitely  higher  than  1.0.  This  maximum  is  heretofore  considered 
to  mark  the  lateral  location  in  the  wake  where  turbulence  production,  shear, 
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Ls  also  Lhe  greatest.  However,  tile  mean-flow  experiments  described  in  Refer¬ 
ence  1  indicate  that  the  region  of  maximum  shear  lies  nearer  the  point  rj  -  0.85. 
At  this  stage  the  conclusion  is,  therefore,  that  turbulence  production  alone 
cannot  account  tor  the  maximum  in  \u/u. 

e.  Density  (Temperature)  Fluctuations.  The  axis  value  of  the  fluctuation 
LP/p  in  the  density  is  shown  in  Figure  21,  while  Figure  22  shows  a  typical 
lateral  distribution.  Because  of  the  insensitivity  of  this  type  of  fluctuation 
to  J,  there  is  no  abnormal  scatter  in  this  data.  It  is  significant  that  this 
fluctuation  is  higher  everywhere  than  uu/u,  and  this  difference  will  later  be 
explained  on  the  basis  of  the  "Strong  Reynolds  Analogy"  for  an  adiabatic  body. 

2.6.4  THF  FLUCTUATION  SPECTRA 

In  addition  to  the  resolution  ol  the  fluctuations  into  the  different  modes 
(e.g.,  velocity  and  temperature),  it  is  necessary  to  carry  out  a  spectral 
analysis  lor  each  mode  separately.  This  analysis  is  doiu  by  the  WEB-VII  com¬ 
puter  program,  which  supplies  additional  information  for  each  point  in  the 
wake  such  as  the  longitudinal  integral  scales  of  velocity  and  temperature. 

These  computations  are  not  completed  as  yet.  In  what  follows,  the  analysis  of 
spectra  is  performed  (from  the  WEB-IV  output)  on  the  modally  unresolved  hot¬ 
wire  output. 

Figure  23  compares  restored  spectra  of  Lhe  signal  at  three  different  positions 
along  the  wake.  The  absolute  level  of  these  spectra  is  irrelevant  at  tbit 
point,  but  the  either  features  shown  are  significant  beyond  tin  context  of  the 
present  experiment. 

First,  observe  tin  intense  peak  (e*n  rgv  concent  rat  ion)  for  tin  pi  tra  at 
168  and  670  diameters  behind  the  model.  This  peak  is  quite  unnatural  for 
longitudinal  (i.e.,  one-dimensional)  spectra  ol  turbulence,  which  should  in 
fact, resemble  the  spectrum  shown  at  1676  diameters.  Secondly,  note  that  the 
peak  frequency  decreases  with  increasing  X;  finally,  the  peak  intensity  also 
decreases  with  increasing  X.  The  interpretation  ol  this  phenomenon  is  as 
follows;  The-  turbulent  front  is  arranged,  geometrically,  intei  a  weakly* 
periodic  manner  so  that  "waves"  of  a  distinct  wavelength  are  formed  by  it. 

To  a  stationary  observer  a  fixed  distance  way  from  the  axis  this  "woviiuss" 
appears  as  periodic  in  time.  If  the  front  was  strongly  periodic  then  the 
periodicity  in  time  (i.e.,  a  peak  in  the  spectrum)  would  be  sharp  and  intense. 

■\  lact,  by  "weakly"  periodic  it  is  meant  that  the  stationary  variable  repre¬ 
senting  the  front  has  a  random  spectrum  (probably  much  like  that  of  turbulence) 
with  a  superimposed  peak.  This  is  called  the  " in termi t tency"  spectrum.  In 
reality  the  hot-wire  responds  to  a  superposition  of  this  spectrum  and  the 
actual  turbulence  spectrum,  and  it  Ls  apparent  that  the  peak  persists  in  this 
superposed  phenomenon.  This  interpretation  thus  states  simply  that  the  peak 
is  due  to  the  intermi ttency ,  and  thus  implies  that  it  should  be  absent  on  Lhe 
wake  axis  and  at  its  edge,  and  prominent  in  between  where  the  intermi ttency  is 
most  intense.  Figure  24  shows  five  typically  restored  spectra  taken  at  the 
same  X  position  and  wire  current.  On  the  axis  and  at  the  edge  of  Liu  wake 
(Y  =  0  and  0.048,  respectively)  the  peak  is  hardly  discernible,  but  indeed 
becomes  very  prominent  in  the  intervening  region. 


-47- 


8- 


SPECTRAL  INTENSITY 
(ARB  UNITS) 


FIGURE  23. 


FREQUENCY  (KHz)  B0866 U 

EVOLUTION  OF  THE  HOT-WIRE  MEAN -SQUARE  VOLTAGE  OUTPUT  SPECTRUM 
ALONG  THE  WAKE  AXIS 


9- 


40 


X  STATION  7 
I  -  7.40  MA 

Y  *«  LATERAL  DISTANCE  (INCH) 


SPECTRAL 
INTENSITY 
(ARB.  UNITS) 


f (KHz) 


FIGURE  24.  TYPICAL  LATERAL  EVOLUTION  OF  THE  HOT-WIRE  MEAN-SQUARE  VOLTAGE 
OUTPUT  SPECTRUM 


-50- 


This  interpretation  was  first  thought  likely  in  the  axi.symmetr  ic  wake  work 
(see  References  4  and  5)  and  was  in  fact  used  to  compute  the  front  wavelength, 

Ap,  in  that  experiment,  giving  Ap  9L  —  2d ;  here  L  is  the  transverse  scale  and 
d  the  approximate  wake  diameter.  The  interpretation  can  be  further  strength¬ 
ened  if  a  closely  similar  Ap  is  found  for  the  present  experiment.  Using  a 
[ieak  frequency  of  360  kHz  and  a  stream  velocity  of  600  meters/sec,  a  wave¬ 
length  Ap  =  1/6  cm  =  0.17  cm  is  computed.  Furthermore,  d,  as  defined  above, 
is  twice  the  distance  from  the  axis  to  the  root -mean-square  position  of  the 
turbulent  front,  represented  by  Ypp  and  shown  versus  x  on  figure  47,  page  91 
oi  Reference  6.  At  X-Station  7,  Y2  =  0.4  x  10' *(  inch)-  and  YpT  =  2  x  10"  . 
Therefore,  d  -  2YpT  =  0.04  inch  =  0.1  cm,  so  that  Ap/d  =  0.17/0.1  =  1.7, 
whereas  the  axisymmetr ic  wake  \p/d  =  2  as  noted  above.  Although  exact 
identity  is  not  expected,  it  is  clear  that  the  ratio  ot  the  front  wave¬ 
length  to  the  wake  diameter  is  nearly  independent  of  geometry,  as  expected. 

2.7  SCALING  LAWS  AND  SIMILARITY 

It  remains  to  show  how  the  present  turbulence  results  car  be  used  to  derive 
or  support  scaling  laws  which  are  useful  in  predicting  hypersonic  wake 
properties.  In  the  absence  of  a  theory  of  turbulence,  the  present  results 
can  be  used  to  test  and  support  the  most,  if  not  the  only,  reasonable  empirical 
scaling  laws  grouped  under  the  title  "Dynamic  Equilibrium  Hypothesis  (DEH) . 

The  DEH  emerged  from  the  original  work  of  Townsend  (Reference  8);  its 
credentials  include  a  satisfactory  demonstration  of  its  applicability  to  com¬ 
pressible  flows  by  Morkovin  (Reference  10)  and  equally  satisfactory  prediction 
of  the  results  of  the  only  detailed  compressible  wake  experiment  to  date 
(Reference  5).  In  the  absence  of  heat  transfer  or  sound  waves  these  scaling  law 
detailed  in  the  latter  reference,  give  the  following  prediction  for  the  self¬ 
preserving  compressible  wake: 


P  \u 

P(0)  Au(0) 

ten 

(23) 

P(0)3  /  u(0)^ 

I2  (  Ap  ) 

2 

=  fOl) 

(34) 

p  3  \  u  ) 

\  \p(0)  / 

=  a-  1)  m2  ~ 

T  u 

(35) 

r  =  -1 

a  r 

(36) 

M _  = 

Uqo  -u(0) 

a(0>  «0)  j 

1/2 

(37) 

\P  (0)  _ 

Poo-P(0) 

fg(0)  ^ 

r  oo 

r 

(38) 

-51- 


-  ^°(0)  u op  -  u(0)  _  (1-r)  (1-w) 

Poo-  P(0)  Au(0)  "  S 


(39) 


U  the  mean  flow  properties  are  know  as  a  function  of  x,  the  above  formulas 
can  describe  fully  the  fluctuation  magnitudes  and  their  variation  from 
point  lo  point  in  the  compressible  wake.  Since  the  physical  quantities 
represented  by  the  above  symbols  are  now  known  (for  example  r*  w,  p,  etc.,  are 
known  from  the  mean  measurements  of  Reference  1  and  the  turbulence  properties 
arc  presented  m  Paragraph  2.6)  the  objective  of  this  paragraph  will  be  to  see 
if  these  formulas  indeed  hold.  Note  that  S  has  been  determined  to  be  0.84 
toi  the  present  wake,  while  the  DEH  gives  the  value  of  g(0)  as  0.14. 

Program  WEB-VII  gives  the  measured  turbulence  properties  in  a  form  exactly 
as  shown  above.  Since  the  WEB-VII  output  is  not  available  yet,  this  compari¬ 
son  will  not  be  made  as  complete  as  desired.  Thus,  in  Figures  25  and  26  the 
lateral  distributions  shown  are  in  the  form  Au/(up-u(0))  and  AP/(C-p(0)) 

thfcomnr^  tHhC,  fr  s“ssested  by  '-nations  (33)  and  (34)  which  account  (or 
Dcnwn  p  blllty’  the  ratLonall2ati°n  here  is  that  at  these  far  distances 
r( 0 ) / P oo  is  very  close  to  unity  anyway.  Note,  first  that  the  scaling  with  T) 

1Si?0jrrUnta  about  X-Station  11 .  Beyond  that  station  the  data  collapse  on 
well-defined  curves;  lateral  similarity  is  thus  obviously  attained  both  for 

MMVCl^ltyi  and  tlrC’  density  fluctuation,  as  predicted  by  Equations  (33)  and 
U^).  The  shape  of  the  correlating  curves  will  be  discussed  later. 


The  all-important  constant  g(0) 
and  Ap(0)/(poo  -  p(0))  versus  x. 
normalized  fluctuations,  drawn 
(37)  and  (38)  in  the  limit  Poo  = 


can  be  checked  by  plotting  Au(0)/(uqc  -  u(0)) 
In  Figure  27  the  asymptotic  values  of  these 
in  solid  lines,  are  obtained  from  Equations 
P(0)  and  tD=  0,  that  is,  in  the  far  wake: 
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The  data  show  remarkably  good  agreement  with  the  value  of  g(0)  used  above, 
the  slight  deviation  implying  a  g(0)  a  little  greater  tha  0.14.  In  this 
figure  three  distinct  regions  are  distinguished:  Up  to  about  X-Stat ion  3 
the  transition  process  is  still  continuing;  from  there  to  about  X-Stat ion  8 
for  the  velocity  and  X-Stat ion  12  for  the  density  the  fluctuations  deviate 
considerably  from  the  DEH-pred ic ted  asymptotic  value,  which  sets  in  beyond 
those  stations.  The  disparity  in  the  approach  to  equilibrium  between  the  two 
fluctuation  modes  is,  of  course,  predictable  because  of  the  generally  slower 
decay  of  temperature  in  high-speed  wakes.  What  is  interesting  is  that  from 
work  done  earlier  the  initial  (relaxation)  values  of  Au(0) /(u«  -  u(0))  were 
predicted  to  be  higher  than  the  asymptotic  value,  as  indeed  verified  in 
Figure  3/  but  the  density  fluctuations,  predicted  to  be  small  until  they  reach 
the  asymptote,  actually  overshoot  the  latter,  reaching  a  level  of  about  0.75 
at  about  X-Stat ion  6 . 
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FLUCTUATIONS  (FAR  WAKE).  BELOW:  NEAR  (RELAX  INC )  WAKE. 
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I' rom  .i  cl  imens  ional  standpoint,  the  velocity  "scale"  at  any  distance  behind  the 
body  is,  ot  course,  the  local  difference  uCXJ  -  u(0)  ;  that  is  why  liquation  (3  7) 
asymto  t  ica  lly  connects  \u  and  u^,  -  u(0)  .  The  question  arises  whether 
Vi  can  ever  exceed  u & -  u(0) .  In  Figure  28,  the  maximum  values  of  the  fluctu¬ 
ations  Au/  U  00  u(0)  and  \p/(Poo  -  P  (0))  are  plotted  as  obtained  from  the 

ott-axis  maxima  of  Figures  25  and  26.  Asymptotically  these  reach  values  of 
U .  70  and  0.80,  respectively.  In  the  relaxing  region,  however,  \u/(uoo-  u(0)) 
teaches  a  value  equal  to  one,  and  near  X-Station  6  the  density  fluctuations 
reach  125  percent  of  the  mean  density  difterence  across  the  wake.  the  precise 
physical  rationale  allowing  such  large  fluctuations  is  not  clear. 

Sine  1  g(0)  controls  L  it  h  the  velocity  and  t  lie’  density.  fluctuations  ,1  ar  1  rom 
the  body,  a  simple*  check  of  the  DE11  is  suggested  I  y  Equation  (3Lt).  Figure  29 
P  lots  the  1  e  1  t  -  hand  -  s  idc  of  this  equation  as  measured  dircctlv,  jnd  indepen¬ 
dently  the  r  ight  -hanci-s  iue  (by  a  sol  id  curve)  with  the  w,  r,  and  S  supp!  ied 
tom  the  mean  measurements.  As  expected,  the  agreement  is  poor  in  the  relax¬ 
ing  portion  of  tlie  wake,  but  be  y one!  900  diameters  the  data  are  close  t  i  the 
expected  behavior.  In  tact  the  data  beyond  1400  diameters,  it  corrected  tot 
the  difficulties  met  in  computing  J,  continue  in  u  i"  good  agreement  with  "• 
the  "theoretical"  prediction. 

The  t  e  1  a  t  i  \  c  magnitude  o  t  density  ( t  cm  pc*  rat  lire)  and  velocity  t  I  uc  t  ua  t  i  c  ns  v  i  h 
respect  to  the  stationary  hot-wire,  arc  given  by  Kquifi  on  (59).  Figure  50 
shows  hat  the  expectations  are  met  (especially  with  data  corrected  for  d  ar 
discussed  above)  in  the  sel  t -preserving  wake.  in  the  latter,  in  tact,  the 

Strong  Reynolds  Analogy  is  cl  iselv  Loll  >wed  he  the  data  in  the*  lateral  coor- 

dniate  except  v  rv  near  the  wake  edge.  Here  the  temperature  thutuati  mis 
aie  uncommonly  1  a  rge ,  as  was  a  I  so  observed  in  the  axisvmnc  trie  wake*  (figure  35, 
pa.Lc'  78,  Re  tore  nec*  5).  This  is  explained  as  ooc  urr  i..g  due  to  intermit  toiuy, 
implving  that  t  lit  temperature  "jt  nps"  across  the-  turbulent  trout  are  much 
lai.ei  than  tin*  correspond  ing  velocity  jumps. 

I:  siiuild  Iw  added  that  the  Dynamic  Equilibrium  Hypothesis  predicts  that  lor 

ll  abat  ic  wakes  \1L,  -  0  and  r(j  =  -1  .  Doth  '  I  these  prc’cl  i.c  t  i  mis,  and  espe- 

ciail'  tin  latter,  ae  Lully  borne  out  in  the  present  experiment,  as  shown 
h\  f  igu  i  c  s  15  and  .50  .  1  he  result  is  that,  oiu  e  tin  ved  h  it\  t  1  uc  t  uat  i  on  in 
m  adiabatic  wake  is  known,  the  temperature  (density)  fluctuations  can  lie 
predicted.  in  i  ii  c  t  ,  the  ratio  o  t  density  t  >  velocity  I  1  uc  L  ua  t  i  o  n  s  is  a  simple 
t  unc  t  ion  at  the  Mach  number.  This  is  well  borne  out  w  Figures  51  and  32. 

A  final  point  concerns  the  rather  high  velocity  liuctuaticms  me  mured  in  the 
wake.  Figure  jj  compares  these  measurements  with  data  taken  hv  Townsend 
(Reference  8)  at  low  speeds.  It  should  be  noted  that  Townsend's  nean- 
vclocity  measurements  are  correlated  in  tin  following  empirical  formula 
(Reference  8,  page  162): 

2 

~  _  u  on  -  u  -14.4a"" 

II  —  1  —  »  ‘ 

u  ,  -  u( 0) 
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FIGURE  28.  AXIAL  DEV ELOl’MENT  01-  HIK  MAXIMUM  ( - \T  FaCH  X)  OF  THE  VELOCITY  AND 
DENSITY  F'  DCTl  AT  IONS  SCALED  TO  THE  LOCAL  L'AKI  DEFICITS 
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FIGURE  30. 


LATERAL  VARIATIONS  OF  THE  VELOCITY- TEMPERATURE  CROSS -CORRELAT I 
COEFFICIENT 
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FIGURE  H.  IF  ST  OF  TUI  STRONG  REYNOLDS  ANALOGY  ALONG  THE  WAKE  AX  I 


FIGURE  32. 


TEST  OF  THE  STRONG  REYNOLDS 
AT  CONSTANT  X 


kNALOGY  IN  THE  LATERAL  DIRECTION 
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B 143  51' 

VELOCITY,  SHEAR,  AND  VELOCITY  FLUCTUATION  PROFILES  (ABOVE) 
COMPARED  WITH  SIMILAR  PROPERTIES  OF  LOW- SPEED  WARES  (BELOW) 


where 


a  = 


(\-\0)  d 


1/2 


(43) 


and  d  is  the.  cylinder  diameter.  This  lateral  coordinate  a  can  be 
to  the  r)  used  herein  since  Cor  cylinder  the  urag  coefficient  CD  - 
since,  in  the  present  experiment 


relati  d 
1  and 


Y 


~  Y 


a)  = 


(X-Xo)^ 


H/2 


Since  in  Townsend's  case  one  obtains 

|cDh|1/2  =[h]'/2  -VS" 


The  result  becomes 

1/2 


L  = 


and 


a  = 


(x  -  x  )  d 

_ o 

4 . 36 


Y 


(-4) 


Y 


Y 


[(x-xo)d]TT2  4.56L  4.561.  ->.56 

With  this  transformation  from  a  to  T),  Equation  (-*-)  becomes 


(->5) 


~  -0.69 1)1  (  1  +  0.05 !»»' 

u  =  e 


(46) 


which  is  in  good  agreement  with  the  mean-velocity  distribution  lound  tn 
the  r'resent  work: 


u  -  e 


-o.eg1")2 


In  Figure  33  the  u  curves  are  thus  identical,  as  at  c  the  coi  i  espond  ins; 
distribution  of  shear  du/3Y.  The  significance  of  the  Utter  is  tha  l  s 

},,  the  present  wake.  The  net  result  is  that  the  fluctuation  thickness  of 
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the  present  wake  is  considerably  thicker  than  at  low  speeds.  No  explanation 
of  this  phenomenon  can  be  offered  at  this  time. 

2.8  CONCLUSIONS 

The  following  conclusions  can  be  drawn  from  the  present  work  on  two-dimensional 
wake  turbulence: 

(1)  General  Conclusion:  The  data  obtained  support 
1  irmly  the  empirical  body  of  predictions  known 
as  the  Dynamic  Equilibrium  Hypothesis. 

(2)  The  axis  values  of  temperature  and  longitudinal 
velocity  fluctuations  reach  asymptotically 

Llie  sane  percentage  of  the  mean  axis  values  as 
encountered  in  previous  experiments,  including 
those  at  low  speeds. 

(3)  The  wake  total  temperature  fluctuations  are 
relatively  insignificant  when  no  heat  is  ex¬ 
changed  between  the  flow  and  the  body. 

(-0  The  local  velocity  in  the  wake  is  almost  per- 
tictly  an t i -e or re  1  a  Led  with  the  local  tempera¬ 
ture;  then  is  no  correlation  between  total 
temperature  and  mass-flux  fluctuations. 

(5)  lor  adiabatic  wakes,  the  temperature  (density) 
fluctuations  can  be  obtained,  by  the  Strong 
Reynolds  Analogy,  from  a  knowledge  of  the  dis¬ 
tribution  of  mean  Mach  number  and  velocity 
flue  tuat ions . 

(6)  Near  the  wake  edge  the  Reynolds  Analogy  fails; 
it  is  conjectured  that  the  front  corrugations  of 
the  wake  can  be  thought  of  as  "jumps"  in  the  mean 
properties,  especially  the  temperature. 

(/)  Similarity  in  the  distribution  of  fluctuation 

occurs  long  after  the  mean  properties  themselves 
become  self-similar. 

(8)  By  comparing  almost  all  measured  mean  and  turbulent 
properties  with  their  low-speed  counterparts,  it  is 
shown  that  the  proper  scaling  length  is  not  the 
body  diameter  but  rather  the  appropriate  momentum 
thickness  of  the  wake. 

(9)  A  notable  exception  shows  in  the  lateral  distribu¬ 
tion  of  velocity  fluctuations;  the  maximum  of  these 
occurs  farther  from  the  axis  and  is  higher  than  its 
counterpart  at  low  speeds. 
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SECTION  III 


INTERM ITTENCY  OF  THE  TWO-DIMENSIONAL  COMPRESSIBLE  WAKE 


3.1  PURPOSE 

In  addition  to  the  mean-flow  investigation  of  the  two-dimensional  wake 
(Reterence  1)  and  ol  its  turbulence  Lield  (Section  II),  its  intermit  lent-} 
characteristics  have  also  been  measured  and  will  be  reported  in  this 
Section. 

The  intermittent  flow  plays  an  extremely  important  role  in  the  characteris¬ 
tics  of  the  hypersonic  wake  as  viewed  by  a  ground  observer.  Lu  Lhe  v  1 r k 
with  the  ax  isymme  trie  compressible  wake  (Reference  4)  the  intermi  t  tenc  \ 

measurements  gave  the  rms  location  3  >i  the  turbulent  front  and  its 

standard  deviation  n.  These  quantities  were  found  to  bear  numerical  1> 

the  same  relationship  to  the  transverse  scale  L  as  in  the  two-d  irnens  L  >na  1 

wake  of  Townsend's  (see  Reference  5  ,  page  42  ).  This  important  finding, 
implying  that  neither  the  geometry  of  the  body  nor  compressibility  affects 
Y/L  and  0/L,  would  so  greatly  simplify  hypersonic  wake  calculations  that 
it  merited  additional  study;  hence  the  interm  1 1  tenc  y  ncasui  ements  with  the 
present  two-uimens ional  wake  were  made. 

3.2  TECHNIQUE  AND  PROCEDURE 

The  technique  necessary  for  this  type  of  measurement  is  dtset  ibed  in  great 
detail  in  Reference  5.  Briefly,  the  hot-wire  anemometer  is  held  aL  various 
lateral  positions  Y  and  at  each  X-Station.  The  heating  current  of  the 
element  is  increased  to  the  maximum  sensitivity  possible  and  it  ac  signal 
is  processed  through  the  intermittency  meter  built  especially  tor  such 
work  at  Aeronutron ic .  Through  a  series  of  rectifying  and  clipping  act.  h  s 
thi«  instrument  distinguishes  between  the  time  the  wire  is  immersed  in  dc 
turbulence  and  the  time  it  is  not;  its  output  in  a  fraction  of  one  volt  ac 
then  gives  the  "intermittency  factor"  0<y<l  which  is  the  fraction  of 
turbulence  immersion  time  to  the  total  time.  In  a  fully  turbulent  tluiu 
<-y  =  whereas  whenever  the  sensor  is  so  far  from  the  a_xis  that  it  never 
encounters  turbulence,  7=0.  If  7(7,  x)  is  known  then  3  (x)  and  a  (x) 
can  be  calculated  by  statistical  methods  (Reference  11,  page  10). 

In  addition  to  7,  the  average  wavelength  of  the  front  is  important.  A 
measurement  of  the  so-called  null  frequencies,  N  sec'  ,  can  give  this  sort 
of  wave  length  in  the  form  of  statistical  scale  of  the  front  it  the  latter 
is  truly  random.  For  the  ax isymme tr ic  wake  and  also  from  work  with  the 
jet  (Section  IV),  it  appears  that  although  the  front  is  truly  random  in 
amplitude  it  may  not  be  truly  random  in  time  and  the  existence,  if  any, 
of  its  wavelength  (rather  than  the  purely  statistical  scale)  became 
another  objectives  of  this  work. 

Hot-wire  No.  8-2/1  was  used  Lor  measuring  the  two-dimensional  wake  intermit¬ 
tency.  Data  were  taken,  beginning  at  Station  2,  for  al 1 X-Stat ions  used 
during  the  mean  and  turbulence  measurements. 
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3.3  RESULTS 


Figure  34  shows  the  front  lateral  radius  Y  and  its  standard  deviation  a, 
normalized  with  the  transverse_scale  L,  plotted  along  the  wake.  It  is 
quite  evident  that  the  ratios  Y/L  and  Q/L  are  constant  within  the  scatter 
.(Note  that  L  increases  by  a  factor  ot  3  in  the  range  200<x<l800).  For 
\/L  this  constant  value  is  about  1.96,  which  is  only  7  percent  lower  than 
the  value  expected  on  the  basis  of  the  known  behavior  of  wake  flows  (see 
Table  III,  Reference  <-*)  •  The  findings  of  the  WEB  experiment  are  therefore 
further  substantiated. 


The  value  of  Ct/L  however,  shown  on  Figure  34  to  be  about  0.38,  is  con¬ 
siderably  lower  than  expected,  by  a  Lactor  of  2  or  more.  This  is  shown  in 
another  manner  on  Figure  35,  where  the  lateral  distribution  of  "Y  is  shown  for 
all  XSTATIONS .  The  data  are  scattered,  so  much  so  that  a  distinction 
between  near-wake  and  far-wake"  behavior  is  made.  What  is  important, 
however,  is  that  the  region  of  intermittent  activity  is  much  shorter  (and 
the  extent  of  fully  turbulent  fluid  much  greater)  than  Townsend  found  at 
lov  speeds,  this  means  that  O  should  be  also  much  smaller,  as  Figure  35 
already  showed.  _In  the  mean-time,  the  distribution  of  intermittent 
activity  around  Y_is  fairly  symmetric;  this  is  also  illustrated  on 
Figure  3b.  That  \  grows  as  x^/-  is  also  well  illustrated  in  Figure  37, 
where  the  Y  corresponds  to  the  maximum  of  the  probability  distribution 
d7/dY. 

To  explain  these  findings,  two  alternatives  seem  likely:  (1)  there  is  a 
substantial  difference  between  this  wake  and  Townsend;  (2)  there  were 
inadequacies  and  errors  in  the  experimental  method.  At  this  stage,  the 
second  alternative  appears  more  likely.  Since  the  corrections  applied  b) 
computer  to  compensate  tor  the  hot-wire  trequency  response  are  applicable 
to  the  rms  sense  only,  there  is  much  oi  the  intermit tency  behavior  which 
the  wire  does  not  see.  The  turbulence  caused  by  the  iront  its,  in  tact, 
proport  ional  to  "Y(l-Y)  and  thus  is  maximum  at  7  =  1/2,  where  Y  is  located. 
Near  the  wake  axis  (Y  -  1)  and  at  the  edge  ( Y  =  0)  this  turbulence  is 
small  and  the  wire  does  not  sense  that  it  is  crossing  into  and  out  of  the 
front.  Under  these  circumstances  the’  location  ot  Y  will  be  given  correctlv 
but  a  will  be  greatly  underestimated,  as  indeed  round. 

The  c  one  1  u.  ion  then,  is  as  follow’s:  the  position  Y  ot  the  turbulent  front 
is  found  with  fair  accuracy  for  WED,  but  the  standard  dev iat ion  i^  under¬ 
estimated  due  to  inadequacies  in  the  technique.  Improveme nt s  in  the  hot¬ 
wire  ane  ometer,  presently  under  way,  will  enable  a  more  accurate  repetition 
of  this  measurement. 
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FIGURE  35. 


DISTRIBUTION  OF  INTERMITTENT Y  FACTOR  IN  WED  COMPARED  WITH 
LOW- SPEED  RESULTS 
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1-1(1  R!  (  *  LATKRAL  DISTRIBUTION  OF  Tin:  FRONT  LOCATION  PROBABILITY  (ABOVF)  . 
BFLOW,  POINTS  DFNOTF  Till  MOST  PROBABLF  FRONT  LOCATION. 


SECTION  IV 


TURBULENCE  BEHAVIOR  OF  THE  PLASMA  JET 


4.L  CRITIQUE  OF  CURRENT  TURBULENT  PLASMA  RESEARCH 

The  work  described  so  tar  deals  exclusively  with  oertect-gas  wakes  devoid 
ot  chemical  effects,  and  from  which  geometric,  mechanical,  and  to  some 
extent  thermal  features  ot  wake  turbulence  can  be  extracted.  A  sec  ind  area 
of  effort,  described  in  this  Section,  deals  with  tree  turbulence  in  wli  ic  i 
the  wake  feature  is  discarded  in  favor  of  producing  hyperthermal  and 
electronic  phenomena.  The  turbulent  plasma  jet,  such  as  employed  here,  is 
a  unique  and  inexpensive  tool  tor  producing  such  phenomena. 

It  is  fair  to  say  that,  in  view  ot  the  present  shortage  >t  and  demand  tor- 
information  on  electron  behavior  in  turbulent  plasmas,-''  accurate  and  we 
resolved  data  obtained  from  laboratory  plasmas  are  significant  regardless 
oL  their  direct  applicability  in  the  hypersonic  wake  problem.  Thus,  the 
data  already  obtained  under  the  present  contract  (Reference  b)  as  well  as 
those  of  Guthart  (Reference  12),  Granatstein  (Reference  U),  and  Johns  >n 
(Reference  14)  have  for  the  tirst  time  yielded  information  on  the  electiin 
density  fluctuation  intensity  and  spectra,  the  correlation  scales,  etc. 
Phenomena  peculiar  to  electron  benavior  (contrasted  t>  ordinary  turbulence) 
have  also  been  observed  and  are  reported  in  Paragraph  4.0  below. 


For  these  data  to  be  useful  tor  the  present  purpose,  they  should  Ih  com¬ 
piled  into  scaling  laws.  Presently  these  scaling  laws  appear  very  d  i  i  t  lcul t 
to  obtain  theoretically.  Even  in  the  case  where  they  are  obtainable 
experimentally  in  a  jet,  their  simple  f  ranspos  it  ion  to  the  case  ot  tu  cum 
is  not  straightforward.  It  the  consensus  u  Retercnccs  o,  1  - ,  I  h  and  1* 
was,  for  example,  that  the  turbulence  scale  ot  electrons  is  a  constant 
fraction  of  the  local  jet  radius,  there  is  still  no  assurance  that  the 
same  relation  applies  to  the  wake.  Even  tor  totally  inert  chemistry,  as 
another  example,  the  laws  of  axial  electron  density  decay  assume  dmerent 
forms  depending  on  compressibility  (high  or  1  n,-  speeds)  geometry  (two- 
dimensional  or  ax  isymir.etr  ic)  and  type  ot  flaw  (jet  oi  wake  >  . 

It  appears  that  the  major  usefulness  it  the  jet  experiments  is  n  t ind  an 
empirical  connection  between  the  electron  turbulence  behavior  ai  1  the 
turbulence  behavior  ol  the  overall  gas.  For  instance,  it  a  unique  relaii 
between  the  electron  and  the  temperature  integral  scales  is  tound,  then 
the  application  to  the  wake  would  only  require  knowledge  ot  the  temperature 
scales  in  the  wake  -  knowledge  which  is  now  he  ng  accumulated  under  this 
contract  (Sections  l  through  HI).  It  ls  therefore  necessary  to  measure 


*A  substantial  body  oL  knowledge  exists  on  the  so-called  "plasma  turbulence 
which,  however,  involves  phenomena  on  the  atomic  scale  with  electric  unc 
magnetic  fields  and  is  of  no  consequence  hen  . 
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the  fluid-dynamical  turbulence  in  the  jet,  as  well  as  the  electron  turbulence; 

nrnhl  laS-ber  th°  PaSt  yGar  Under  this  “’"tract.  The  next  crucial 

1  1cm  is  to  tine  the  one  or  more  variables  on  which  this  empirical  connec- 

or' \t?eZ  ,  HaS  already  been  done  ^nd,  in  tact,  the  problem  solved) 

or  the  two  extremes  m  chemistry:  equilibrium  flow  (References  15  and  lb) 

c"ntrolMn  Jr  (Reference  13)‘  For  thG  Present,  the  choice  of  parameter 
controlling  the  connection  between  gasdynamic  and  electron  fluctuations  is 

nv.  lear ,  as  the  discussion  of  turbulence  results  will  show. 

There  is  one  phenomenon  which  does  not  depend  overtly  on  the  chemistry: 

werJ  til  tUrUlenCevintermittency-  Intermittency  data  in  the  plasma  jet 
rc  taken  during  the  previous  contract  period  (Reference  b,  pages  17b-180). 

n  the  present  contract  period  these  data  were  reduced  and  interpreted  and 

U  be  dlSCUS  flrSt  Pri°r  L°  the  turbulence  results,  in  the  following 
paragraphs.  ^ 

4.2  INTERMITTENCY  IN  THE  PLASMA  JET 

The  intermittency  properties  of  wake  flows  have  already  been  obtained  and 
iscussed  in  Reference  4  for  the  axisymmetric  wake  and  in  Section  III  of 
this  report  for  the  two-dimensional  wake.  Scaling  laws  have  been  given 
(Reference  4,  page  47)  by  which  the  geometric  structure  of  the  hvpers  mic 
wa  e  front  can  be  predicted.  The  question  now  arises:  will  these  scaling 
laws  appear  identical  to  a  "temperature"  and  an  "electron"  observer?  In 
other  words,  does  the  front  surface,  which  is  known  to  abruptly  separate 

gases  of  dissimilar  temperature  also  separate  the  electron-rich  from  the 
e lectron- poor  gas? 

4.2.1  PROCEDURE 

the  present  experiment  the  front  properties  were  measured  with  the  hot¬ 
wire  anemometer  and  the  Langmuir  probe.  The  signal  from  the  hot-wire 
amplifier  (or  the  ac  amplifier  connected  across  the  Uad  oi  the  Langmuir 
probe)  was  led  to  an  intermittency  cirenit;  this  circuit,  described  in 
eference  4  is  not  unlike  similar  devices  used  by  Bradbury  (Reference  17) 
and  others,  and  supplied  the  intermittency  factor  as  a  friction  of  one 
volt  dc  The  factor  7  was  recorded  at  each  of  about  20  radial  positions 

f  at  each  x  The  rms  front  position  V  and  its  standard  deviation  a  were 
then  obtained  from  the  usual  formulas: 


00 

=  fx  — 

J  dY 


dY 


(47) 


-(/  (v-'b2  ^d(v-Y>| 

o 


1/2 


(48) 


In  these  formulas  computational  convenience  requires  that  7  is  expressed 
points“with°the  ^computcr^  ^  ^  — mcnln, 


-72- 


In  addition  to  Y  and  a  the  autocorre  Lat  Lon  stale  of  the  front  is  of  some 
interest.  If  the  front  behaves  like  a  stationary  random  variable  it  is 
possible  to  find  this  scale  from  a  measurement  of  the  so-called  null- 
crossings  (or  null  frequencies)  oi  the  tront,  that  is  by  the  frequency 
N0  '  N(Y)  at  which  a  sensor,  located  at  Y  =  Y,  enters  the  front.  To  tins 
end  the  distribution  of  null  frequencies  N(Y)  with  rudiut-  wut  recorded  by 
electronically  counting  the  rectangular  pulses  at  the  Schmid t- t r igge r  output 
of  the  interm  it  tencyr  circuit.  These  distributions,  which  of  course  peak  at 
or  neav  Y,  were  also  computer- f  itted  by  polynomials.  The  microscale  of 
the  front  is  given  by 


o 


(49) 


where  u  (Y)  is  the  jet  velocity  at  Y.  An  alternate  anJ  possibly  more 
meaningful  front  "wavelength"  is  simply 


o 


(50) 


4.2.2  RESULTS  AND  DISCUSSION 


Figure  38  shows  a  typical  comparison  of  the  distributions,  with  radius, 
of  the  intermit  tency  factor,  the  temperature,  and  t  lie  electron  density. 

As  expected,  the  fluid  is  strongly  intermittent  at  distances  from  the  axis 
at  which  the  mean  temperatures  is  very  close  to  the  ambient.  Also  irn 
Figure  38,  observe  that  the  corresponding  relation  between  electron 
density  and  the  intermit tency  factor  is  analogous  but  highly  inconsistent 
from  one  x  to  another:  at  a  fixed  value  ot  7  =  0.5,  for  example,  rf  varies 
considerably  in  the  x-range  investigated. 

Important  questions  to  be  settled  in  analyzing  the  in  term  i  ttency  data 
include  the  randomness  of  the  interface,  the  coincidence  between  electron- 
and  temperature- interface  and  the  comparison  ot  the  front  geometry  and 
behavior  to  that  of  low-speed  jets.  The  growth  of  the  tront  ^radius"  Y 
is  shown  dimensionally  in  Figure  39  and  non-d  imens iona 1 1 y  as  Y/L,  together 
with  the  standard  deviation  (j/L,  in  Figure  40.  There  is  little  doubt  that 
these  two  quantities  grow  linearly  and  scale  with  L,  and  that  they  cer¬ 
tainly  are  independent  of  the  method  ot  measurement.  With  Y  and  a  so 
computed  it  is  possible  to  make  a  test  ot  the  randomness  ot  the  tront  such 
as  is  shown  in  Figure  41;  we  see  that  the  Gaussian  cl  is  t  r  ibut  ion_o  L  the 
intcrmittency  tactor  about  the  front  position,  that  is  with  (Y-Y)/o,  again 
obtains  in  the  tashion  characteristic  ot  stationary  random  variables.  In 
the  same  figure  further  evidence  ot  randomness  is  supplied  by  the  distribu¬ 
tion  ot  the  normalized  crossing  frequencies  N/Nniax  about  the  radial  location 
Ymax  =  T(Nmax).  The  scatter  is  appreciable,  but  enough  data  points  are 
available  to  illustrate  the  Gaussian  behavior  of  N/Nmax  characteristic  ot 
randomness.  As  before,  there  is  no  systematic  difference  between  dat.i 
taken  with  the  hot-wire  and  data  taken  with  the  Langmuir  probe . 
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FIGURE  *9.  MOST  PROBABLE  TURBULENT  FRONT  POSITION  IN  THE  JET 
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FIGURh  40.  AXIAL  VARIATION  Oh  Tilt;  ITRBULKM  FRONT  POSITION  AND  ITS  STANDARD 
DKVIATION  IN  THE  J  KT ,  NORMALIZED  WITH  THE  TRANSVERSE  SCALE 


Th«  null  frequencies  N0>  defined  as  the  crossing  frequencies  measured  at 
the  position  of  the  interface,  were  found  to  be  very  ,  1  ,se  to  the  maximum 
frequency  Nmax  measured  (at  each  x)  at  Ymax;  thus,  N,  a  N(Yraax>  =  N(?) 

max  as  expected.  Figure  42  shows  the  variation  of  N„  with  x;  avain 
he  frequencies  measured  with  the  hot-wire  and  the  Langmuir  probe  are  about 
the  same  at  each  x.  The  decrease  of  N0  shown  in  .is  figure  immediately 
suggests  the  lengthening  of  the  front  scale  XF  (or  Ap)  expec  ted  as  the 
urbulent  flow  widens  with  distance.  However,  the  computation  of  these 
scales  from  Equations  (49)  and  (50)  cannot  be  performed  with  the  velocity 
shown  so  tar ;  these  latter  are  mean  velocities,  whereas  what  is  of 
course  needed  is  the  velocity  ol  the  front  itself  or  at  least  the  mean 
velocity  within  the  front  -  the  so-called  eddy  velocity. 

Measurement  of  eddy  velocities  within  "corrugations"  in  the  jet  boundary 
re  made  with  the  double  Langmuir  probe  (two-u.  c  correlation  probe) 
mentioned  previously  and  pictured  on  page  126  of  Reference  6.  Eddy 
velocities  with  such  a  probe  can  best  be  measured  i rom  the  cross-correlation 
unction  ol  the  probe  signals;  in  the  present  work  the  mire  primitive 
method  was  adopted  of  obtaining  a  finite  number  of  dual-trace  oscillograms, 
the  turbulence  pattern  as  it  was  converted  by  the  wires.  From  these 
samples  a  rudimentary  statistical  distribution  of  the  local  eddy  velocity 
was  constructed  at  each  point  from  which  the  mean  velocity  (and' its  fluctua- 

tion)  could  be  extracted.  The  (mean)  eddy  velocity  at  Y  was  then  used  in 
place  of  u  in  Equations  (49)  and  (50). 

Reference  4  describes  experiments  done  on  the  intermit tenev  features  oi 
compressibie  wakes  and  advises  the  use  of  the  "front  wavelength"  Ay  (rather 
than  AF)  as  a  more  significant  measure  of  the  front  geometry  In  do  mg  so, 

TnliTll  lSf  n  tr°m  ob£frvations  (Reference  5)  of  the  turbulence  spectra 

but  r  Hp  0n^P°Sltl0n  Y)  th3t  the  front  ma>'  not  bo  lruiy  random , 

but  rather  weakly  periodic  at  a  wave  length  Ap.  The  asymptotic  tv  m 

.f/L.  7  “/N°L  01  th ’s  wave length  for  the  axisymmetric  wake  is  drawn  in  Figure  . 
in  which  the  jet  front  wavelengths,  obtained  per  the  previous  paragraph,  are 
also  presented.  The  agreement  between  these  two  independent  experiments 
is  quite  revealing  The  wavelength  of  the  jet  front  appears  to  bear  the 
same  relation  to  L  (about  Ap/L  **  9)  that  it  does  for  the  wake. 

It  snould  be  noted  that  the  physical  dimensions  of  the  iront  are  given 
here  in  physical  coordinates,  and  this  should  be  accounted  for  comparing 
e  present  results  with  other  experiments.  In  the  latter  the  radius  is 
usually  normalized  with  the  so-called  half-radius  rl/2.  Thus,  Becker  et  al 
(Reference  18)  gives  Y/r1/2  -  1.78,  Bradbury  (Reference  19)  1.73  and 
orrsin  and  Kistler  (Reference  11)  values  ranging  from  l.b5  to  2.00.  From 
the  present  experiment  (see  Figure  40) ,_we  obtain  Y/L  =  3.35  on  the  average 
an  in  terms  of  the  half-radius  we  get  Y/r1/2  =  2.0.  it  therefore  appears 
that  quite  a  few  jet  experiments  generally  agree  as  to  Y/n/9.  It  should 
be  remembered,  however,  that  heating  the  fluid  generally  increases  the  jet 
radius;  the  end  result,  shown  on  Figure  39,  is  that  our  Y  extends  farther 
out,  b>  about  30  percent,  than  Corrsin's  Y  for  example. 
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FIGIRE  VJ .  FRONT  WAVELENGTH  TOR  JETS  AND  WAKES  NORMALIZED  WITH  THE  TRANSVERSE 
SCALE  FOR  EACH  F*LOW 


The  value  of  o/ Y  obtained  here  is  0.23,  not  far  from  Bradbury  s  value  oL 
0.22  or  Corrsin's  range  of  0.23  to  0.29.  However,  Becker  et  al  obtained 
CJ/Y  =  0.16.  This  might  be  interpreted  as  a  difference  in  intermittent 
properties  between  an  inert  contaminant  and  the  temperature  lie  Id,  but  it 
is  more  likely  that  Becker's  method  of  measurement  and  perhaps  differences 
in  the  turbulent  Reynolds  number  are  responsible.  Beckci  and  his  co-workers 
have  at  any  rate  obtained  a  unique  relation  between  the  intermittency  tactor 
and  the  contaminant  density  distribution  with  radius.  In  the  reacting 
jet  employed  here  the  electrons  are  distributed  radially  in  an  irregular 
manner  because  of  reactions  (see  Figure  38)  and  this  relation  does  not 
mater  ial ize . 


4.2.3  CONCLUSIONS 

The  main  question  posed  in  this  section  has  been  answered:  the  front 
characteristics  are  identical  for  the  electron-sensitive  and  the  temperature 
sensitive  probe.  The  implication  is  that  the  intermittency  characteristics 
of  the  hypersonic  wake,  as  predicted  from  gasdynamical  studies  alone,  are 
also  the  intermittency  characteristics  relevant  to  the  radar. 


4.3  ELECTRON  DENSITY  FLUCTUATIONS 

Data  reduction  has  been  completed  on  the  magnitude  of  the  electron  density 
fluctuations  in  the  plasma  jet,  raw  data  on  which  were  presented  in  last 
year's  report  (Reference  b) .  The  final  data  reduction  was  accompanied  by 
further  experimental  measurements  needed  to  "fill  in"  portions  el  interesting 
electron  behavior.  In  this  section,  the  following  is  discussed:  (1)  Theore¬ 

tical  calculations  on  the  validity  of  the  Langmuir  probe  technique,  (2)  the 
final  form  of  the  rms  electron  density  fluctuation  magnitude  as  distributed 
axially  and  radially,  (3)  data-reduct ion  of  the  electron  density  spectra, 
and  (4)  studies  and  interpretations  of  the  very  large  fluctuations  observed 
in  the  transition  region.  Comparison  with  other  measurements  will  also  be 
g  i  ve  n . 


4.3.1  LANGMUIR  PROBE  SHEATH  CALCULATIONS 


The  radius  of  the  electrostatic  sheath  surrounding  the  Langmuir  probe  used 
in  this  work  is  crucial  to  the  validity  of  the  data.  Specifically,  the 
sheath  ra  '  ius 


a«Ae 


(51) 


where  A ,  is  the  electron  mean- free- path .  The  earliest  computations  of  a, 
under  this  contract,  were  based  on  a  combination  of  the  Langmuir-Mott-Smith 
and  Langmuir-Blodgett  theories  (Reference  20).  More  recently  tins  computa¬ 
tion  was  replaced  by  one  performed  by  Bettinger  and  Walker  (Reference  21) 
whose  attractiveness  rests  of  its  use  of  integral  and  c  uiservation  laws. 
Sutton  (Reference  22)  has  proposed  that  the  sheath  is  computed  on  the  basis 
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of  work  done  by  Laframboise  which  is  basically  an  extension  of  original 
work  by  Langmuir.  To  compare  these  choices  for  computing  a,  the  following 
formulas  were  plotted: 


(1)  Langmuir,  Mott-Smith  and  Blodgett: 


3/2 


(2)  Bettinger  and  Walker: 


(52) 


(r  '  ‘)  ‘"if  -  >)  “  '•“  (?)  ^ 

(3)  Laframboise: 


(53) 


(54) 


where  r  and  h  are  the  probe  and  Debye  radius,  respectively,  and  T)  the  non 
dimensional  probe  voltage  Tj  =  qv/kT.  The  former  two  equations  are  valid 
for  h/r  >1,  whereas  the  latter  presumably  for  0.1<h/r<2. 


These  relations  have  been  plotted  on  Figure  44.  They  indicate  that  the 
Langmuir-Mott-Smith-Blodgett  approach  gives  sheaths  which  are  in  agreement 
with  the  other  theories  for  h/r  =  1,  but  much  higher  than  the  latter  for 
h/r>>l.  Because  this  theory  is  also  the  least  sophisticated  (and  pre¬ 
sumably  least  accurate)  of  the  others,  it  is  concluded  that  the  original 
estimates  of  a  were  unduly  pessimistic;  the  results  on  electron  fluctuations 
discussed  in  this  report  are  thus  more  valid  than  originally  thought. 

4.3.2  MAGNITUDE  OF  ELECTRON  DENSITY  FLUCTUATIONS 


Figures  45  and  4b  summarize  the  final  results  of  the  electron  fluctuation 
magnitude  in  the  turbulent  plasma.  From  these  two  figures,  it  is  clear  that 
the  electron  fluctuations  undergo  two  phases:  up  to  about  24  inches  from 
the  nozzle  the  fluctuations  are  extremely  large  on  the  axis  (about  5  times 
the  mean)  whereas  beyond  that  distance  they  decrease  to  a  small  percentage 
of  the^  mean  and  develop  off-axis  peaks.  The  former  will  be  called  "transi¬ 
tional"  or  "heterogeneous"  fluctuations  whereas  the  latter  will  be  called 
"asymptotic"  or  "homogeneous" . 


Perhaps  the  single  most  interesting  phenomenon  observed  during  the  turbulence 
measurements  was  the  very  large  magnitudes  of  the  electron  density  fluctua¬ 
tion.  It  is  statistically  possible  for  the  standard  deviation  of  a  stationary 
variable  to  exceed  the  mean  value  of  the  variable  when  the  amplitude  distri¬ 
bution  of  the  latter  is  extremely  skew.  Such  situations  exist  in  the 
intermittent  zones  of  ordinary  turbulent  jets  for  any  flow  variable  Q  which 
fluctuates  between  nil  in  the  ambient  ("irrotat ional")  environment  and  some 
mean  level  within  the  turbulent  front.  If  the  latter  does  not  decrease 
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FIGURE  44. 


COMPARISON  OF  THREE  THEORETICAL  PREDICTIONS  OF  THE  SHEATH  RADIUS 
Rs  IN  TERMS  OF  THE  PROBE  RADIUS  Rp  AND  DEBYE  DISTANCE  H. 
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FIGURE  45.  AXIAL  VARIATION  OF  THE  RMS  ELECTRON  DENSITY  FLUCTUATIONS , 

NORMALIZED  WITH  THE  LOCAL  ELECTRON  DENSITY,  ON  THE  JET  AXIS 
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too  last  with  increasing  radius  then  il  ls  well  known  that  AQ/(j  can  diverge 
with  radial  distance.  In  tact  Becker,  liottel  and  Williams  (Reference  18) 
have  observed  this  divergence  by  utilizing  an  inert  contaminant  concentra¬ 
tion  for  an  observable. 

The  possibility  that  the  large  fluctuations  on  the  axis  pictured  on 
Figure  are  caused  by  the  intermit  tency  is  excluded  by  the  simple  fact 
that  the  flow  on  the  axis  is  Lully  turbulent  as  already  demonstrated.  From 
the  results  thus  far  shown  in  fact,  it  appears  that  we  have  here  a  fully 
turbulent  tluid  with  normal  temperature  fluctuation  but  with  highly-skewed 
electron  fluctuations.  An  osc il loscop ic  study  of  the  turbulent  jet 
immediately  showed  large  "spikes"  in  the  electron  density  time  history  at 
a  point;  these  spikes  are  indeed  what  statistics  would  require  to  produce 
the  large  An/n  of  Figure  45 .  Figure  47  shows  a  comparative  iscillnscope 
study  of  the  Langmuir  probe  and  hot-wire  output  at  various  radial  positions 
but  at  the  same  distance  from  the  jet  nozzle.  By  proper lv  adjusting  the 
scope  gain  it  is  possible  to  show  conclusively,  tirst,  the  differences 
between  intermittent  and  fully  turbulent  l  1  >w  and,  secondly,  the  differences 
between  Langmu ir-probe  and  hot-wire  output  in  the  lully  turbulent  region. 
Attention  is  drawn  to  the  large  spikes  in  the  Langmuir  probe  traces  which 
are  clearly  absent  t rom  the  hot-wire  traces. 

The  appearance  of  the  large  spikes  persisted  lor  several  inches  in  the 
upstream  end  ot  the'  turbulent  zone  studied,  but  decayed  and  eventual  lv 
vanished  as  the  probe  was  moved  downstream;  the  lluctuation  level  An/n  alsi 
decayed  bel  )w  unity  as  Figure  o5  shows.  To  complete  the  picture,  an 
investigation  ot  the  laminar  and  transitional  portions  of  the  jet  was  made 
using  a  Langmuir  probe  thicker  chan  described  previously  in  this  paper,  so 
that  it  could  withstand  the  much  higher  heat  t  luxes  in  the  laminar  jet.  For 
obvious  reasons  the  output  of  this  probe  was  not  quantitatively  interpretab  le , 
but  it  made  possible  to  piece  together  a  reasonable  description  it  the  events 
leading  t>  the  generation  ot  the  large  fluctuations. 

Figure  48  shows  this  picture.  The  Langmuir  probe  observation  made  is 
typified  by  an  instantaneous  "snapshot"  of  the  electron  density  on  the  ,a 
axis  from  the  nozzle  exit  well  into  the  turbulent  zone.  In  the  laminar  jet 
the  mean  electron  density  is  high  with  an  unsteady  component  consisting  of 
the  rectifier-driven  360  Hz  modulation  whose  rms  magnitude  is  about  10  percent 
ot  the  mean.  At  the  flow'  velocities  prevalent  here  the  wavelength  it  this 
sinusoidal  Lluctuation  is  actually  very  much  larger  than  the  length  it  the 
laminar  j_..  (anuu1  10  feet  compared  to  about  1  toot).  In  the  transition 
region  the  breakdown  ’Zocess  mixes  ambient  fluid  with  the  ionized  tluid; 
the  former  is  <  ^ectron-poor  so  that  large  "negative"  pulses  (ii  one  thinks 
of  the  electron  current  as  producing  a  positive  displacement  of  the  scope 
beam)  appear  as  shown  on  Figure  -t8:  the  signal  is  thus  "negative  lv  one¬ 
sided  The  mean  electron  density  thus  decreases  be  1  ov  the  maximum  level, 
which  also  decreases  because  of  recombinations.  Farther  downstream  the 
spikes  widen  (hut  do  not  necessarily  deepen),  and  still  farther  the  signal 
begins  reversing  itself,  i.e.,  becoming  "positively  me-sided".  In  this 
latter  case  the  mean  electron  density  drops  considerable  below  tin  tips  ot 
the  spikes,  the  fluctuation  distribution  is  very  skew  and  the  rms  fluctuation 
level  is  considerably  higher  than  the  mean  level.  It  is  at  this  stage  that 
the  oscillograms  shown  on  Figure  47  were  taken. 
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FIGURE  47.  OSCILLOSCOPIC  STUDY  OF  TEMPERATURE  AND  ELECTRON  DENSITY  FLUCTUATIONS 
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FIGURE  48.  SCHEMATIC  OF  TYPICAL  INSTANTANEOUS  ELECTRON  DENSITY  DISTRIBUTE 
ALONG  THE  JET  AXIS  AS  INFERRED  FROM  PROBE  MEASUREMENTS 
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It  should  agaui  be  emphasized  that  at  this  stage  the  lluul  is  iully  turbu¬ 
lent  from  the  standpoint  of  intermit  tency,  and  t  ha  t  turbulent  fluctuations 
in  the  electron  density,  of  order  20  to  40  percent  of  the  mean,  exist  in 
addition  to  the  large  electron  density  spikes.  The  latter  are  obviously 
th'’-  remnants  of  the  laminar  flow  upstream  of  the  transition  zone,  in  which 
the  electron  density  level  is  very  high  when  compared  to  the  level  in  the 
turbulent  fluid  itself.  Eventually  these  spikes  also  disappear  and,  as 
Figure  45  shows,  the  ratio  \u/n  settles  to  values  oi  order  10  percent  in 
the  far  jet. 

In  Figure  48  it  is  schematically  shown  that  the  maximum  ot  the  normalized 
fluctuations  is  attained  where  the  mean  electron  level  is  much  below  the 
maximum  level  1%,  that  is,  the  locus  of  the  spike  maxima.  If  the  spikes 
are  truly  the  results  of  the  process  of  engulfment  oi  the  ambient,  electron 
poor  fluid  by  the  electron-rich  fluid  of  the  jet,  and  il  tin  electron 
density  of  the  former  is  nc ,  then 


where  P/fu  is  the  ratio  of  the  total  mass  to  the  jet  mass.  Here  it  has  been 
assumed  that  both  fiuid  components  are  laminar  so  that  the  signal  can 
locally  oscillate  between  nm  and  nL  only,  in  the  manner  suggested  some  time 
ago  by  Feldman  and  Proudian  (Reference  2J)  and  Lin  (Reference  24).  The 
type  of  fluctuations  of  Equation  (55)  are  nil  be  tore  entrainment  begins 
(P/m  =  1)  and  also  when  P/m  tends  to  infinity  in  the  tar  jet,  and  they 
peak  in  between  with  a  magnitude  depending  on  nm/nc .  For  equilibrium 
chemistry,  when  nin  and  nc  are  directly  obtainable  from  the  corresponding 
temperatures  Tm  and  Tc ,  it  is  shown  by  Reference  15  that  An/ n  of  order  10 
or  higher  can  obtain  for  practical  values  of  Tm/Tc.  For  finite  reaction 
rates  the  ratio  An/n  becomes  extremely  sensitive  to  the  chemistry. 

Nevertheless  the  simplicity  of  the  jet  flow  brought  about  by  gasdynamical 
self-preservation  allows  certain  broad  tests  of  Equation  (55)  to  be  made 
easily.  For  example  the  entrainment  rate  is  given  theoretically  and  has 
also  been  measured  directly.  For  recombination  reactions  the  nni  might  also 
be  obtained  as  a  function  of  x.  The  normalized  electron  density  fluctuations 
of  Equation  (55)  thus  become  dependent  only  on  the  ratio  nmo  (ot  the  maximum 
electron  density  in  the  laminar  flow)  to  the  "ambient11  value  nc .  The  former 
is  actually  obtainable  from  the  experiment  but  the  latter  is  not,  since  the 
electron  density  at  the  base  of  the  spikes  was  detected  to  be  turbulent  and 
oi  a  mean  value  much  higher  than  that  outside  the  jet.  Computations  made 
with  Equation  (55)  using  estimated  values  of  nc  ga\e  reasonable  numerical 
agreement  with  the  measured  axial  variation  of  An/n  shown  on  Figure  45. 

Equation  (55)  thus  appears  to  be  a  reasonable  descriptor  of  the  physical 
mechanism  giving  rise  to  the  large  electron  density  fluctuations  -  on  the 
strength  of  the  mixing  model  it  describes.  It  does  not  necessarily  endorse 
the  "laminar  mixing"  ideas  of  References  (23)  and  (24)  since  the  partici¬ 
pating  lumps  of  fluid  (whether  electron-rich  or  electron-poor)  are  turbulent 


-89- 


in  character.  At  this  time  it  also  appears  to  affect  the  transitional 
zone  more  so  than  the  turbulent  activity  in  the  tar  jet;  in  t lie  present 
experiment  this  zone  is  only  a  few  jet  thicknesses  long. 

Once  beyond  the  region  of  very  high  live  l,  the  electron  density  fluctuations 
decrease,  according  to  Figure  45,  to  values  more  typical  of  homogeneous 
mixing.  The  following  points  are  significant:  first,  the  axis  values  ot 
\n/n  do  not  appear  to  attain  an  asymptote.  This  may  be  attributed  either 
to  the  tact  that  the  overall  jet  fluctuations  themselves  have  not  relaxed 
to  the  self-preserving  state;  or,  which  is  more  likely,  that  no  sclL- 
preservation  of  the  electron  fluctuations  attains  where  all  other  types  of 
fluctuations  have  become  self-preserving  (e.g.,  constant  fractions  oL  the 
mean  values).  Secondly,  the  oti-axis  maxima  of  these  tluctuations  (see 
Figure  4b)  occur  in  the  general  vicinity  of  T)  =  1.3.  By  comparison,  the 
region  of  maximum  shear  in  the  jet  occurs  at 


which,  using  the  sell-preserving  velocity  distribution 

,2 


u  =  e 


-0.24T 


is  seen  to  occur  at 

1/2 


T)  ,, 


1 


0.48 


=  1.45 


(5b) 


(57) 


(5b) 


It  thus  seems  that  the  maximum  of  electron  density  tluctuation  occurs  at 
or  near  the  same  radial  position  in  the  '  t  occupied  by  the  maximum  shear 
region.  This  is  an  important  finding  with  obvious  applications  to  the 
uypersonic  wake  problem. 


At  radial  positions  beyond  7)  =  1.45,  the  tluctuation  level  does  not  decay 
as  rapidly  as  expected.  In  fact  it  appears  that  I rom  the  standpoint  of 
fluctuations  the  jet  is  much  wider  than  evidenced  bv  the  mean- flow  measure¬ 
ments.  There  is  presently’  no  explanation  ot  this  phenomenon  available  other 
than  the  effect  of  intermit tency .  As  is  already  shown  in  Paragraph  -*.2, 
the  intermittent  electron  flow  extends  considerablv  larthor  Lrom  the  axis 
than  the  distribution  of  electron  density.  For  example,  at  x  =  30  inches, 
the  flow  is  intermittent  still  at  T)  =  5.5,  whereas  the  electron  density  has 
reached  ambient  (zero)  value  by  T)  =  4.5.  Now  in  intermittent  zones  the 
"turbulent'  signal  is  made  up  o t  a  combination  of  real  turbulence  and 
intermit  tency  in  a  way  given  functionally  by  Corrsin  and  Kistler  (Reference  11). 

In  the  example  cited  the  interm  it  tency  contribution  should  be  prevalent  m 
the  outer  region,  say  in  4<7)-  b.  To  test  this,  Figure  49  plots  the  \n/n 
distribution  at  x  =  30.  Also  plotted  is  the  signal  expected  it  it  consisted 
of  the  intermit  tenev  contribution  alone,  as  given  by  the*  tormulas  developed 
in  Appendix  B  of  Reference  4:  \ 


FIGURE  49  TYPICAL  RADIAL  VARIATION  OF  NORMALIZED  RMS  ELECTRON  DENSIT\ 
FLUCTUATIONS  ALSO  SHOWING  THE  SUSPECTED  CONTRIBUTION  OF 
INTERMITTENCY 
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where  7  is  the  intermittency  factor.  Near  the  axis,  where  7=  1  the  signal 
is  obviously  all  due  to  turbulence,  but  the  curves  cross  around  TJ  =  2.3.' 
This  is* surpr is ing  since  it  says  that  rear  the  jet  edge  the  signal  is 
smaller  than  the  minimum  expected.  The  question  is  open  at  this  juncture 
A  likely  explanation  is  that  Equation  (39)  is  oversimplified  because  it 
assumes:  (1)  no  electrons  outside  the  jet,  (2)  an  intermittency  signal 

i  ea  ized  by  flat  top  distribution,  and  (3)  a  jump  in  electron  density 

across  the  front  which,  at  a  point  fixed  near  the  jet  edge,  is  constant 
in  time. 


4.3.3  SPECTRA  OF  ELECTRON  DENSITY  FLUCTUATIONS 

Ih  the  previous  reporting  period,  spectra  of  the  electron  density  fluctua¬ 
tions  were  taken  at  ail  x  and  Y  positions  with  a  wave  analyzer  whose  minimum 
bandwidth  was  200  Hz.  These  data,  reduced  during  the  current  period,  showed 
two  deficiencies  connected  with  the  very-low  frequency  range.  First,  many 
of  the  spectra  were  recorded  with  an  x-Y  reorder  gain  such  that  the  low- 
frequency  spectral  density  went  off  scale.  Secondly,  it  was  evident  that 
the  spectral  density  dropped  off  rapidly  in  the  first  few  hundred  cycles, 
so  that  measurements  in  that  region  were  inaccurate  due  to  the  200  Hz 
bandwidth.  Thus,  the  data  could  only  be  considered  useful  above  a  certain 
frequency,  which  was  set  at  2  kHz. 

To  fill  the  low- frequency  gap  in  each  spectrum,  all  spectra  were  re-taken 
in  the  0-2  kHz  range  using  the  setup  of  Figure  50.  A  Hewlett-Packard 
Model  No.  30 2 A  Wave  Analyzer  was  used  with  a  range  of  50  kHz  and  bandwidth 
fixed  at  7  Hz.  This  bandwidth  could  resolve  the  spectral  density  at  much 
lower  frequencies  than  before,  but  it:  also  introduced  the  classic  output 
unsteadiness  typical  ot  sharply  filtered  low-frequency  spectral  analysis. 

To  damp  this  out,  a  vacuum  thermocouple  was  interposed  between  the  wave 
analyzer  dc  output  and  the  plotter.  Figure  51  shows  the  transfer  function 
ot  the  thermocouple  which  was  used  in  the  computerized  data  reduction 
process  to  connect  the  ut  appearing  at  the  plotter  input  to  the  filtered 
output  of  the  wave  analyzer.  With  this  equipment  spectra  were  recorded  from 
0  to  2  kHz.  Each  ot  these  spectra  was  then  joined  with  its  corresponding 
predecessor  (taken  with  the  previous  arrangement)  at  the  2  kHz  point.  In  the 
process,  account  was  taken  of  the  fact  that  the  preceding  spectra  were  in 

terms  of  the  root-mean-square,  while  the  most  recent  ones  in  terms  of  the 
mean-square  spectra. 

The  data  were  reduced  by  the  JEA-II  (SRS  0024)  computer  program  which  is 
described  in  Appendix  A.  This  program  gives  an  independent  check  of  the 
total  rms  fluctuation  by  integrating  the  area  under  each  spectrum.  After 
the  data  was  first  reduced  it  developed  that  this  cross-check  did  not  con¬ 
form  with  expectations.  It  was  found  that  the  Runge-Kutta  integration 
scheme  used  was  not  accurate  for  functions  changing  rapidly  between  adjacent 
requenc les .  The  data  were  therefore  processed  again  through  the  program 
using  a  Simpson-type  integration. 

The  electron  density  fluctuation  spectra  arc  shown  on  Figures  52  through  58 
The  first  two  figures  show  the  spectra  on  the  axis;  the  remainder  show 
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FIGURE  50.  BLOCK  DIAGRAM  OF  THE  JET  TURBULENCE  ELECTRON  ICS 
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FIGURE  51.  TRANSFER  FUNCTION  OF  VACUUM  THERMOCOUPLE  USED  TO  CONVERT  THE 
SPECTRAL  DENSITY  SIGNAL  INTO  DC  VOLTAGE 


-94- 


he  f 

U 


B1446U 


FIGURE  54.  NORMALIZED  SPECTRAL  DENSITY  OF  ELECTRON  FLUCTUATIONS  AT  18  INCHES 
FROM  THE  NOZZLE 
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FIGURE  55. 
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NORMALIZED  SPECTRAL  DENSITY  OF  ELECTRON  FLUCTUATIONS  AT  22  INCHES 
FROM  THE  NOZZLE 
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FIGURE  56.  NORMALIZED  SPECTRAL  DENSITY 
FROM  THE  NOZZLE 


OF  ELECTRON  FLUCTUATIONS  AT  26  INCHES 


representative  spectra  at  off-axis  positions.  The  ordinate  in  these  graphs 
is  in  the  square  of  the  electron  fluctuation  spectral  density  and  is 
normalized  to  the  asymptotic  value  at  zero  frequency.  The  abscissa 
ij>  =  f^e/u  is  the  frequency  made  dimensionless  with  the  local  integral  scale 
of  electrons  Ae  and  the  local  velocity*  u.  For  example,  f^  =  0.1  relers 
to  an  eddy  size  ten  times  larger  than  the  integral  scale. 

The  axis  spectra  evolve  along  a  characteristic  pattern.  According  to 
Figure  52,  nearer  the  nozzle  the  decline  of  spectral  density  with  frequency 
is  quite  gradual,  but  becomes  most  precipitous  at  the  four  XSTATIONS.  At 
the  latter,  the  spectra  appear  to  become  "locked"  to  an  invariant  shape, 
for  which 


f 
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beyond  *.^>1.  .t  is  further  demonstrated  by 
thi  spectral  shapes  occurs  mainly  between  x 
is,  in  the  region  where  the  transition  from 
turbulence  occurs. 


Figure  53  that  the  change  in 
=  24  inch  and  x  =  28  inch,  that 
"heterogeneous"  to  "homogeneous" 


In  the  radial  direction  the  change  in  spectral  shapes  is  generally  gradual 
and  small,  especially  at  large  X.  In  the  "heterogeneous"  portion,  the  spectral 
density  tends  to  decrease  at  the  higher  frequencies  as  on®  moves  away  trom 
the  axis.  At  larger  X  a  slight  but  repeatable  reversal  occurs,  so  that  the 
higher  frequencies  are  more  prominent  on  the  axis  and  again  near  the  jet 
edge  and  less  so  in  between. 


It  will  be  noted  that  in  all  graphs  shown  there  is  no  portion  of  any  spec¬ 
trum  where  the  spectral  density  decays  with  the  (-5/3)  power  of  frequency 
as  expected  for  an  inertial  subrange.  This  matter  will  be  discussed  in 
detail  in  Paragraph  4.5. 


The  longituuinal  scales  (macroscales)  ol  electron  iluctuations  are  derived 
in  the  JEA-II  program  irom  the  measured  spectra  by  the  usual  formula: 


1 1 

e 


u  i^Sli 

—  9 

4  (\n) 


(60) 


Here  u  is  the  local  velocity  and  (\n)  *"  is  the  total  ( i requenev- integrated) 
mean-square  fluctuation  in  units  of  cm’b.  The  (\n(0))~  is  the  contribution 
to  (Xn)^  which  comes  from  a  one-cps-wide  pussband  at  t  =  0;  it  has  the 
units  (cm*^)  x  (cps)”^  *  (cm'*1)  (sec”^)~^  =  (cm  )  (sec),  and  has  been 
obtained  irom  the  spec*,  al  density.  Thus,  **c  has  the  units  of  length,  as 
it  should.  To  double-che^k  the  validity  of  (\n(0))~,  the  program  pertorms  the 
following  integration  from  the  measured  spectrum  (\n(f))-  versus  t  at  each 
po int : 


/*  °°  2  d  f 

EINT  =  J  (\n( f ) )  — 

B 


(61) 


Since  Ae  is  proportional  to  u  also,  this  nond  ime ns  tonal  frequency  is 
independent  of  u. 
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voltmeter  and  shown  on  Figure  . 

The  scales  so  measured  are  shown,  °» 

following  observations  are  ma  •  heterogeneous  region,  decreases 

(2)  initially  the  scale  i.  higl tn  th ^  ^  reRlon.  Bccause  of  the 

abruptly  and  increases  slowl*  ™  .  after  the  spectral  measurements, 

scat  .er  a  cross-check  was  Per  f ormed  ,^te" Darate  iy .  The  cross-check  data 
where  (\n)2  and  (  n(0))  we*e  generally  higher  values,  somewhat 

(tilled- in  symbols  on  Figure  5  ),  g  turns  out  that  the  average  of 

less-scattered,  than  the  spectra  TaklnR  thc  last  six  XSTATIONS 

both  measurements  has  stl_  1  £oUowinB  straight  line  was  fitted  y 

and  the  point  (x  =  0,  “  u' 

least  squares: 


(62) 


a  =  0.034  x 
e 


,  and  rhe  transverse  scale  L 

and  the  following  proportionality  between  *e 

of  the  jet  was  found: 

(63) 

»  =  0 .47  L 

e 

i  ♦- t-hp  scaLc  magni.- 

As  will  be  seen  below,  these  formulas  _comc  the  scatter 

tude  for  jets  in  the  absence  °  [  \  he  oxperiment  (see  Reference  25), 

typical  of  the  technique  rather  than  P  correlation  measurements, 

is  much  worse  than  expected  tram,  say ,  t  , 

i  ,i rnadv  shown  in  last  )ear  s 

Radial  distributions  of  electron  = -[-^^ity  measurers  are 
report  (Reference  6)  are  oeuu 
c omp lete . 

4.4  the  gas  temperature  fluctuations 

i.4.1  FLUCTUATION  MAGNITUDES 

The  frequency—  integrated  ™aSnitudes  of  the  “^^^^^tn'Reference  6, 
together  with  a  discussion  °£  ,tht  '  jnal  data  were  taken  to  corroborate 

page  181.  in  the  present  pe.iod  ^‘^a,,,  compared  with  ear  let 

t„e  -tuer  results  t^  ratories,  ^  certain  rmprovements  in  techntqu 

measurements  in  otnci 

were  made. 

The  temperature  fluctuations  are  shown  on  flR"CS  fjecrease  of  the  mean 
these  decrease  rn  a  rate  nearly  pro po «  ^tlon  the  familiar  off- 

-=s:i“nuct^ln‘nsi^  is  noted,  to  a  degree  dependent 
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FIGURE  60.  AXIAL  VARIATION  OF  VELOCITY  AND  GAS  TEMPERATURE  FLUCTUATIONS 
NORMALIZED  WITH  THEIR  CORRESPONDING  AXIS  "DEFICITS" 
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strongly  on  the  manner  by  which  AT  is  normalized.*  There  are  expected 
similarities  between  these  results  and  those  of  Corrsin  and  Jbero 
e,peciaUy  in  the  »Uo  AT/ (T-T-)  which  attatn^^lue^f  about^.^ 

(corre  la  t  ion)  ‘  probe^ar^a  ls^  shown  on  figure  bO.  These  peaks  originally, 

but  eventually  ““  I  ”  ^^.““^.“.2  and  the 

Hr-  £rom  a 

Hi: 

d  iame  ters . 


A. 4. 2  THE  TEMPERATURE  FLUCTUATION  SPECTRA 


During  the  current  period,  repeated  measurements  of  the  temperature 
^t^/spectra  were  made  In  these  “^2.“  ^ 

r^SS^t^S:  hit-wire  ^"constant  that  is  with  the  change 

^he 'transfer  function  (frequency  response)  of  the  detectinS  equipment  n 
wav  difficult  to  detect  and  compensate  for.  in  practice  t  e  pro 
U»t  -»«-  outside  the  j-l  ^  Jherma.  Unmeasured  by 

*  -  a. 

compensation,  deviations  of  the  time  constant  Iron  .  se“  u 

resistance  through  the  jet. 

In  contrast  to  results  obtained  with  the  wake,  the  technique  is  not  easy 

to  use  in  the  jet.**  For  the  0 .00014- inch-d iameter  wire  used,  a  time- 

narant  oi  4  5  milliseconds  was  measured  outside  the  jet.  Initially  the 
constant  oi  4.0  nuuiac'-u'  _  :11  •  ...-He  in  the  flow. 

than  about  5  milliseconds,  and  that  there  was  no  increase  at  x  -  30 
32?  34  and  36  inches.  Anomalies  at  x  =  28  inches  were  noted  and  these 

data  are  not  discussed  herein. 

. » ;-,rr;s  nx,r 

the  X-Y  plotter.  In  the  range  from  0  to  2  kHz,  the  sig 


**The  difficulties 
at  extremely  sma 
res  pec  t ive ly) . 


arise  in  the  problems  of  heat  transfer 
11  Mach  and  Reynolds  numbers  (typically) 
A  technical  report  on  this  problem  is  co 


from  cylinders 
0.01  and  0.001 
ntemplated . 


> 
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outside  the  7  Hz  band  ol  a  J02A  Wave  Analyzer,  whose  restored  output  is 
channeled  to  a  vacuum  thermocouple  and  from  there  to  the  plotter.  Capaci¬ 
tor  damping  at  the  plotter  input  was  used  tor  stabilizing  the  signal 
further.  The  circuit  is  shown  in  Figure  50. 

As  with  the  electron  spectra,  spectra  were  recorded  every  two  inches  along 
the  jet  axis  from  x  =  24  inches  to  x  =  Jb  inches.  The  Lorward  limit, 
x  =  24  inches,  was  chosen  for  Lear  that  closer  to  the  nozzle  the  hot-wire 
would  melt.  At  each  XSTATION,  spectra  were  recorded  at  radial  distances 
snaced  about  0.75  inch  apart.  The  data  were  reduced  with  the  aid  ol  the 
SPECTR  program  written  for  the  General  Electric  Time-Sharing  System.  An 
access  terminal  to  this  computer  is  available  on  an  almost  Lull-time 
basis  for  this  work.  Details  of  the  program  appear  on  Appendix  B. 

A  summary  of  the  temperature  spectral  results  appear  on  Figures  62  through 
65.  As  with  the  electron  spectra,  the  ordinate  in  these  graphs  is 
normalized  with  the  spectral  density  at  zero  frequency.  The  frequency  has 
been  non-d imens iona l ized  with  the  integral  scale  ' j  (which  will  be 
discussed  later)  and  the  local  axial  velocity  u.  Figure  62  shows  the 
evolution  of  these  spectra  on  the  axis,  while  the  other  figures  show  their 
radial  development  at  selected  XSTATIONS.  The  spectra  are  remarkably 
similar  to  each  other  and  not  too  different  from  expectations.  They  seem 
to  fall  into  three  broad  categories: 

Type  I:  Nearer  the  nozzle  and  close  to  the  axis  the 
high-frequency  content  is  intense.  The  spectrum  goes 
as  fjj'5/3  where  fjj  =  Ajf/u. 

Type  II:  Over  the  rest  of  the  entire  jet  but  not  too 
close  to  the  edge  the  data  collapse  into  a  behavior 
showing  a  gradual  decrease  of  the  spectral  density 
without  any  distinct  regions  where  it  follows  the 
(-3/3)  slope.  Typically  at  f^  =  4,  the  spectral  con¬ 
tent  is  down  a  factor  of  10  below  the  Type  I  spectra. 

Type  III:  At  the  edge  of  the  jet,  the  high-frequency 
content  is  still  lower  and  the  spectral  decay  some¬ 
what  steeper. 

The  evolution  of  the  spectra  from  one  type  into  another  is  best  illustrated 
on  Figure  63.  It  is  very  clear  that  on  and  near  the  axis  one  obtains  the 
(-5/3)  slope.  At  Y  =  1.47,  the  "break"  away  from  this  behavior  (Type  I) 
occurs  at  earlier  frequencies.  Also,  the  departure  of  the  data  at  the 
jet  edge  (Y  =  6.64)  from  the  Type  II  data  is  obvious.  Tentatively,  the 
following  conclusion  is  reached  that  Type  I  spectra  are  associated  with 
the  "heterogeneous"  jet,  where  the  large  electron  density  fluctuations 
also  occur.  Type  II  are  apparently  connected  with  the  "homogeneous"  jet, 
where  the  electron  spectra  also  assume  a  constant  shape. 

An  attempt  was  made  to  compare  these  data  with  spectra  taken  in  other  jets. 
Of  the  latter  there  is  only  the  data  by  Corrsin  and  Uberori  (Reference  26). 
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FIGURE  64. 


normalized  spectral  density  of  the  temperature  fluctuations 
AT  30  INCHES  FROM  THE  NOZZLE 
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FIGURE  65.  NORMALIZED  SPECTRAL  DENSITY  OF  THE  TEMPERATURE  FLUCTUATIONS 
AT  34  INCHES  FROM  THE  NOZZLE 
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In  F inures  66  and  67,  the  data  are  instead  compared  with  more  Bemral 
turbulent8  flows  fr™  which  some  conclusion  about  the  present  data  can  be 
drawn.  It  is  seen  that  the  Type  1  spectra  arc J  comp  able  wc  ^ 

Reynolds  number  flows  where  isotropy  cs  expected.  At  lower  Ke  n l 
grid  turbulence  compares  well  with  Type  II  spectra.  The  tnterpolat 
formula  giving  a  dependence  f„**  does  not  correlate  well  either  type 

data. 

ThP  longitudinal  integral  scale  of  the  temperature  iluctuations  is  shown 
on  Figure  68.  As  with  the  electron  scales  the  method  ol  measurement, 
which  utilizes  the  spectra,  results  in  rather  large ■  st®  "•  tit  by 

are  supposed  to  grow  linearly  with  distance,  and 

the  straight  line  si  own  has  been  forced  to  go  through  the  virtual  origin 
of  the  jet  (x  =  12  inches).  The  linear  growth  is  actually  ^Ftwrte 

fact  by  the  near-proportionality  of  *T  with  the  transverse  scale  L  also 
shown  in  Figure  68.  The  following  relations  emerge. 


0.079 


0.79  L 


(64) 


(65) 


A  rjn 

It  should  be  noted  that  these  two  relations  imply  that  L  =  0.076  x,  which 
is  not  far  from  the  result  found  lor  this  jet  in  Reference  2,  page  78. 


L  =  0.079  x 


(66) 


A  comparison  of  this  result  with  previous  experiments  is  also  shown^on^ 
Tv,e  excellent  experiments  of  Wygnanski  and  Fiel  ( 

BB  in  ;  have  b^n  chosen  L  the  purpose.  These  data  are  chosen  because 

of  their^clear  superiority  In-... « log  techn ^^"f^sc^io!, 

and  also  because  they  were  obtained  in  t  J  Wvtmanski  and 

has  been  attained.  The  velocity  scale  growth  observed  b>  Uygnanski 

Fielder  is  summarized  as  follows: 


L  =  0.053  x 


u 


=  0.0385  x 


(67) 

(68) 


in  comparing  the  present  scale  growth,  Equation  (64)  with  the  Boeing  result, 
Equation  (68),  it  should  be  pressed  that  the  P«.,nt  Je, ‘^highly ^heat^, 

A  i^the^ve  loc  ity  -  if  levels  easily  explain^ 

the  apparent  and  contrast  of  Figure  68.  First,  the  prc  -ent  J 
because  of  the  initially  very  high  temperature,  an  e  feet  pre  “ 
observed  by  Corrsin  and  Uberoi  (Reference  26)  and  analyzed  by  L 
(Referenced)  and  also  immediately  apparent  from  Equation  (66)  an  (  >• 

f  "corrected”  for  the  temperature  effect  the  scale  Au  obtained  by  Boeing 
Ion  Id  move  up  by  a  factor  of  0.079/0.053,  from  Equations  (66)  and  (67). 


A  = 


0 ,0—  0.0385  x  =  0.057  x 


(69) 


0.053 
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COMPARISON  OF  TYPE  I  TEMPERATURE  SPECTRA  WITH  HIGH- REYNOLDS 

™“ru^S0TR0PIC  VELOcm  SPECTRA  AND  the  karhan  interpolation 


PRESENT  DATA  WITH 
-5/3  SLOPE  (TYPE  I) 


PRESENT  DATA  (TYPE  II) 


o  FAVRE  GRID 
(ReM  -  21,500) 
VELOCITY  SPECTRA 


u  B1459U 

FIGURE  67.  COMPARISON  OF  TYPE  U  TEMPERATURE  SPECTRA  WITH  TYFICmL 
LOW- REYNOLDS -NUMBER  VELOCITY  SPECTRUM 
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In  the  second  place,  however,  it  is  known  (Reference  26)  that  the  tempera¬ 
ture  s<  ale  is  smaller  than  the  velocity  scale,  so  that  value  would  again 
decrease  somewhat.  In  any  way,  the  agreement  with  the  Boeing  results  is 
good  in  v iew  of  the  scatter  in  the  present  data  and  the  possibility  that 
the  present  jet  has  not  yet  become  se 1 t - prese rv ing . 

Two  features  of  Figure  68  are  important  to  the  present  problem:  in  both 
the  Boeing  and  the  present  jet,  A/L  is  about  constant  at  0.75;  in  other 
words,  the  longitudinal  integral  scale  appears  to  be  a  constant  fraction 
of  the  jet  width.  Secondly,  this  ratio  is  apparently  nearly  the  same  in 
jets  and  wakes;  note  that  about  the  same  ratio  A/L  was  obtained  lor  the 
WEB  experiment  (Reference  4,  page  42). 

4.5  COMPARISON  OF  ELECTRON  AND  TEMPERATURE  FLUCTUATIONS 

As  it  was  noted  in  the  introduction  to  this  section  (Paragraph  4.1),  the 
ultimate  aim  of  this  task  is  to  establish  a  connection  between  electron 
and  temperature  fluctuations  in  a  turbulent  plasma.  Certain  results  in 
this  direction  are  now  available  and  arc  summarized  in  Figures  69  through  72. 

On  Figure  69,  the  ratio  (An/n) / (\T/T)  is  shown  on  the  jet  axis  beyond 
x  =  26  inch,  that  is,  in  the  region  in  which  temperature  fluctuation  data 
are  available.  The  electron  density  fluctuation  here  predominates, 
especially  nearer  the  nozzle,  where  ratios  as  high  as  5  are  found.  Note 
that  at  x  =  26  inches,  the  jet  flow  is  already  beyond  the  so-called 
"heterogeneous"  stage,  implying  that  ratios  much  higher  than  5  are  probable 
in  the  range  18  inch  v  x<  26  inch.  Toward  the  iaithest  point  of  the  axis 
(x  =  36  inch),  the  ratio  approaches  unity.  Near  the  jet  edge,  as 
Figure  70  shows,  the  ratio  is  as  high  as  10  or  so. 

Figures  71  and  72  plot  the  ratio  of  electron-to-temperature  fluctuations 
resolved  according  to  wave  number.  Physically,  these  graphs  indicate  the 
desired  ratio  separately  for  each  eddy  size.  The  two  temperature  spectral 
types,  I  and  II,  are  shown  with  their  electron  fluctuation  counterparts 
as  drawn  from  Figures  such  as  52,  62,  etc.  The  abscissa  lias  been  adjusted 
to  wave  number  units  by  essentially  removing  the  scales  Ae  and  AT  from  the 
abscissas  of  the  cited  figures.  Note  that  ir  each  eddy  size  the  electron 
density  fluctuation  also  predominates  over  the  corresponding  temperature 
fluctuation.  The  ratio  of  the  two  is  plotted  at  the  bottom  of  eacli  Figure. 
The  ratio  of  the  relevant  fluctuations  \n/n  and  \T/T  can  then  be  obtained 
from  the  following  formulas: 
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FIGURE  69.  AXIAL  VARIATION  OF  THE  RATIO  OF  TOTAL  NORMALIZED  ELECTRON  DENSITY 
FLUCTUATION  TO  THE  TOTAL  NORMALIZED  GAS  TEMPERATURE  FLUCTUATION 
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FIGURE  71. 


RELATIVE  IMPORTANCE  OF  ELECTRON  AND  TEMPERATURE  SPFCTR  \I 
DENSITIES  IN  THE  "HETEROGENEOUS"  PORTION  OF  THE  JET 
(ABSCISSA  IN  ARBITRARY). 
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FIGURE  72.  RELATIVE  IMPORTANCE  OF  ELECTRON  AND  TEMPERATURE  SPECTRAL 
DENSITIES  IN  THE  "HOMOGENEOUS”  PORTION  OF  THE  JET 


here  the 


nomenclature  is  as  follows: 


VKO  Spectral  density  of  Q  at  frequency  f 

XQ(0)  Spectral  density  of  Q  at  zero  frequency 

Q  Frequency- integrated  (total)  fluctuation  of  Q 


Q  Mean  value  of  Q. 

From  the  definition  of  the 


integral  scale 


we  obtain 


(\n(f)/nr 

( \T ( f ) /T) 2 


\n(f)/\n(o) 

\T(f)/\T(o) 


All/.Q 

At/t 


(72) 


(73) 


The  first  ratio  in  the  bracke 
the  second  ratio  is  typically 
ratio  Ae/AT  also  varies  from 
can  be  given  to  the  degree  of 
Equation  (65)  are  obtained: 


t  is  the  one  given  in  Figures 
given  in  Figure  69  on  the  jet 
point  to  point  in  the  jet,  but 
accuracy  with  which  Equation 


71  and  72,  while 
axis.  The 
on  the  axis  it 
(63)  and 


_  0.47L 
At  “  0 . 79L 


0.60 


(74) 


is  cons  iderably3^  igher  tha^uni^v"  ^  ^  SqU3red  bracket  in  Equation  (73) 

fore,  the  electron  density  f  luctnnt- ^  the  present  tlow  conditions.  There- 
larger  (occasional!  y "  r^u  ^  |  In  XTE.**  "  P°  ^  ^  T 

the  jet  i yT:;.iubr:“at,on<  with  wven“b«  ^ 
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SECTION  V 


SUPERSONIC  WAKES  AT  ANCLE  OF  ATTACK 


5.1  INTRODUCTION 

Drastic  changes  are  introduced  in  the  wake  of  a  body  flying  at  supersonic 
speeds  when  the  body  is  inclined  at  some  incidence  (angle  of  attack)  with 
respect  to  the  flight  vector:  first,  the  increased  total  drag  generates 
an  intense  inviscid  wake,  i.e.,  large  gradients  appear  in  the  flow  surround¬ 
ing  the  inner,  or  turbulent  wake.  Secondly,  these  gradients  are  genera  1 1 \ 
asymmetric  about  the  wake  axis.  Thirdly,  changes  in  the  viscous  drag  change 
the  characteristics  of  the  viscous  wake.  Finally,  with  the  generation  of 
lifting  forces  vortical  flows  might  also  appear,  causing  added  complexities 
to  the  flow  field. 

Although  the  new  character  of  the  turbulent  wake  is  hard  to  assess  in  thes<j 
circumstances,  its  role  in  the  mixing  process  can  be  predicted  along  general 
lines.  Its  main  complexity  should  be  its  growth  in  and  interaction  with 
the  asymmetric  inviscid  wake.  The  flow  thus  consists  of  a  "shear  layer" 
contributed  to  by  the  inviscid  wake,  with  a  superposed  wake  flow.  It  is 
quite  clear  that  the  behavior  of  this  composite  flow  cannot  always  follow 
simple  wake  rules  of  lateral  growth,  for  example,  or  of  velocity  decay, 
since  shear  layers  and  wakes  are  generally  governed  bv  different  rules. 

Thus,  a  pure  shear  layer  supplies  a  constant  reference  velocity  scale  by 
definition,  while  the  corresponding  velocity  scale  of  the  wake  decreases 
with  distance.  In  limiting  cases,  the  competition  between  the  shear-laver- 
like  and  wake- like  behavior  can  be  foreseen:  for  small  angles  of  attack 
and  large  shear  drag,  the  wake-like  behavior  should  dominate  early,  while 
lor  inclination  angles  large  enough  to  cause  very  persistent  inviscid  asym¬ 
metries  the  flow  is  probably  shear- 1  aver- 1  ike  for  large  distances. 

The  present  experiment  was  designed  mainly  in  a  way  so  as  to  accentuate 
this  competition  between  the  inviscid  and  the  viscous  (turbulent)  wake. 
Calculations  of  the  supersonic  inviscid  flow  behind  a  two-dimensional 
slender  body  at  incidence,  verified  experimentally,  showed  a  considerable 
persistence  of  asymmetries  in  a  pattern  approaching  an  ideal  shear  layer. 
Thus,  the  problem  is  not  only  relevant,  but  was  done  in  a  way  capable  of 
fundamental  conclusions  on  the  basic  behavior  of  fluids. 

5.2  MODEL  AND  FACILITIES 

This  work  is  performed  in  the  supersonic  wind-tunnel  at  a  continuous  air 
flow  at  Mach  3.0,  a  stagnation  pressure  of  730  mi ,  'g  absolute  and  a  total 
temperature  of  39°C.  The  model  consisted  of  a  0  -inch  thick,  0.116-inch 
wide  stainless-steel  ribbon  stretched  across  the  .nel  test-section, 
identical  to  that  used  in  the  WED  experiment  (two-dimensional  adiabatic 
wake  at  zero  incidence).  A  major  modification  to  the  model  support  con¬ 
sisted  of  the  addition  of  a  fixed  steel  protractor  and  pointer  by  which 
incidence  angle  could  be  set  and  resolved  to  better  than  0.25  degree.  Twist- 
angle  resolution  is  about  0.1  degree.  The  protractor  is  shown  on  Figure  73. 
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FIGURE  73. 


FXPhR IH.XI’AL  ARRAXGICMICXT  FOR  ANGLIC- OF- ATTACK  STUDIKS 
HU  PROTRACTOR  tX  TIIK  TIC  ST  SICCTIOX. 


SHOWING 
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1  igure  74  sliows  the  arrangement  o  f  the  model  and  the  nomenclature  of  the 
low  phenomena  and  coordinate  axes. 

5.3  EXPERIMENT  DES1CN 

5.3.1  INVISCID  WAKE 

Detailed  calculations  of  the  inviscid  wake  of  ^ he  flat  plate  were  made  using 
the  method  of  two-dimensional  characteristics  at  an  angle  of  attack  of 
20  degrees;  for  different  angles,  the  flow  field  should  not  alter  qualita¬ 
tively  from  this  calculation,  provided  that  the  incidence  angle  for  shock 
de tacuaient  at  the  leading  edge  (about  34  degrees  at  this  Mech  number)  is 
not  exceeded.  In  this  calculation  the  upstream  boundary  conditions  on  the 
windward  side  were  t ho  oblique  shock  at  the  leading  edge;  for  the  leeward 
side,  the  trailing  edge  shock  as  modified  by  the  leading-edge  expansion 
fan.  Figure  75  sliows  graphically  the  result.  Attention  should  be  drawn 
to  the  two  nearly  parallel  streamlines  shown ing  on  either  side  of  the  wake 
axis.  It  turns  out  that  these  demark  abruptly  changing  conditions  in  the 
sense  of  flow  gradients  (not  in  the  sense  of  flow  discontinuities).  The 
net  result  is  that  between  each  of  these  two  streamlines  and  the  axis  there 
appears  a  uniform  "plateau"  in  flow  conditions;  the  two  adjacent  "plateau" 
constitute,  then,  a  shear-layer  situation  which  apparently  does  not  decay 
much,  either  ir  extent  or  intensity,  as  one  moves  downstream.  Typical  such 
plateaus  are  seen  in  cross-sectional  views  of  Figure  75  shown  in  Figure  76. 

The  calculations  outlined  above  demonstrate  that,  although  quant i La t ivel \ 
different  lor  the  axisvmmetric  hypersonic  wake,  an  essential  feature  of  the 
wake-a t- inc idence  problem  could  indeed  be  obtained  in  our  wind-tunnel.  As 
described  in  the  following  paragraphs,  the  shear-layer  feature  was  clearly 
evident  from  the  data  collected  so  far. 

5.3.2  CHOICE  OF  INOIDINCF  SKTTINCS 

To  study  the  eftoct  ot  incidence,  more  than  one  incidence  angle  should  be 
chosen,  if  possible;  a  judicious  choice  of  angles  is  needed  in  order  to 
illuminate  the  effect  without  undue  labor.  It  is  known  on  theoretical 
grounds  that  beyond  the  detachment  angle  (34  degrees)  the  inviscid  wake 
"freezes"  to  a  configuration  not  very  sensitive  to  angle  ot  incidence. 

Tests  verified  this  effect,  and  in  the  detached  case  it  was  decided  to  work 
only  with  a  =  90°  (plate  normal  to  the  stream)  so  that  the  blunt-body  tur¬ 
bulent  wake  could  at  least  be  probed. 

In  the  pre-detachment  range  of  awake  "pitot  maps",  such  as  shown  for 
CL  =  0°  on  Figure  30  of  Reference  1,  were  taken  at  a  =  10°,  15°,  20°  and  30°. 
Such  maps  are  easy  to  construct  and  are  invaluable  in  giving  an  overall 
picture  of  the  flow  field.  Since  the  shear- layer- l ike  behavior  was  present 
in  each  case,  the  angles  CL  =  10°  and  20°  were  chosen  for  further  study. 

In  designing  the  experiment,  the  farthest  downstream  distance  where  measure¬ 
ments  can  be  made  is  important  to  know.  In  the  wind-tunnel  this  distance 
is  shortened  when  the  model  is  pitched  because  of  the  earlier  reflection 
of  the  strengthened  shocks  from  the  tunnel  wall.  In  the  present  set-up  the 
farthest  downstream  distance  for  a=  0°,  10°,  20°  and  90°  is  7.8  inches, 
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7  inches,  b.2  inches  and  A. 7  inches,  respectively.  Note  that  lor  a-  0  , 
a  distance  ot'  7.8  inches  implies  a  distance  oi  about  1,800  viitual  body 
d i jraeters . 

5.4  DIAGNOSTIC  INSTRUMENTATION 

The  diagnostic  tools  used  for  this  experiment  were  identical  to  those 
employed  for  the  WED  (see  Reference  1).  The  0.00005  inch  diameter  P  L  10 
percent  Rh  hot-wire  used  for  this  work  managed  to  stay  intact  through  the 
mean-flow  and  turbulence  measurements  and  through  most  of  the  spectral 
measurements.  This  remarkable  wire  survived  45  calendar  days  oi  use, 
during  which  time  it  was  operated  for  about  70  hours  in  supersonic  1  low. 
Through  frequent  oven-  and  flow-calibrations  small  changes  in  its  charac¬ 
teristics  could  be  monitored  and  accounted  for. 


5 . 5  PROCEDURE 

This  experiment  was  designated  WEE  (Wake  Experiment  1 )  and  the  CX  =  10°, 
20°,  and  90°  sub-experiments  were  designated  WEI  -10,  WEI  -20,  and  WEI  — f  U  ’ » 
respectively.  As  done  with  the  WED  experiment,  the  distance  between  the 
model  mid-chord  and  the  end  (farthest  downstream)  point  was  marked  by  X- 
stations  spaced  by  0.500  inch  increments.  Continuous  analog  traverses 
were  taken  at  each  X-Station  with  the  pitot  and  static  probes  and  the  hot¬ 
wire  anemometer,  the  latter  operated  either  in  tl  e  mean-1  low  or  the  tuibu- 
1  on  t  mode.  From  these  data  the  mean  and  turbulent  properties  can  he 
extracted. 


5.6  RESUETS  OF  MEASUREMENTS 


As  already  mentioned,  the  actual  viscous  wake  measurements  were  preceded  by 
tracing  the  ir.viscid  wake  in  order  to  observe  the  overall  tlow  field  and  to 
examine  the'  ex  stance  ot  the'  " slie’ar- 1  aver"  el  Ic'ct.  figure1  7/  compares  a 
typical  pitot  traverse  with  the  pitot  pressure  txpected  I rom  the  character¬ 
istics  computation  shown  on  1  ig, vires  75  and  76  lor  a  -  20°.  It  is  obvious 
from  the  measurcjment  that  two  sho 'k  ■>  appear  on  each  side  cM  the  inviscid 


waki 


instead  of  one  as  given  b\  the  charac t e r i s t i c s  method. 


On  tin*  wind¬ 


ward  side,  a1 .  in  addition  to 
shock  is  thought  to  originate 
as  illustrated  tentatively  on 


neck" 


the  how  shock  at  the  leading  edge,  a  second 
in  the'  recompression  region  of  the  wake 
Figure  74.  In  the  leeward  side  a  shock 
appears,  of  course,  at  the  model  leading  edge  in  addition  to  the  expansion 
fan.*  Nevertheless  the  agreement  shown  on  figure  76  between  the  calculatec 
and  measured  pitot  pressure's  is  adequate;  both  demonstrate,  in  tae  t,  tlu 
"plateaus"  on  either  side  ot  the  turbulent  wake. 


A  second  preliminary  measurement  made  before1  detailed  studies  were  attempted 
concerned  the  existence  oL  turbulence.  Qualitative1  traverse's  with  the  hoL- 
wire  showed  that  angle  of  attack  did  not  inhibit  tlu  appearance  oi  turbu¬ 
lence  in  the  inner  wake.  Proceeding  to  the1  latter,  Figure  78  shows  pitot- 
pressure  traces  at  typical  X-S  tat  ions  tor  all  three1  experiments,  i  .  e .  , 

*Photographv  of  the  two-dimensional  model  itselt  is  impossible  because  ot 
obstructions.  The  identification  of  the  additional  two  shocks  is  lett  to 
reasonable  conjecture. 
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W I£ !•  -  10,  WEF-20,  and  WEI-  90.  A  trace  for  WEI)  (correspond i ng  to  WEF-0)  i -> 
also  shown.  The  object  of  this  figure  is  to  demonstrate  the  wake  symmetry 
which  understandably  obtains  for  the  symmetric  body  configurations,  namely 
lor  WED  and  WEF-90.  Note,  also,  the  expected  change  in  " f ree- s t ream"  pitot 
level  for  the  WEF-10  and  WEF-20  traces,  illustrating  the  " shear- layer" 
ef  fee  t . 

Typical  static  pressure  traverses  are  shown  on  Figure  79.  Although  these 
traces  appear  reasonably  smooth,  it  should  be  recalled  that  the  static- 
probe  resolution  is  quite  poor,  since  the  diameter  of  the  probe  was  never 
smaller  than  five  times  or  so  the  wake  thickness.  On  the  other  hand  there 
is  strong  belief  that  the  mixing  in  the  far  wake  occurs  at  constant  pressure, 
or  else  the  whole  flow  system  should  shift  progressively  in  the  direction  of 
decreasing  pressure.  In  fact,  this  probe  gave  seemingly  good  resolution  in 
regions  of  sensible  pressure  gradients,  such  as  at  or  near  shock  waves. 

Hot-wire  measurements  consist  ol  the  usual  three  sets:  one  set  consists  of 
live  curves  at  each  X-Station,  necessary  lo  complete  the  inputs  for  the 
mean  flow  measurements.  A  second  set,  with  15  curves  (each  at  different 
current)  traces  the  variation  of  mean  voltage  across  ihe  wake;  a  companion 
set  of  15  curves  plots  the  mean-square  of  the  wire  ac  (turbulent)  component. 
Here  one  should  note  that  the  turbulent  fluctuations  are  highest  on  the 
side  ot  the  wake  where  the  velocity  gradients  are  the  highest.  Spectral 
measurements  have  not  been  taken  as  yet. 
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SECTION  VI 

SUPERSONIC  WAKE  WITH  HEA1  TRANSFER 


6.1  INTRODUCTION 

There?  are  several  distinct  reasons  why  the  presence  ot  heat  transfer  is 
important  to  understanding  the  fluid  mechanics  of  the  high-speed  wakes. 

In  the  first  place,  it  is  well  kno\/n  from  stability  arguments  (Reference  7) 
that  the  location  of  wake  transition  to  turbulence  is  sensitive  to  heat 
exchanged  between  the  bodv  and  the  flow.  Secondly,  the  extent  to  which 
heat  transfer  changes  t lie  self-preservation  behavior  is  not  well  understood. 
Third,  the  fluctuations  in  the  wake  are  much  more  complex  with  heat- t rans for 
tor  example,  the  simple  relationship  between  velocity  and  density  fluctua¬ 
tions  valid  tor  the  isoenergctic  woke  does  not  apple. 

In  the  present  series  ot  experiments,  it  is  planned  to  study  the  mean, 
intermittent  and  turbulent  flow  component  behind  an  axisymnetric  cooled 
model  and  a  two-dimensional  heated  model.  'llu  latter  experiment  has  been 
initiated  and  most  of  the  tests  have  actually  been  completed.  It  will  be 
described  in  this  section. 

6.d  MODEL  AND  FACILITIES 

This  experiment  is  performed  in  the  same  wind-tunnel,  at  the  same  conditions 
with  the  same  model  and  concurrently  as  the  angl e-of-j t tack  experiment  dis¬ 
cussed  in  Section  V.  Provisions  to  the  model  support,  additional  to  those 
described  in  the  latter  section  as  necessary  tor  the  angl  e-of-a  t  tack  measure¬ 
ment,  were  incorporated  specifically  for  the  heat- transfer  measurement. 

The  steel  ribbon  is  heated  electrically  by  a  Harrison  Laboratories  Model 
No.  81oA  regulated  dc  power  supplv.  Load  regulation  is  held  to  bi  low  0.05 
percent,  and  power-line  ripple,  already  low,  is  made  et tec  lively  impercep¬ 
tible  because  of  high  thermal  inertia  ot  the  heated  ribbon.  The  ribbon 
current  can  be  adjusted  from  0  to  25  amperes  dc  and  is  continuously  dis¬ 
played  on  a  Speedomax-W  strip-chart  recorder  running  at  a  speed  ot  1  inch 
of  paper  per  minute.  Aided  by  the  chart  resolution  ind  a  fine  variable 
resistor  added  externally  to  the  circuit,  the  operator  can  thus  readjust 
and  control  the  ribbon  current  to  about  0.015  percent  ot  the  desired  value, 
an  error  lying  much  below  the  current  change  needed  lor  perceptible  changes 
in  the  heated  wake  characteristics.  The  power  settled  upon  for  this 
measurement  will  be  given  farther  below. 

Expansion  ot  the  model  by  thermal  rdaxition  necessitated  a  scheme  by  which 
the  ribbon  tension  could  be  adjusted  automatically.  This  was  done  by 
spring- loading  one  end  of  the  ribbon  so  that  the  model  could  expand  when 
heated  and  retract  to  a  stretched  position  when  the  electrical  power  was 
Lu  rned  o  t  t . 
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b.J  ADJUSTMENT  OF  Ol'EKA'l  INI  COM)  11 IONS 

lo  bring  about  easily  observable  lien  t-  t  rans  1  er  etucts,  it  was  desired  to 
dissipate  the  highest  possible  amount  oi  power  in  the  ribbon.  Several 
other  conditions  should  he  satisfied  s  imu  1  t.ineous  1  v ;  the  ribbon  should  not 
/oaken  structural ly  or  deform,  the  turbulence  should  not  Ik  inhibited  to 
th0  p0lnt  whorc  the  wnko  became  wholly  laminar,  and  the  ribbon  temperature 
non-uni  t  ormi  tv  because  of  heat  conducrion  into  tin  supports  should  he  kept 
immal.  There  was  furthermore  some  fear  that  the  increased  wake  tempera¬ 
tures  could  restrict  the  operating  temperature  of  the  hot-wire  or  even 
destroy  it  when  both  the  model  and  the  hot-wire  operating  temperatures  were 

ugh.  As  it  developed,  all  these  criteria  were  satisfied  once  the  ribbon 
maximum  power  was  set. 

ibis  was  accomplished  by  gradually  increasing  the  ribbon  current  until  the 
mo  .  glowed  faintl\  in  the  supersonic  flow.  The  current  was  then  reduced 
unUl  the  glow  disappeared.  At  this  point  the  ribbon  current  was  set  at 
1  '.  *'  amperes  giving  a  voltage  drop  across  it  ot  ..  jV*  volts  dc  for  a 
total  power  dissipation  ot  W  =  86.8  watts.  It  is  interesting  to  compare 
this  power  to  the  "Kinetic”  power  of  the  flow  P  =  (drag)  >.  (How  velocity) 
Uiliging  the  viscous  drag  coefficient  of  the  ribbon  from  the  zi  ro-incideiu’t 
aU)  experiment  (Section  II)  a  ratio  of  U’/P  =  l.j  js  obtained. 

Alter  the  necessary  tests  for  finding  tin  ini,  zero-twist  and  zero- iiu  idenc  < 
angles  (see  Keterence  1),  the  hot-wire  was  enp loved  qua  1 i tet i vf 1 v  for 
s  tidying  the  transition  to  turbulence  with  and  tU  transverse  (span-  is,-) 
hi  oi  the  heating  process.  '1  rans  it  ion  was  cK  tested  hv  the  usual 

ice!  ot  noting  the  point  oi  maximum  turbuhna  inlensitv  along  tat  cenur 
’  me  (Keterence  28).  Figure  80  shows  tin  remarkable  diiureiice’s  in  transi¬ 
tion  distance  observed  with  heating*  this  distance,  located  about  1  <>0 
virtual  diameters  (momentum  thicknesses)  behind  tin  adiabatic  bode,  mol 
to  about  800  diameters  when  beating  was  applied.  This  phenomenon'  is  con¬ 
sistent  with  the  predictions  oi  stability  theory  by  which  tin  laminar  wake 
IS  stabilized  it  the  model  or  vehicle  producing  ii  is  heated  (Reference  7). 
Oik  result  o!  this  t rans i t i on  del  ay  was  that  the  available  range  in  the 
axiai  direction  over  which  diagnostics  could  he  mack  was  only  about  hall 
tin  range  available  for  wake  studies  with  the  adiabatic  body.  in  exchange, 
however,  one  obtains  a  long  transitional  wake  iron  which  information  impor¬ 
tant  to  the  transition  process  could  be  extracted.  lor  this  reason  probings 
began  again,  as  in  WFI),  shortly  downstream  of  the  bod\  rather  than  at  the 
transition  distance. 

Visual  observation  of  the  ribbon  following  its  overheating  to  the  glowing 
point  showed  thru  its  rear  hall  (from  about  mid-chord  to  the  trailing  edge) 
had  been  discolored  in  the  way  typical  ot  metals  after  exposure  to  hi  h 
temperatures.  That  the  ribbon  temperature  increases  along  its  chord  toward 
the  trailing  edge  is  fully  expected  because  of  the  more  efficient  heat 
transier  tor  the  thinner  boundary  layer.  The  discolored  region  extended 
uniformly  along  the  span,  but  it  was  observed  to  stop  abruptly  about  on. 


••Henceforth  ’beating”  refers  to  the  86.8  watts  settled  upon  as  note, 
previ ous 1 v. 
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FIGURE  80.  EFFECT  OF  BODY  HEATING  ON  TRANSITION  DISTANCE  (WEG) 


ribbon  chord  before  it  reached  the  tunnel  sidewalls.  The  ends  of  the 
ribbon  were  therefore  cooler  since,  of  course,  heat  conduction  to  the 
soUd  supports  increased  at  those  points. 

Hie  observations  described  above  are  crucial  to  the  span-wise  uniformity 
ot  heat  transfer  from  the  ribbon  to  its  wake,  that  is,  the  question  of 
true  two-dimensionality  of  the  experiment.  This  question  was  further  exa¬ 
mined  by  taking  hot-wire  profile  of  the  local  total  wake  temperature  at 
fixed  distances  behind  the  model  and  difference  spanwise  stations.  Within 
an  inch  of  the  spanwise  direction  around  the  center-span  position  (covering 
about  J5  percent  of  the  total  span),  the  wake  axis  temperature  varied  o\ 
about  10  percent;  at  least  part  of  this  variation  is  attributable  to  slight 

flow  non-uniformities.  The  two-dimensionality  was  thus  judged  to  be 
adequa  te . 


6.4  PROCEDURES  AND  INSTRUMENTATION 


Hie  instrumentation  used  for  this  work  was  identical  to  that  used  for  the 
WED  and  WEF  experiments.  Although  the  ribbon  itself  was  heated  highly, 
turbulent  diffusion  reduced  the  local  stagnation  temperatures  to  'he  point 
where  no  question  of  damage  to  the  probes  by  high  temperatures  ever  arose. 
In  tact  the  0.00005  inch  diameter  wire  was  operated  in  the  same  full  range 
ot  heating  currents  (maximum  of  about  7.5  mi  1 1 iamperes)  without  ever  fail- 
ing  b\  causes  attributable  to  high  wake  temperatures. 


Ihc  present  experiment,  done  with  the  heated  two-dimensional  bodv  at  zero 
incidence,  is  designated  WEG  (Wake  Experiment  G) .  The  measurements  were 
performed  at  each  ot  23  axial  positions  along  the  wake.  (X  STATIONS  1 
through  23)  beginning  at  0.658  inch  from  the  ribbon  mid-chord  point.  At 
each  station  one  pitot-pressure  transverse,  one  static-probe  transverse, 
and  five  hot-wire  traverses  were  performed  to  collect  data  sufficient  to 
(.escribe  fully  the  mean  1 1  ow  field.  Each  wire  traverse  was  ot  course  per- 
t ormed  at  different  heating  currents.  For  the  turbulence,  15  traces  ot 
mean  wire  voltage  and  rms  voltage  were  collected  at  each  station.  lor  the 
spectra,  ten  distinct  radial  locations  were  chosen  at  whiJi  spectra  were 
recorded  at  each  ot  six  different  heating  currents.  Since  the  acquisition 
an  reduction  ot  the  spectral  data  is  extremely  time-consuming,  data  were 
taken  only  at  alternate  X  STATIONS  (the  odd-numbered  stations'  1  3  5 

23).  Thus  a  total  of  10  x  6  x  12  =  720  spectra  were  recorded. 


It  is  perhaps  advisable  to  reiterate  here  the  data  reduction  process  through 
which  the  WEG  data  are  now  progressing.  The  mean-flow  data  have  been 
digitized  and  put  into  the  WEb-V  and  WEB-V1U  computer  programs  which  have 
for  some  time  been  prepared  for  the  l’hilco-2000  digital  computer.  The 
spectral  data  are  first  used  into  the  WEIJ-II  program  from  which  one  extracts 
the  so-called  error  ratios,  that  is,  a  correction  factor  necessary  to  use 
in  correcting  the  system  frequency  response.  The  output  of  this  program 
and  the  turbulence  data  are  then  used  into  the  WEB-311  and  WEB-IV  programs 
to  deduce  the  turbulence  properties  and  further  into  the  WEI1-V11  program 
from  which  turbulence  properties  in  scaled  form  are  obtained. 
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SECTION  VII 


PROJECT  LOGISTICS 


7  .  1  PERSONNE1 

This  work  is  performed  for  the  Space  and  Re-entry  Systems  Division  under 
the  general  supervision  of  John  Knudsen.  The  technical  work  is  done  at 
the  Fluid  Mechanics  Department  of  the  Aeronutronic  Advanced  Development 
Operation  by  A.  Demetriades,  Supervisor  of  tne  Experimental  Aerodynamics 
Section,  assisted  by  Ernest  L.  Doughman,  Research  and  Development  Engineer 
and  Lee  Von  Seggern,  Engineer.  Computation  services  are  rendered  by  the 
SRS  Computation  Department  under  N.  llabibe,  Supervisor. 

7.2  DEVELOPMENT  OF  FACILITIES,  INSTRUMENTS  AND  TECHNIQUE 

7.2.1  FLOW  FACILITIES 

Tile  work  described  in  this  report  is  performed  in  the  supersonic  wind  tunnel 
and  plasma  wind  tunnel  constructed  and  maintained  by  Philco-Ford.  In  the 
present  period,  improvements  have  been  made  by  the  Corporation  to  these 
facilities  for  the  purpose  of  enhancing  the  efficiency  of  the  work  and  the 
quality  of  data.  An  automatic  current  control  has  been  installed  in  the 
plasma  wind  tunnel  for  the  purpose  of  avoiding  excursions  in  the  arc  current 
of  t'ne  plasma  generator  and  the  resulting  changes  in  the  turbulent  jet 
characteristics .  An  improved  plasma  generator  has  also  been  purchased  with 
which  improvements  in  the  precision  and  repeatability  of  the  measurements 
will  be  sought. 

7.2.2  INSTRUMENTATION 

The  following  specialized  instrumentation  has  been  added  by  the  c oar' any  in 
the  current  period  for  the  purpose  of  facilitating  and  improving  current 
and  projected  work: 

a.  Heat -Flux  Probe.  A  heat  flux  probe  has  been  constructed  for  cross¬ 
checking  mean- flow  data  acquired  in  the  plasma  jet  by  measuring  the  local 
heat  flux.  The  sensor  lias  a  flat  circul  face  0.25 -inch  diameter  and  a 
calibrated  range  of  500  Btu/ft-/sec.  Its  output  is  1  nv  per  50  Btu/ft  /sec. 
This  is  shown  in  Figure  81. 

b.  Ultra-Small  Hot-Wire  Probes.  The  smallest  hot-wire  probes  used 
in  the  supersonic  wind  tunnel  to  date  have  been  0.01- inch  long  and 
0.00005- inch  diameter,  and  have  a  limiting  (zero-current)  time  constant  of 
0.100  msec.  Although  the  frequency  response  of  these  wires  is  greatly 
augmented  by  electronic  and  computational  methods,  there  is  continuing  need, 
in  this  work,  for  higher  response  instruments.  Attempts  were  made  in  the 
current  period  to  mount  wires  oi  0.00002  and  0.00001  inch  diameter  (20  and 
10  microinches,  respectively)  which  -potentially  have  time  constants  of  25 
and  10  microsec  uids.  Reportedly,  such  wires  have  been  mounted  elsewhere  ui  1  \ 
rarely  and  have  never  been  before  used  in  supersonic  flow.  So  far, 
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FIGURK  81. 
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0  00002- inch  diameter  Pt  10  percent  rh  wires  have  been  successfully  mounted 
in  lengths  of  0.004  inch,  annealed  and  calibrated  in  the  controlled  ovei  . 
They  have  also  been  immersed  in  the  supersonic  flow  without  breakage  and 
there  heat- transfer  characteristics  in  that  flow  have  been  measured, 
are  pictured  in  Figures  82  and  83. 

c.  Wire  Conditioning  Circuit.  A  hot-wire  conditioning  circuit  for 
mounting ,  annealing  and  calibrating  hot-wire  probes  has  been  designed  and 
built  by  Phi lco-Ford  personnel.  It  has  greatly  facilitate  tie  io  w 

work  reported  herein. 

^  iHah-Frenuencv  Amplifier.  The  need  to  cover  the  frequency  range 
of  turbulence  extant  in  presented  projected  experiments  without  overtaxing 
present  techniques  of  electronic  and  computer  compensation,  has  beei 
frustrated  by  the  lack  of  commercially  available  hot-wire  amplifiers  f  c 
to  2  MHz.  A  new  amplifier  has  therefore  been  designed  and  is  now  being 
built  at  the  Company.  Frequency  response  is  designed  to  2  MHz  at  a  gam  o 
1,000.  Compensation  to  10  microseconds  will  be  available,  as  we 
innovations  such  as  pushbutton  wire  current  selection  Iwo  such  nmpiif 
are  being  built  since  two-point  correlation  measurements  in  the  wake  aie 

being  planned. 

e  Correlation  Probes.  For  two-point  correlation  work,  two  novel 
probes  are  being  readied.  For  the  supersonic  wind  tunnel  a  two;dcgree-o 
freedom  (five  degree  if  the  translation  of  the  entire  probe  ^ d^cc-of- 
two-wire  probe  will  be  used*  for  the  plasma  wind  tunnel  a  si..-degiee 
freedom,  two-wire  (or  langmuir)  probe  is  being  adapted  from  an  en  lie 
design.  These  probes  pictured  in  Figures  84  and  85  wi  e  ac  u 
remotely  so  that  the  operator  will  be  able  to  change  the  spacing  of  the 
sensors  rapidly  and  at  will  along  any  of  two  (or  three)  Cartesian  a..e  . 

f.  Processing  Electronics.  Other  electronic  circuits  now  planned 
include  a  constant- temperature  anemometer  amplifier  and  a  circuit  oi 
performing  one  or  more  of  the  following  operations  on  two  separate  signals, 
addition,  subtraction,  multiplication,  division,  averaging,  an  m  egia  i 


7.2.3 


ADDITIONS  TO  COMPUTATIONAL  FACILITIES 


In  addition  to  the  regular  computing  facilities  available  from  the 
Aeronutronic  Computing  Services  Center,  Time-Si  aring  services  are  now 
conveniently  available  for  this  work.  A  terminal  to  the  Los  Angeles 
General  Electric  235  System  is  available  on  full-time  basis.  The  mo 
common  language  is  BASIC  and  FORTRAN  II.  in  adtUion  to  L  u  -  ‘  \PSE1  1) 

program  shown  in  Appendix  B,  two  BASIC  language  programs  (LIM  IT  and  M 
have  been  written  for  calibrating  the  hot-wire  in  the  controlled  oven  and 
in  the  supersonic  flow.  These  are  shown  m  Figures  8b  r-ul  o/. 
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AND  "LARGE"  HOT-WIRE  PROBES 


SYMBOL  WIRE  NO.  DIAM. (INCH) 


143 


HEAT  TRANSFER  CHARACTERISTICS  OF  TYPICAL  0.00002-INCH  WIRE 
OF  SAME  ASPECT  RATIO  AS  0.00005- INCH  WIRE 


FIGURE  83.  TWO-WIRE  CORRELATION  PROBE  FOR  WORK  IN  THE  PLASMA  TUNNEL 
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7.2.4  THE  NEW  DATA  ACQUISITION  AND  A/D  CONVERSION  SYSTEM 

In  the  present  period  a  review  was  made  of  needs  to  improve  the  flow  of 
work  from  the  data  acquisition  to  the  data  analysis  stage.  Four  major 
needs  were  delineated: 


(1) 

Accelerate  the  analog-to-digi tal 

conversion 

proces  s . 

(2) 

Accelerate  certain  steps  in  data 
such  as  recording  the  spectra. 

acquisition 

(3) 

Store  the  raw  signals  of  a  complete  experiment 

(A) 

Enable  certain  measurements  such 
point  correlations. 

as  two- 

in-  ’  u  I1"  f°5/the//D  conversion  scheme  were  (1)  continuous  monitor¬ 
's1"  a  °8  V ^  dl8ital  staSe>  and  <2>  compatibility  with  the 
ihilco-.OOO  digital  computer.  Candidates  for  conversion  were,  first 

vnif1Varat!  Si8"a]?  mean-square  wire  voltage  versus  probe  position 

ltage)  and  secondly,  the  high-frequency  time-dependent  raw  hot-wire 
signal  containing  frequencies  up  to  2  to  3  MHz.  The  conversion  of  each  dc 
signal  is  accomplished  by  channeling  the  analog  voltage  through  a  digital 

IBM  mTII  ^!eW^t“PaSkard  M°del  N°-  344°A)  t0  3  COU'^er  a"d  then  toban 
IBM  Model  No  526  card-puncher  which  is  continuously  monitored  by  a  digital 

recorder  (Hewlett-Packard  Model  ;,o.  5050A) .  The  output  o£  the  litter  alio 

p  lL’°,OPtl°n,?f  ra,'lnu:l1  keypunch  if  desired.  Finally,  a  Mosely  Model 
No.  F-JB  Line  Follower  System  will  be  available  for  non-real-tine  A/D 
conversion  of  data  already  accumulated  in  analog  form. 

Spectra  have  been  hitherto  recorded  in  analog  form  at  five-minute  intervals, 
n  ie  supersonic  tunnel  as  many  as  1000  spectra  per  experiment  are  needed 
resulting  in  a  great  deal  of  elfort.  In  tlu  new  system,  the  sweep  with  ’ 
frequency  is  done  electronically  in  real  time,  so  that  an  acceleration  of 
the  spectral  measurements  by  a  factor  of  3  to  5  is  thought  possible. 

A  Honeywell  Model  No.  7«00  seven-channel  magnetic  tape  recorder  has  been 
acquired  ma.nly  in  the  purpose  >c  staling  the  raw  signals  tor  future  pr  .- 
cessing  as  the  need  arises  and  as  tedn  iqcies  become  available.  A  second 
important  use  ot  the  playback  will  be  to  "splice-on”  the  tape  electronically 
t  r  leasurement  ol  the  eddy  properties.  A  third  1  unction  will  be  to  bring  ' 
the  signals  within  the  frequency  range  >,:  processing  electronics  having 
^tate-ol-tlu-ait  bandwidth  limitations.  For  example,  the  state-of-the-art 

blllw^ho  2  00f)'akll  ’  "’"dCla-  tr°SS*eorr*  ,ators  r  urn  nt  I  y  is  23(1  kite,  muc  h 
below  the  2,000  kite,  range  ot  turbulence.  By  recording  the  latter  signal  at 

ips  an  p  ayinc  back  at  15  ips  the  Signal  comes  within  the  ^50  kHz 
range. 
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7.3  INFORMATION  EXCHANGE 


Contact  with  other  ARPA  contractors  has  been  maintained  and  accelerated  in 
the  current  period  through  personal  exchanges  (visits,  telephone  conferences 
etc.)  and  technical  reports  issued  by  Philco-Ford  and  distributed  via  SAMSO 
channels.  Especially  useful  have  been  the  exchanges  with  the  following: 

(1)  Cal  Tech  Hypersonic  Group  (Lees,  Kubota,  Behrens). 

(2)  University  of  Southern  California  (Lauter) . 

(3)  Jet  Propulsion  Laboratory  (Kendall)  . 

(4)  University  of  California,  La  Jolla  (Lin, 

Gibson) . 

(5)  University  of  Colorado  (Sanborn) . 

(6)  Bell  Telephone  Laboratories  (Granatste in) . 

(7)  TRW  Systems  (Webb,  Lewis). 

(8)  RCA  Montreal  (Johnson) . 

(9)  Stanford  Research  Institute  (Cuthart) . 

(10)  Aerospace/El  Segundo  (Ross). 

(11)  Aerospace/SBO  (Renau) . 

The  three  Technical  Reports  issued  in  the  current  period  are: 

(1)  "Turbulent  Front  Structure  of  An  Axif.ymmetric 
Compressible  Wake",  Philco-Ford  Publication 
No.  UG-4259 ,  SAMSO  TR  68-44,  15  November  1967. 

(2)  "Mean  Flow  Measurements  In  A  Self-Preserving 
Turbulent  Plasma.  Jet",  SAMSO  TR  68-166, 

February  1968. 

(3)  "Turbulent  Mean-Flow  Measurements  In  A  Two- 
Dimensional  Compressible  Wake",  SAMSO  TR 
68-369,  July  1~58. 

Publication  in  the  technical  literature  of  selected  portion  of  this  work, 
as  approved  by  SAMSO,  has  stimulated  valuable  feedback  from  the  technical 
community.  In  the  current  period  publication  activity  was  as  follows: 

(1)  "Mean  Flow  Measurements  in  An  Axisymmetric 

Compressible  Turbulent  Wake",  AIAA  I.,  Vol.  6, 

No.  3,  p.  432,  March  1968. 
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(2)  1  Turbulence  Measurements  la  An  Axisymmetric 
Compressible  Wake",  Physics  of  Fluids, 

Vol.  11,  No.  9,  p.  1841,  September  1968. 

(3)  Turbulent  Front  Structure  of  An  Axisymmetric 
Compressible  Wake",  to  be  published  in  the 
Journal  of  Fluid  Mechanics,  1968. 

(4)  "The  Mean  and  Intermittent  Flow  of  A  Self- 
Preserving  Plasma  Jet",  accepted  for  publication 
in  the  AIAA  Journal. 

(5)  "Turbulent  Mean-Flow  Measurements  In  A  Two- 
Dimensional  Supersonic  Wake",  accepted  for 
publication  in  the  Physics  of  Fluids. 

(6)  "Structure  of  Compressible  Turbulent  Wakes", 
Proceedings  of  the  25th  AMRAC  Meeting,  Washington, 
D.C.,  24-26  April  1968,  p.  327. 
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APPENDIX  A 


THE  JEA-II  COMPUTER  PRf  '.RAM 


The  purpose  of  the  JEA-II  program  was  to  reduce  the  raw  data  of  the  electron 
density  spectra  measurements  taken  in  the  turbulent  plasma  jet.  The 
program  also  calculates  the  total  mean  square  values  and  the  longitudinal 
integral  scale.  Inputs  to  this  program  were  the  wave  analyzer  center 
frequency,  F,  signal  level,  Vs,  noise  level  Vn  per  200  Hz  bandwidth,  the 
mean  current  E^an,  directly  from  the  Langmuir  probe  and  the  mean  velocity 
U  taken  from  the  JEA-I  computer  program.  Input  consisted  of  data  at  zero 
Hz,  500  Hz,  1000  Hz,  and  then  in  steps  of  1000  Hz  up  to  frequencies  where 
the  signal  and  noise  levels  were  the  same.  The  highest  frequency  was 
2b  kHz,  Therefore  the  input  consisted  of  less  than  30  input  frequencies 
per  spectrum. 

The  program  output  listed  the  following  data:  the  probe  position,  x  in 
inches  from  the  jet  nozzle,  Y  in  volts  (1  volt  =  1/4  inch),  the  total  mean 
square  value  equal  to  the  area  under  each  spectrum  curve  E  INT,  the  ratio 
of  the  mean  square  spectal  density  at  zero  frequency  to  the  total  mean 
square  value  under  the  spectra  curve  E(0)2/E2,  and  the  longitudinal  integral 
scale  LAMDAF.  Listed  as  a  function  of  frequency  are  the  following:  rms 
signal  value  per  200  Hz  passband  E(F)S,  rms  noise  value  per  200  Hz  passband 
E(F)N,  mean  square  value  per  200  Hz  passband  of  the  true  signal  (signal- 
noise)  E(F)2,  the  ratio  of  the  true  rms  per  200  Hz  to  the  total  mean  value 
for  all  frequencies,  Al,  the  ratio  of  the  mean  square  per  cycle  to  the  total 
mean  square  for  all  frequencies  A2 ,  ratio  of  the  mean  square  per  cycle  to 
the  normalized  mean  square  at  zero  frequency  A3,  and  finally  the  normalized 
wave  number,  FN.  The  program  computational  procedure  appears  o-»  Figure  A-l, 
its  listing  on  Figure  A-2  and  a  typical  output  page  on  Figure  A-3. 
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FIGURE  A- 1 .  PROCEDURAL  DIAGRAM  OF  THE  JEA  PROGRAM 
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APPENDIX  B 


THE  SPECTR  COMPUTER  PROGRAM 


Tie  SPECTR  computer  program  was  written  for  the  GE  Time  Sharing  computer 
to  reduce  the  hot  wire  anemometer  temperature  spectra  data  taken  in  the 
turbulent  plasma  jet.  Inputs  to  the  program  as  a  function  of  every  333.3  Hz 
were  signal,  S,  and  noise,  N,  levels.  The  program  was  limited  to  41  sets. 
FORTRAN  II  language  was  used.  Other  input  for  each  spectra  consisted  of 
mean  velocity,  u,  from  the  JEA-1  program,  and  hot  wire  time  constant,  TW. 

A  number  of  other  inputs  were  initially  integrated  into  the  program  such 
as  wave  analyzer  bandwidth,  compensated  time  constant  setting,  amplifier 
gain,  and  the  equations  describing  the  hot  wire  response  curve  and  tiie 
characteristic  roll-off  curve  of  the  hot  wire  amplifier,  as  well  as  the 

conversion  factor  for  changing  inch  readings  of  data  measurements  into 
millivolts . 

The  program  listed  the  following  outputs  as  a  function  of  frequency: 

The  signal  and  noise  D|  levels;  true  signal  H,  (signal-noise);  correction 
factor,  OT,  to  adjust  for  drop-off  of  hot  wire  response  and  amplifier  with 
frequency;  corrected  mean  square  of  the  true  signal,  E2 ;  normalized  corrected 
mean  square  of  the  true  signal  DN;  and  finally  the  nondimensiona lized  wave 
number  FL.  Also  included  is  the  total  mean  square  value  under  the  spectral 
curve  before  using  the  correction  factor  II,  and  after  using  the  correction 
factor  12;  the  ratio  of  the  two  is  given  as  J;  and  the  longitudinal 
temperature  integral  scale  is  given  as  LT.  The  computational  procedure 
of  this  program  appears  on  Figure  B-l,  the  program  appears  on  Figure  B-2, 
and  a  typical  output  (at  a  point  in  the  jet)  on  Figure  B-3. 
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FIGl’KL  15-1.  mxiKin'RAL  DIAGRAM  OK  TIIII  Sl’KCTR  PROGRAM 
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